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Conversion  Factors,  Non-SI  to 
SI  Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units  as 
follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.01745329 

radians 

feet 

0.3048 

meters 

inches 

25.4 

millimeters 

kips  (force) 

4.448222 

kilonewtons 

kips  (force)  per  square  inch 

6.894757 

megapascals 

miles  (U.S.  statute) 

1.609347 

kilometers 

pounds  (force)  per  square  inch 

6.894757 

kilopascals 

square  inches 

6.4516 

square  centimeters 

tons  (short,  2000  lb) 

907.1847 

kilograms 
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1  Introduction 


1.1  Background 

One  of  the  most  frequent  loads  applied  to  the  locks  of  the  inland  waterway 
system  is  the  impact  made  by  a  barge  flotilla  as  it  transits  the  lock.  Consequently, 
this  load  case  represents  one  of  the  primary  design  loads  considered  by  the  U.S.  ' 
Army  Corps  of  Engineers  for  lock  walls,  approach  walls,  guide  walls,  and  guard 
walls.  The  primary  focus  of  engineers  performing  these  impact  computations  has 
been  the  lock  approaches,  where  the  worst-case  loads  are  likely  to  occur. 

In  1993,  the  Headquarters,  U.S.  Army  Corps  of  Engineers,  issued  the  first 
formal  Corps-wide  analysis  procedure,  providing  guidance  for  analyzing  the 
effects  of  barge  impact  loading  on  navigation  structures  in  the  form  of  an 
engineer  technical  letter  designated  as  ETL  1 1 10-2-338.  This  ETL  gives  the  basic 
equations  of  an  engineering  procedure  for  the  collision  of  a  barge  flotilla  with  a 
rigid  structure.  According  to  the  ETL  1 1 10-2-338  engineering  procedure,  the 
magnitude  of  the  impact  forces  generated  by  a  particular  collision  event  is 
dependent  on  the  mass  (including  hydrodynamic  added  mass  of  the  barge 
flotilla),  the  approach  velocity,  the  approach  angle,  the  barge  flotilla  moment  of 
inertia,  damage  sustained  by  the  barge  structure,  and  friction  between  the  barge 
and  the  wall. 

A  major  distinction  between  this  procedure  and  the  traditional  Navy  method 
for  determining  berthing  forces  is  the  estimation  of  collision  energy  dissipated  in 
deformation  of  the  barge  structure  and  transferred  to  the  rotation  of  the  barge 
flotilla.  The  analytical  method  uses  the  structural  interaction  mechanism  of 
Minorsky  (see  ETL  1 1 10-2-338),  which  provides  an  empirical  relationship 
between  the  (nonrecoverable)  hull  deformation  and  the  energy  absorbed  in  a 
collision.  The  relationship  between  kinetic  energy  lost  in  a  collision  and  the 
volume  of  in-plane  (barge)  material  damaged  is  used  to  determine  impact  force 
as  a  relationship  to  instantaneous  contact  area  of  damaged  structure. 

Minorsky  used  the  conservation  laws  of  momentum  and  energy  and  the 
principles  of  rigid  body  mechanics  to  estimate  the  kinetic  energy  lost  during  a 
collision  between  two  vessels.  He  then  calculated  a  resistance  factor  that  is 
essentially  the  volume  of  material  damaged  in  the  bow  of  the  striking  ship  and  in 
the  side  of  the  struck  ship.  Minorsky  reasoned  that  the  principal  resistance  to 
collision  penetration  is  provided  by  deep  structure  that  suffers  in-plane  damage. 
For  the  case  of  a  barge  striking  a  fixed  wall,  the  main  deck,  the  bottom  plate,  the 
head  log,  and  the  transverse  frames  would  offer  resistance  to  damage.  Minorsky 
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selected  and  analyzed  26  actual  ship  collisions  and  correlated  the  energy 
absorbed  in  the  collision  with  the  Minorsky  resistance  factor.  Using  the 
equivalency  between  energy  absorbed  and  the  work  performed  in  deforming  the 
structure,  a  constant  described  as  the  force  per  unit  of  damaged  surface  area  was 
defined  (=  13.7  ksi).1 

The  Minorsky  structural  interaction  mechanism  is  a  constant  pressure  process 
operating  with  a  pressure  of  13.7  ksi  acting  over  the  instantaneous  face  area  of 
the  damaged  element.  This  allows  for  the  definition  of  an  equivalent,  linear 
spring  constant  representing  the  crushing  of  the  barge  structure  in  the  ETL  1110- 
2-338  analytical  formulation.  It  is  important  to  note  that  the  entire  structural 
interaction  mechanism  is  modeled  as  a  linear  spring  in  the  direction  of  collision 
corresponding  to  the  energy  absorption  in  the  crushed  barge  structure.  The 
formulation  becomes  one  of  an  initial  value  problem  for  barge  flotilla  collision 
with  a  rigid  wall,  representing  a  lock  wall  in  this  case,  and  leads  to  the  solution 
given  in  ETL  1 1 10-2-338. 

Two  significant  concerns  have  been  raised  since  the  ETL  1 110-2-338 
procedure  was  released.  First,  a  key  aspect  of  the  ETL  1 1 10-2-338  engineering 
formulation  is  computation  of  collision  energy  dissipated  in  nonrecoverable, 
plastic  hull  deformation  of  (i.e.,  damage  to)  the  comer  of  the  barge  where  impact 
with  the  wall  occurs.  However,  the  majority  of  the  impacts  made  by  barge 
flotillas  transiting  Corps  locks  do  not  result  in  damage  to  the  barge  structure  or  to 
the  walls.  Second,  several  engineers  who  have  used  the  ETL  1 110-2-338 
engineering  procedure  have  questioned  the  accuracy  of  the  computed  results. 

To  investigate  these  issues  as  well  as  to  provide  a  basis  for  the  development 
of  an  improved  numerical  impact  model,  a  fully  instrumented,  full-scale  impact 
experiment  was  devised  to  directly  measure  the  impact  forces.  This  report 
addresses  the  interpretation  of  the  resulting  December  1998  full-scale,  low- 
velocity,  controlled  impact,  barge  flotilla  impact  experiments  conducted  at  the 
decommissioned  Gallipolis  Lock  at  Robert  C.  Byrd  Lock  and  Dam.  A 
comparison  between  measured  impact  forces  for  these  full-scale,  barge  flotilla 
impact  tests  and  corresponding  computations  made  using  the  ETL  1 1 10-2-338 
engineering  procedure  is  included  in  this  report. 

In  addition,  a  easy  to  use  empirical  correlation  is  derived  that  reports  the 
maximum  impact  force  (normal  to  the  wall)  as  a  function  of  the  linear 
momentum  normal  to  the  wall  (immediately  prior  to  impact),  using  the  results 
from  the  impact  forces  measured  during  these  full-scale  impact  experiments.  The 
authors  envision  that  this  new  empirical  correlation  will  be  used  for  impacts  that 
do  not  involve  damage  during  impact  to  either  the  comer  barge  of  a  barge  flotilla 
or  to  the  wall.  An  alternate  empirical  correlation  is  given  for  the  maximum 
impact  force  (normal  to  the  wall)  as  a  function  of  the  kinetic  energy  normal  to  the 
wall  (immediately  prior  to  impact). 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  xi. 
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1 .2  Summary  of  Full-Scale  Barge  Impact 
Experiments 

In  December  1998,  full-scale,  low-velocity,  controlled  impact  barge 
experiments  were  conducted  at  the  decommissioned  Gallipolis  Lock  at  Robert  C, 
Byrd  Lock  and  Dam,  Gallipolis  Ferry,  WV  (Figure  1.1).  The  primary  goal  of 
these  experiments  was  to  measure  the  actual  impact  forces  normal  to  the  wall 
using  a  load-measuring  device.  The  focus  of  the  experiments  was  to  obtain  and 
measure  the  baseline  response  of  an  inland  waterway  barge,  quantify  a  multiple- 
degree-of-freedom  system  during  the  impact,  and  investigate  the  use  of  energy¬ 
absorbing  fenders. 


Figure  1 .1 .  Robert  C.  Byrd  Lock  and  Dam,  Gallipolis  Ferry,  WV 


The  full-scale  experiment  used  a  15-barge  commercial  flotilla  (Figure  1 .2).  Each 
barge  was  a  jumbo  open-hopper  design  (35  by  195  ft)  with  rake  barges  at  the 
front  of  the  tow.  The  barges  were  ballasted  with  anthracite  coal  to  a  draft  of  9  ft. 
The  total  weight  of  the  flotilla  was  30,012  short  tons.  The  total  mass  is  1,865.59 
k-sec  /ft,  equal  to  the  total  weight  divided  by  the  gravitational  constant,  g. 

The  target  area  was  the  rigid  concrete  upper  guide  wall  (see  Figure  1.1), 
lacking  the  friction-reducing  steel  armor  found  on  modem  lock  walls.  A  total  of 
44  impact  experiments  were  successfully  conducted  against  the  unaltered  guide 
wall  and  a  prototype  fendering  system.  The  angle  of  impacts  ranged  from  5  to 
25  deg,  with  velocities  of  0.5  to  4  fps.  Of  these,  21  were  baseline  impact  tests, 
with  no  alterations  to  the  comer  barge.  For  the  remaining  23  impact  tests,  the 
comer  barge  was  replaced  with  an  identical  barge  that  had  been  fitted  with  an 
in-house-designed  load-measuring  beam.  This  device  will  be  discussed 
thoroughly  later  in  this  report. 
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Fourteen  of  the  impact  experiments  were  made  against  the  fendering  system. 
These  tests,  although  interesting,  are  not  the  focus  of  this  report  and  will  not  be 
discussed  further. 


Figure  1 .2.  Fifteen-barge  tow  configured  for  impact  experiments 


Instrumentation  mounted  on  the  barge  tow  consisted  of  55  measurements 
(data  channels)  for  the  baseline  experiments  and  44  measurements  for  tests  with 
the  load  bumper-equipped  barge.  The  instruments  included  15  axial  capacitive- 
accelerometers  on  the  impact  comer;  10  servo-accelerometers  to  monitor  motion 
at  key  locations  on  the  tow;  strain  gages  installed  on  principal  steel  in  the  impact 
area  (17  on  the  baseline  barge,  six  on  the  bumper-equipped  barge);  two  pressure 
cells  to  monitor  any  hydrodynamic  loading  as  the  tow  approached  the  wall;  nine 
clevis-pin  load  cells  in  the  lashings  that  secure  the  tow;  and  two  additional  high- 
range,  clevis-type  load  cells,  integral  to  the  load  bumper. 

The  tow-mounted  instrumentation  was  cabled  to  an  instrument  enclosure 
located  on  the  center-front  barge.  This  structure  housed  the  signal-conditioning 
electronics,  data  acquisition  computer,  and  system  operators  throughout  the 
project.  In  addition  to  the  tow-mounted  instrumentation,  servo-accelerometers 
were  located  on  either  side  of  the  joints  of  the  guide  wall  impact  monolith.  Data 
from  the  shore-mounted  instrumentation  will  not  be  discussed  in  this  report. 

The  precise  tow  position  during  the  44  impact  experiments  was  determined 
using  kinematic  differential  global  positioning  system  (DGPS)  instruments.  Two 
DGPS  units  were  deployed  along  the  port  side  of  the  comer  barge,  with  a  third 
unit  on  the  stem  area  of  the  center  barge.  These  units  recorded  tow  position  each 
second  to  an  accuracy  of  2  cm.  A  fourth  DGPS  unit  was  used  in  the  pilothouse  of 
the  towboat  to  monitor  and  control  the  tow’s  approach  speed.  The  accuracy  of 
this  unit,  operated  in  real-time  mode,  was  limited  to  3  m. 

These  experiments  were  conducted  over  a  period  of  3  days  in  December 
1998  on  a  highly  compressed  schedule.  The  first  day  consisted  of  installing  cable 
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and  instruments,  calibrating  systems,  and  surveying  instrument  locations  and 
other  key  points  on  the  tow.  The  baseline  experiments  were  begun  early  on  the 
second  day  and  were  completed  by  late  evening.  These  consisted  of  12  well- 
varied  impacts  (Experiments  1-12)  conducted  on  the  unurmored  lock  wall  (see 
Figure  1.3),  and  nine  impacts  (Experiments  13-21)  on  the  fender-equipped 
section.  Instruments  were  removed  from  the  comer  barge  at  the  end  of  the 
workday.  During  the  night,  the  tow  crew  replaced  the  standard  comer  barge  with 
a  second  barge  that  had  been  equipped  the  previous  week  with  the  comer  load 
beam. 


On  the  morning  of  the  third  day, 
instru mentation  was  remounted  to  the  new 
comer  barge,  with  testing  commencing 
shortly  before  noon.  Ten  more  impacts 
(Experiments  22-31)  were  made  against  the 
same  unarmored  section  of  wall  that  was  the 
focus  of  the  tests  conducted  the  previous  day. 

The  tow  was  moved  to  another  location  along 
the  wall  where  the  five  fender  impact  tests 
(Experiments  32-36)  were  conducted.  Then, 
the  tow  was  repositioned  to  complete  the  test 
matrix  at  the  primary  target  area  (Experiments 
37-44).  Testing  was  completed  in  the  early 
evening  of  the  third  day  of  testing. 

The  operational  procedures  followed  by 
the  tow  and  helper  boat  are  crucial  to 
understanding  the  force  system  applied  in 
these  experiments.  For  each  experiment  the 
tow  pilot,  assisted  by  the  helper  boat  positioned  at  the  starboard  bow,  would  align 
the  tow  with  the  target  at  the  desired  angle  of  attack.  This  angle  was  verified  with 
a  wall-based  survey  transit.  Once  positioned,  the  towboat’s  engines  were 
powered  up  to  bring  the  flotilla  to  the  required  speed,  which  was  verified  with  the 
GPS  unit  in  the  pilothouse  of  the  towboat.  Just  prior  to  impact,  the  tow  pilot 
would  cut  power  to  idle  and  drift  the  flotilla  into  the  wall.  Several  seconds  after 
impact,  the  towboat’s  engines  were  reversed,  while  the  helper  boat  would  throttle 
up  and  push  the  comer  against  the  wall  to  further  reduce  the  forward  momentum. 
The  tow  was  then  repositioned  for  the  next  experiment. 

As  stated  previously,  the  primary  goal  of  this  project  was  to  quantify  the 
normal  force  the  barge  tow  imparts  to  the  lock  wall.  The  design  load  for  this 
system  was  estimated  (prior  to  experiments)  to  be  less  than  1,000  kips.  Two 
approaches  to  making  this  measurement  were  considered.  The  first  option 
considered  was  to  build  a  load-measuring  platform  into  the  lock  wall.  This  option 
would  optimize  testing  efficiency,  as  no  swapping  of  the  comer  barge  was 
needed;  however,  it  was  unclear  how  reliably  the  1,100-ft  tow  could  be  steered 
onto  a  small  target.  A  further  concern  was  that  a  second  system  would  have  to  be 
built  into  the  fender  section. 


Figure  1 .3.  Barge  impacting  on  lock  wall 


without  armor  (note  concrete  dust) 
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The  second  option  was  to  fit  the  load-measuring  system  to  the  comer  barge 
itself,  using  a  heavy  curved  beam  as  a  load  bumper,  and  two  clevis-type  600-kip 
load  cells  configured  to  measure  the  normal  load  transferred  to  the  barge.  This 
second  option  was  chosen  because,  not  only  did  it  eliminate  the  concern  about 
hitting  a  target  zone,  it  was  also  much  less  complex  and  less  expensive  to  install. 
The  initial  load  bumper  design  called  for  a  pin  connection  at  the  front  of  the 
beam  with  a  roller  connection  at  the  rear  attachment.  With  this  design,  no  shear 
loads  would  be  transferred  to  the  barge  at  the  roller  support.  As  construction 
began  on  this  design,  it  became  clear  that  the  load  beam  would  be  poorly 
supported  at  the  rear,  and  that  this  weakness  would  most  likely  damage  the  front- 
mounted  load  cell.  The  solution  (since  the  investigators  were  not  overly 
interested  in  the  shear  component)  was  to  weld  both  the  front  and  rear  pin 
supports  to  the  barge.  This  created  an  indeterminate  pin-pin  connection  with  two 
reactions  at  each  pin. 


Figure  1 .4.  Load  bumper  mounted  to  the 
front  barge  corner 


Figure  1 .5.  Load  bumper  as  configured  for 


Experiments  22-27 


The  load  bumper,  installed  on  the 
barge  comer,  is  shown  in  Figures  1 .4  and 
1.5.  The  load  beam  was  constructed  of 
mild-steel  with  an  outer  radius  of  72.6  in., 
outer  arc-length  of  43.6  in.,  a  cross  section 
measuring  9  in.  in  width  by  5  in. 
in  height,  with  a  separation  between  the 
6-in.-diam  load  pins  of  35.5  in. 

It  should  be  noted  here  that  a  second 
load-sensing  device,  intended  as  a  backup, 
was  added  to  the  outer  surface  of  the  load 
beam.  This  device  used  three 
polyvinylidene  fluoride  (PVDF)  thin-film 
pressure  sensors  configured  in  a  1/8-in.- 
thick  flat-pack,  sandwiched  between  two 
steel  plates.  This  system  proved 
ineffective,  as  it  caused  excessive  shear 
drag  at  the  rear  brackets  that  secured  the 
PVDF  material  to  the  bumper.  For  this 
reason,  it  was  removed  after  Experi¬ 
ment  27,  creating  a  smooth,  continuous 
impact  surface.  Load  data  from 
Experiments  22-27,  especially  for  the 
lower  angles  of  impact,  was  adversely 
affected  by  this  drag,  and  should  be  used 
with  caution.  These  data  were  therefore 
discarded  for  the  subsequent  discussion  of 
the  forces  derived  from  the  F10  and  Fn 
clevis-pin  load  cells.  The  PVDF  sensors 
will  not  be  discussed  further  in  this  report. 
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The  approach  angle  and  velocity  for  the  12  most  credible  bumper 
experiments  (Experiments  28-31  and  Experiments  37-44)  are  summarized  in 
Table  1.1.  Impact  velocity  for  these  experiments  ranges  from  0.88  to  2.87  fps, 
with  approach  angles  ranging  from  8.8  to  21.1  deg.  This  table  was  constructed 
using  the  DGPS  data  from  the  two  instmments  located  on  the  comer  barge.  The 
time  and  point  of  impact  were  obtained  through  the  positioning  plots  of  the 
forwardmost  DGPS  station  (see  Appendix  A,  Figures  A.20-A.30),  This  station 
followed  the  general  course  of  the  tow  up  to  the  time  of  impact.  After  impact,  the 
front  of  the  barge  tracked  parallel  to  the  lock  wall. 


Table  1.1 

Impact  Velocity/ Angle  Data  for  Bumper  Experiment 

Experiment 

Number 

Impact 

Angle,  deg 

Velocity 

Velocity  Normal  to  the  Wail 

fps 

fps 

mph 

28 

9.7 

2.41 

1.64 

0.41 

0,28 

29 

12.7 

2.21 

1.50 

0.48 

0,33 

30 

12.2 

2.35 

1.60 

0.50 

0.34 

31 

10.6 

1.62 

1.10 

0.30 

0.20 

37 

10.3 

1.96 

1.33 

0.35 

0.24 

38 

11.9 

1.84 

1.25 

0.38 

0.26 

39 

14.1 

1.62 

1.10 

0.39 

0.27 

40 

17.5 

1.91 

1.30 

0.57 

0.39 

41 

8.8 

2.87 

1.95 

0.44 

0.30 

42 

17.5 

1.84 

1.25 

0.55 

0.38 

43 

21.1 

0.88 

0.60 

0.32 

0.22 

44 

20.9 

1.22 

0.83 

0.44 

0.30 

Once  the  time  of  impact  was  identified,  the  impact  angle  (the  angle  formed 
by  the  port  side  of  the  comer  barge  with  the  lock  wall)  was  determined  from  the 
DGPS  data  corrected  for  the  relative  positions  of  DGPS  stations  1  and  2  on  the 
barge  (see  Appendix  A,  Table  A.3).  This  angle  is  critical  to  the  bumper  geometry 
and  resulting  force  system.  Velocity  (actually  speed)  is  simply  calculated  from 
the  displacement  of  the  front  comer  GPS  unit  per  unit  time  (1  sec).  Data  typical 
of  the  Fj0  and  Fn  load  cells  are  presented  in  Figures  1.6  and  1.7.  During 
Experiment  41  (Figure  1.6),  the  angle  of  attack  is  very  shallow,  8.8  deg,  with  a 
substantial  velocity  of  2.87  fps.  The  point  of  impact  on  the  load  bumper  is 
slightly  in  front  of  the  F n  (rear)  load  cell.  As  expected,  the  loads  are  distributed 
with  the  major  portion  of  the  force  on  the  rear  cell.  Peak  loads  for  the  front  and 
rear  cells  were  measured  at  29.5  and  315.5  kips,  respectively  (see  Appendix  A, 
Table  A.1).  The  rise  time  for  this  measurement,  defined  here  as  Att  between 
contact  and  peak  load,  was  measured  at  164.8  msec  (see  Appendix  A,  Table  A.2). 
Total  contact  time  with  the  lock  wall  for  this  experiment  and  all  others  was 
approximately  9  sec. 

For  Experiment  43  (Figure  1.7),  the  impact  angle  is  21.1  deg,  but  with  a 
reduced  speed  of  only  0.88  fps.  For  this  orientation,  the  impact  point  on  the 
bumper  is  several  inches  forward  of  the  midpoint  between  the  load  pins.  In  this 
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significant  portion  of  this  report  will  deal  with  the  resolving  of  the  reactions  in 
the  bumper,  and  transferring  these  reactions  to  the  lock  wall. 

The  initial  orientation  of  the  bumper  relative  to  the  longitudinal  axis  of  the 
barges  was  adopted  to  be  54  deg  from  the  longitudinal  axis  (local  axis  of  the 
model)  of  the  barges.  Figure  1.8  shows  a  near-vertical  view  of  the  impact  site 
captured  from  video,  and  the  initial  orientation  of  the  bumper  and  the  recorded 
forces  at  the  pins,  which  were  assumed  to  be  in  the  radial  direction  of  the 
bumper.  The  force  notation  of  Spa,  Spb  used  in  the  bumper  analysis  is  F n  and  Fj0, 
respectively.  The  designations  “10”  for  forces  at  the  front  pin  and  “11”  for  the 
rear  pin  have  been  carried  through  from  the  measurement  numbers  FI0  and  Fn 
for  these  load-measuring  clevis-pins.  The  precise  orientation  of  the  bumper  on 
the  barge  is  critical  to  this  effort.  The  as-built  orientation  of  the  bumper  was  then 
developed  from  a  combination  of  design  drawings  and  documentary  photos.  The 
survey  data  were  intended  for  this  purpose;  however,  the  uncertainty  caused  by 
the  barges  shifting  and  the  tow  drifting  against  its  moorings  between  sightings 
compromised  the  accuracy  of  these  measurements  sufficiently  to  make  them 
unusable  for  this  purpose. 

Subsequently,  it  was  established  from  the  design  drawings  and  documentary 
photos  that  the  recorded  forces’  orientation  was  not  aligned  in  the  radial 
direction.  Taking  into  account  this  observed  discrepancy,  a  new  recorded  forces 
(Fu  and  FI0)  orientation  was  established.  Figure  1 .9  shows  the  bumper-arc 
geometry,  with  the  angle  and  separation  of  the  load  pins  indicated.  The  system  of 
forces  based  on  this  arc  geometry  is  also  shown  in  the  figure.  This  second 
configuration  was  analyzed,  and  considering  the  magnitude  of  the  angles 
indicating  the  support  reactions  orientation,  an  impossible  geometrical 
arrangement  was  produced. 

A  final  configuration  was  established  based  on  the  range  of  probable  angles 
for  the  forces’  orientations  relative  to  the  radial  direction,  the  location  of  the 
bumper  related  to  the  longitudinal  axis  of  the  barges,  and  the  appropriate 
coefficient  of  friction  between  concrete  and  steel.  This  configuration  is  shown  in 
Figure  1.10.  The  forces  F n  and  Fw  are  out  the  radial  direction  5.5  and  1.5  deg, 
respectively.  The  location  of  the  bumper  with  respect  to  the  longitudinal  axis  of 
the  barge  also  changes,  from  an  initial  approximation  of  54  to  57.5  deg.  It  will  be 
demonstrated  that  this  final  configuration  produces  reasonable  results  based  on 
the  values  of  the  coefficient  of  friction  between  the  wall  and  the  steel  bumper 
found  in  technical  literature,  and  using  the  fact  that  the  bumper  must  be  in 
compression  during  the  impact  process.  (Note  that  it  is  impossible  for  the  bumper 
to  be  in  tension  during  the  impact  process.) 

The  systems  of  forces  for  each  of  the  three  configurations  for  the 
indeterminate  load  bumper  are  shown  in  Figures  1.8-1.10.  A  total  of  six  forces 
are  present: 

Two  known  forces:  Fw  and  Fn 

(labeled  in  the  figures  as  SPB  and  SPA,  respectively) 

Four  unknown  forces:  Fw,  Sw  (friction  force),  Fm,  and  FNB 
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a.  Overhead  photo  of  barge  with  bumper  about  to  impact 
the  lock  wall 


b.  Arc  geometry  used  for  force  analysis 


c.  Forces  acting  on  load  bumper — initial  configuration 
Figure  1.8.  Initial  configuration  and  recorded  load 
direction 
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b.  Arc  geometry  used  for  force  analysis 


90-0 


c.  Forces  acting  on  load  bumper 
Figure  1 ,9.  Second  configuration  and  recorded  load 
direction 


a. 


Overhead  photo  of  barge  with  bumper  about  to  impact 
the  lock  wall 


a.  Overhead  photo  of  barge  with  bumper  about  to 
impact  the  lock  wall 


c.  Forces  acting  on  load  bumper 
Figure  1 .10.  Final  configuration  and  recorded  load 
direction 
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The  pair  of  unknown  forces  acting  at  the  point  of  contact  (during  impact) 
between  the  bumper  and  the  wall  is  designated  as  Fw,  Sw.  The  unknown  forces 
acting  in  the  direction  normal  to  the  measured  pin  reaction  forces  are  Fm  and 
Fm  at  each  of  the  pins.  Four  equations  will  be  required  to  solve  for  these  four 
unknown  forces.  The  first  three  governing  equations  are  the  equilibrium 
equations: 


T,Fx=o 

(1.1) 

M 

*”ri 

V; 

II 

o 

(1.2) 

2>„=o 

(1.3) 

Two  methods  will  be  used  to  formulate  the  fourth  equation.  One  of  the 
methods  investigated  assumes  a  value  for  the  coefficient  of  friction,  ju,  between 
concrete  and  steel.  This  equation  takes  the  form 

SW=V*FW  (1.4) 

This  formulation  will  be  referred  to  in  the  report  as  the  equilibrium 
formulation  with  a  fixed  assumed  value  for  the  coefficient  of  friction.  Note  that 
this  method  fixes  the  value  for  Sw  to  Fw  by  this  equation  using  the  value  assumed 
for  fi.  This  relationship,  in  conjunction  with  the  three  equations  of  equilibrium,  is 
developed  in  Chapter  3  to  solve  for  the  four  unknowns. 

An  alternative  formulation  investigated  during  the  course  of  this  research  for 
the  required  fourth  equation  is  referred  to  in  this  report  as  the  energy  method. 

The  energy  method  uses  Castigiliano’s  theorem,  which  states  that  the  partial 
derivative  of  the  total  strain  energy  of  the  arc  of  the  beam  with  respect  to  one 
external  load  is  equal  to  the  displacement  in  the  direction  of  the  external  applied 
load.  In  this  case  (a  first-degree  indeterminate  system),  one  of  the  support 
reactions  was  selected  as  the  redundant.  All  the  internal  forces  (axial,  shear,  and 
moment)  were  expressed  in  terms  of  the  redundant  force  located  in  cell  10  (the 
front  load  cell).  The  additional  equation  is  the  displacement  expression  obtained 
from  Castigiliano’s  theorem,  which  is  equal  to  zero  because  the  redundant  was 
taken  in  the  support.  The  relationship  (developed  in  Chapter  4)  results  in  the 
fourth  equation,  which  is  used  in  conjunction  with  the  three  equilibrium 
equations. 


1.3  Impulse  and  Linear  Momentum  Principle 

This  report  presents  the  analysis  of  the  experimental  results  to  obtain  the 
applied  forces  to  the  lock  wall  due  to  the  impact  of  a  barge  flotilla.  The  study  of 
the  motion  of  particles  is  based  upon  axioms  and  laws  of  nature  gained  from 
experience,  first  formulated  by  Isaac  Newton.  Newton’s  laws  assume  that  one 
can  define  a  frame  of  reference,  which  is  fixed  in  space,  called  an  inertial  frame 
of  reference  or  a  Newtonian  frame  of  reference.  Only  in  this  frame  of  reference 
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will  Newton’s  second  law  be  valid.  For  most  purposes,  a  frame  of  reference 
attached  to  a  star  will  serve  as  an  inertial  reference  frame.  However,  any 
nonrotating  reference  frame  moving  at  a  constant  velocity  relative  to  an  inertial 
frame  can  also  be  considered  an  inertial  frame. 

The  concept  of  mass  arises  in  two  of  Newton’s  laws.  In  the  second  law, 
inertial  mass  is  considered  to  be  a  measure  of  a  particle’s  resistance  to 
acceleration.  In  Newton’s  fourth  law,  gravitational  mass  is  defined  as  the 
property  of  the  particle  that  influences  its  gravitational  attraction.  Newton  further 
assumed  that  these  two  concepts  of  mass  were  equivalent.  The  mathematical 
form  of  Newton’s  second  law  states  that  a  resultant  external  force  applied  to  a 
body  is  equal  to  the  mass  of  the  body  multiplied  by  the  absolute  acceleration  the 
body  experiences.  Also,  it  can  be  expressed  in  terms  of  the  absolute  velocity  of 
the  body  by  introducing  the  first  derivative  with  respect  to  time  of  the  velocity, 
which  is  the  acceleration. 

One  useful  tool  that  can  be  derived  from  Newton’s  second  law,  F  =  ma,  is 
obtained  by  integrating  both  sides  of  the  equation  with  respect  to  time.  This 
integration  can  be  done  only  if  the  forces  acting  on  the  particle  are  known 
functions  of  time.  The  external  forces  acting  on  the  particle  change  the  linear 
momentum.  The  mathematical  form  of  the  resulting  expression  after  the  process 
of  integration  states  that  the  impulse  during  a  period  of  time  due  to  the  applied 
impulsive  force  is  equal  to  the  difference  in  linear  momentum  during  the  same 
interval  of  time. 

This  relationship  establishes  the  impulse  and  linear  momentum  principle. 

The  units  of  both,  impulse  and  momentum,  are  force  and  time,  and  therefore, 
impulse  and  momentum  are  expressed  in  N»s  or  kips#s.  The  impulsive  force  is  a 
function  of  time  and,  in  general,  varies  during  its  period  of  application.  A  large 
force  that  acts  over  a  short  period  of  time  is  called  an  impulsive  force  and  occurs 
during  phenomena  such  as  the  impact  of  a  bat  with  a  ball,  collisions  of  cars,  or  a 
barge  impacting  a  lock  wall.  If  the  average  impulse  force  is  zero,  the  linear 
momentum  does  not  change  during  that  interval  of  time. 

The  impulse  and  linear  momentum  principle  will  be  used  in  the  next  chapters 
to  develop  a  tool  to  determine  the  normal  force  in  a  lock  wall  due  to  the  impact 
of  a  barge  flotilla.  It  will  be  done  by  using  the  results  from  the  full-scale,  low- 
velocity,  controlled  barge  impact  experiments  at  Robert  C.  Byrd  Lock  in 
December  1998.  Several  models  were  generated  to  determine  the  load  that  a 
barge  applied  to  a  lock  wall,  based  on  these  experiments.  The  load  between  the 
barge  and  the  wall  was  transferred  by  a  curved  beam  mounted  on  the  comer  of 
the  barge,  as  described  in  the  preceding  section.  This  indeterminate  beam  was 
analyzed  using  two  approaches. 

The  first  model  that  was  developed  considers  the  three  global  equilibrium 
equations — the  equation  that  relates  the  shear  force  in  the  wall  and  the  normal 
force  in  the  wall  due  to  the  contact,  and  the  axis  transformation  at  both  supports. 
Both  terms,  the  normal  and  shear  force  in  the  wall,  were  related  by  the  use  of  the 
dynamic  coefficient  of  friction.  A  review  of  the  technical  literature  indicates  that 
the  coefficient  for  steel-concrete  interface  typically  ranges  in  value  from  0.3  to 
0.7.  That  is,  the  shear  force  is  0.3  to  0.7  the  magnitude  of  the  normal  force. 
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Values  for  the  unknown  forces  were  obtained  by  solving  Equations  1.1  through 
1.4.  As  shown  in  Figure  1.1 1,  the  known  forces  are  the  FI0  and  Fu  ( SPB  and  SPA, 
respectively)  measured  during  the  experiments,  and  the  unknowns  are  the  shear 
force  in  the  wall,  the  normal  force  in  the  wall,  and  the  other  two  reactions  at  the 
pin  supports.  This  model  will  be  developed,  explained,  and  validated  in 
Chapter  3. 

The  second  model  consists  of  using  the  three  global  equilibrium  equations 
(Eqs.  1.1-1 .3),  the  equations  obtained  from  an  energy  method,  and  the  axis 
transformation  at  both  supports.  This  means  that  the  beam  was  considered  as  a 
deformable  body  and  not  as  a  rigid  body  (as  is  the  case  for  the  first  model).  This 
additional  (fourth)  equation  was  obtained  from  the  global  strain  energy  of  the 
beam.  The  methodology  used  is  also  known  in  the  mechanics  of  material  field  as 
the  Castigliano’s  theorem,  which  states  that  the  displacement  in  the  direction  of 
an  external  applied  load  can  be  calculated  by  differentiating  the  internal  force 
equations  with  respect  to  the  external  load  in  the  direction  of  the  required 
displacement.  In  this  problem,  the  selected  force  was  one  of  the  reactions 
providing  that  the  displacement  should  be  equal  to  zero.  In  Chapter  4,  this  model 
is  presented  and  validated  against  results  from  finite  element  models  of  the 
bumper. 

During  the  course  of  this  research,  three  possible  reaction  force 
configurations  acting  on  the  bumper  were  studied,  as  was  mentioned  earlier. 
During  the  course  of  the  study,  these  conclusions  were  reached: 

a.  The  forces  were  not  measured  in  the  radial  direction. 

b.  A  possible  configuration  was  estimated  using  design  drawings  and 
documentary  photos. 

c.  A  third  configuration  was  adopted  based  on  the  reasonable  values  of 
coefficient  of  friction  between  concrete  and  steel  and  with  the  fact  that 
the  bumper  must  be  in  compression  during  impact. 


Chapter  5  presents  the  comparisons  of  these  three  models.  Each  of  the 
configurations  is  distinguished  by 

a.  The  fixed  location  where  the  bumper  was  attached  relative  to  the 
longitudinal  axis  of  the  barges. 

b.  The  orientations  of  the  measured  forces  F!0  and  FIt. 


These  three  configurations  are  shown  in  Figures  1.12-1.14,  respectively. 
Chapter  6  presents  the  use  of  the  Impulse  and  Linear  Momentum  Principle  as 
applied  to  our  case,  and  the  third  configuration  studied,  shown  in  Figure  1.14.  In 
this  case,  structural  indeterminacy  is  related  to  the  three  external  forces  that  were 
applied  during  the  impact  and  were  not  measured.  These  forces  are  the  drag 
force,  the  towboat  force,  and  the  helper  boat  force.  With  this  principle,  we  have 
fewer  equations  than  unknowns — thus,  an  indeterminate  system.  A  simple 
approach  was  developed  to  calculate  the  forces  applied  to  the  wall  based  on  the 
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Structure 


Idealization 
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Figure  1.11.  Idealization  and  free-body  diagram  of  the  bumper 
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Global  Axis 


a.  Load  cell  layout 


b.  Configuration  and  forces  acting  on  load  bumper 

Figure  1.12.  Initial  configuration  and  forces  acting  on  load  bumper 
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b.  Configuration  and  forces  acting  on  load  bumper 
Figure  1 .13.  Second  configuration  and  forces  acting  on  load 
bumper 
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Global  Axis 


a.  Load  cell  layout 


b.  Configuration  and  forces  acting  on  load  bumper 
Figure  1.14.  Third  and  final  configuration  and  forces  acting  on  load 
bumper 
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calculated  force  time-histories  normal  to  the  wall  for  eight  of  the  12  instrumented 
bumper  experiments  presented  in  Chapter  5.  Chapter  7  suggests  the  layout  for 
future  experiments  that  would  avoid  the  highest  degree  of  indeterminacy  of  the 
system.  Conclusions  and  recommendations  for  the  numerical  model  developed  in 
this  study  are  presented  in  Chapter  8. 

Supplemental  information  on  this  research  is  presented  in  the  four  report 
appendixes:  A,  Summary  of  experimental  results;  B,  FORTRAN  source 
programs  and  Maple  worksheets;  C,  Numerical  integration;  and  D,  Statistical 
theory. 
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Impulse  and  Linear 
Momentum  Principle 
Applied  to  Barge  Flotilla 
Model 


The  principle  of  impulse  and  linear  momentum  for  a  system  of  particles 
moving  relative  to  an  inertial  reference  frame,  as  illustrated  in  Figure  2.1,  is 
obtained  from  the  equation  of  motion 

(2.1) 

which  can  be  expressed  as 


n  n 

15 -I* £ 

i=l  i= 1  dt 
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The  term  on  the  left  side  represents  only  the  sum  of  all  external  forces  acting 
on  the  system  of  particles.  The  internal  forces  between  the  particles  do  not  appear 
with  this  summation,  since  by  Newton’s  third  law  they  occur  in  equal  but 
opposite  collinear  pairs  and  therefore  cancel  out.  Multiplying  both  sides  by  dt, 
and  integrating  between  the  limits  time  t  =  t\,  velocity  v(-  =  (v,^  ,  and  time  t  =  1 2, 

velocity  v,  =  (v,)2,  yields 


X  mi  fa  \  +  X  f F‘dt  =  X  mi  fa  )l  (22) 

i=l  i=l  t[  i= 1 

Equation  2.2  states  that  the  initial  momentum  of  the  system  added  vectorially 
to  the  impulse  of  all  the  “external  forces”  acting  on  the  system  during  the  time 
period  t\  to  h  is  equal  to  the  system  final  linear  momentum,  where  t\  is  the  time 
when  impact  of  the  flotilla  with  the  wall  begins  and  t2  is  the  time  when  the  impact 
terminates  (approximately  9  sec  after  tj  for  the  1 1  most  credible  bumper 
experiments).  If  each  of  the  vectors  in  Equation  2.2  is  resolved  into  its  x,  y,  and  z 
components,  we  can  symbolically  write  the  following  three  scalar  equations: 


n  (  \  n  \  n  (  \ 
XmAv4+X  K^=XmAv2 

i=l  1=1  tx  t’=l 

(2.3) 

2>/(»v,  1 + Z  V/' = 2>/(vj,  \ 

(=1  i= 1 i=l 

(2.4) 

Xm*k)i+X  )Fzdt=lLmkz\ 

i- 1  /=1 1  i=l 

(2.5) 

These  equations  represent  the  impulse  and  linear  momentum  principle  for  the 
system  of  particles  in  the  x,  y,  and  z  directions,  respectively. 

The  flotilla-wall  system  consists  of  a  rigid  concrete  wall,  and  15  barges 
joined  together  by  lashings  form  the  assumed  rigid  body  flotilla.  The  dynamic 
event  occurs  during  the  impact  of  the  flotilla  with  the  wall.  Figure  2.2  shows  the 
system  to  be  analyzed  in  the  global  coordinates  system. 

The  free-body  diagram  of  the  system  is  shown  as  Figure  2.3.  Applying  the 
Impulse  and  Linear  Momentum  Principle  in  the  global  coordinates,  parallel  and 
perpendicular  to  the  wall,  respectively,  produces  the  following  expressions: 


m  flotilla  v\x 


2  r  r  r 

j '  F tow  —  J  Swdt  —  j  F  helper  t  dt  ~  J  F drag  ^dt 

^1  h  h  h 


=  m  flotilla  v2x  ® 


(2.6) 
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Figure  2.2.  Flotilla-wall  system 


Figure  2.3.  Free-body  diagram  of  flotilla-wall  system 

h  h  U  1 2 

m flotilla  vly  +  \Ftowydt  ~  \fwdt  +  J F helper ydt  ~  J F drag  dt 

h  h  h  ”  r,  ’  (2.7) 

=  m  flotilla  v2  y  ^  ® 

where 

m flotuia  =  mass  of  the  flotilla  (including  the  tow  boat  and  the  helper 
boat) 

Ftow  =  force  the  tow  boat  applies  to  the  flotilla  during  impact 
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Sw  =  shear  force  parallel  to  the  wall  during  impact  with  the  front 
comer  of  the  barge  as  the  flotilla  slides  along  the  wall 

F helper  =  force  the  helper  boat  applies  to  the  flotilla  during  impact 

Fdmg  =  total  drag  force  on  the  flotilla  as  it  moves  through  the  water 
during  impact 

Fw  =  force  component  normal  to  the  wall  during  impact  with  the 
front  comer  of  the  barge  as  the  flotilla  slides  along  the  wall 

The  force  provided  by  the  towboat  can  be  estimated  using  the  power 
definition.  Power  is  defined  as  the  amount  of  work  performed  per  unit  of  time, 
and  can  be  expressed  as 

P  =  =  F  •  v  =  Fvcos^  (2.8) 

dt 

where  <f>  is  the  angle  between  the  force  and  velocity  vector. 

If  the  power  is  known  during  the  time  interval  from  t,  and  t2 ,  then,  F,ow  can 
be  calculated  as 


Ftow 


—  </>  =  0° 
vi 


(2.9) 


Then,  the  components  of  this  force  in  the  global  coordinate  system  are 

Ftow  =^-cos0  (2.10) 

1  Vi 

Ftow  =  ^sin0  (2.11) 

y  n 

An  additional  modification  to  these  expressions  has  to  be  made  to  take  into 
account  the  fact  that  the  power  that  affects  the  system  being  considered  is  the 
power  provided  by  the  propeller  to  the  water,  and  not  the  power  provided  by  the 
engine  to  the  propeller.  However,  similar  expressions  can  be  used  to  define  the 
force  applied  to  the  water  by  the  propeller.  In  the  same  way,  the  force  provided 
by  the  helper  boat  can  be  estimated  using  this  approach.  Then,  the  force  provided 
by  the  helper  boat  to  the  system  can  be  calculated  by  the  following  expressions: 


F helper  .  Q 

F helper, ;  =  sin^ 

(2.12) 

Fhelper  Q 

F helper  =  COS0 

(2.13) 
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Note  that  to  use  these  expressions  for  the  forces  Ftow  and  Fhdper,  the  power  terms 
P tow  and  P helper  must  reflect  the  power  that  the  propellers  of  the  tow  and  helper 
boats  provide  to  the  water  (and  not  the  power  the  engines  provide  to  the  propeller 
shafts  of  the  tow  and  helper  boats)  during  the  approximately  9  sec  of  impact  with 
the  lock  wall.  Measurements  leading  to  the  definition  of  these  power  terms  were 
not  made  during  the  full-scale,  low-velocity,  controlled  barge  impact  experiments 
at  Robert  C.  Byrd  Lock  in  December  1998. 

The  drag  force  is  a  resistance  force  that  the  water  presents  opposite  to  the 
motion  of  the  flotilla.  To  calculate  this  force,  it  is  necessary  to  introduce 
hydrodynamics  terminology  and  technology  (e.g.,  see  Martin  1989,  p  33). 

All  these  forces  (Flom  Sw,  Fhetper,  Fdrag,  and  Fw)  produce  impulse  to  the  system. 
In  addition,  the  added  mass  of  water  due  to  the  approach  of  the  barge  to  the  wall 
must  be  included.  The  additional  mass  added  to  the  system  due  to  this  effect 
compared  with  the  mass  of  the  flotilla  is  unknown  at  this  time.  This  effect  is 
presented  in  Figure  2.3  as  a  concentrated  mass  in  the  point  of  impact  between  the 
bodies.  This  water  mass  has  the  effect  of  cushioning  the  impact. 

The  hydrodynamic  effect  due  to  the  approach  of  the  flotilla  to  the  wall 
(motion  of  mass  of  water)  can  be  included  as  a  change  in  linear  momentum  of  the 
water  particle.  Now,  considering  that  the  flotilla  travels  almost  with  constant 
velocity  before  and  after  impact  in  these  experiments,  the  Impulse  and  Linear 
Momentum  Principle  for  this  system  of  particles  can  be  illustrated  as  shown  in 
Figure  2.4. 

Following  Figure  2.4,  Equations  2.6  and  2.7  can  be  written  as 

h  h  h 

mfioMiavi  cos0+  \FwwcosQdt-  \Swdt-  } Fhdper  sin Qdt 

'i  '/  t, 

h 

-  \  F drag  cos  Qdt  =  m  flotitta  v2  cos  e 


im flotilla  -  m water  V l  sin0+] Ftow  sin Odt  -  } Fwdt 
*t  h 

h  t2 

+  iF helper  cos  Odt  -  \Fdms  sinddt  =  mflotillav2  sin0 
*t  h 

where  the  effect  of  the  (hydrodynamic)  cushioning  “force”  was  introduced, 
reducing  the  initial  linear  momentum  by  subtracting  the  mass  of  water. 

Finally,  rearranging  terms,  we  obtain 


(2.14) 


(2.15) 


m 


flotilla^  i  cos0-mflotiUav2  cos0+  \Ftowco$9dt-  jF,w,  s  in  Qdt 


<2  h 

-  \F drag  cos  6dt=  jswdt 
• I  t , 


(2.16) 
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{m  flotilla  ~  m  water  K  sind  +  }Ftow  sin  ddt  +  }  Fhelper  cosddt 

\  \  (2-17) 

-  lFdmg  sin  ddt  =  }  F wdt 

*i  >i 

One  important  detail  regarding  Equation  2.17  is  that  the  velocity  normal  to 
the  wall  during  impact  is  equal  to  zero.  Thus,  the  linear  momentum  normal  to  the 
wall  is  also  zero  during  impact  (that  is,  from  f,  to  t2).  Then,  Equations  2.16  and 
2.17  provide  the  impulse  due  to  the  shear  and  normal  force  in  the  wall  during 
impact,  respectively. 

Measurements  leading  to  the  definition  of  the  forces  Ftow,  Fhelper,  and  Fdrag 
and  the  hydrodynamics  added-mass  term  mwater  were  not  made  during  the  full- 
scale,  low-velocity,  controlled  barge  impact  experiments  at  Robert  C.  Byrd  Lock 
in  December  1998.  However,  measurements  were  made  during  these  experiments 
leading  to  the  definition  of  forces  Sw,  and  Fw  and  will  be  discussed  in  subsequent 
chapters. 
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3  Equilibrium  Formulation 
with  Assumed  Coefficient 
of  Friction 


3.1  Background 

To  apply  the  Impulse  and  Linear  Momentum  Principle  as  presented  in 
Equations  2.16  and  2.17,  it  is  necessary  to  determine  the  normal  force  time- 
history  applied  to  the  wall  and  the  shear  force  acting  along  the  wall  during  the 
impact.  The  determination  of  these  time-histories  can  be  done  using  the  time- 
histories  recorded  during  the  experiments.  These  time-histories  were  obtained 
from  the  pin  load  cells  located  at  the  supports  of  the  curved  beam  bumper.  The 
recorded  force  time-histories,  designated  as  F,0  and  Fu,  are  presented  in 
Appendix  A.  The  model  presented  herein  considers  the  recorded  load  in  the 
radial  direction  of  the  circular  arc  beam.  It  was  the  initial  specification,  from 
those  in  charge  of  the  experiments,  of  the  first  assumed  geometrical  configuration 
to  be  analyzed  in  this  study.  Using  this  first  model,  a  procedure  to  determine  the 
unknown  external  loads  for  the  indeterminate  circular  arc  is  presented  in  this 
chapter.  Later,  in  Chapter  5,  the  analysis  using  this  procedure  with  two  other 
configurations  will  be  presented. 

The  structural  “bumper”  system  consists  of  a  circular  arc  beam  with  a 
constant  rectangular  cross  section  made  of  A-36  steel.  The  instrumented  bumper 
is  shown  in  Figure  3.1.  This  pin-pin  connection  produces  an  indeterminate 
structure  to  the  first  degree.  As  shown  in  Figure  3.2,  the  four  unknowns  are  the 
normal  force  in  the  wall,  the  shear  force  in  the  wall,  and  two  support  reactions 
that  were  not  recorded  during  the  experiments  ( FNI0  and  FNII).  The  two  known 
forces  are  the  FI0  and  Fn>  which  were  assumed  in  this  first  model  to  act  in  the 
radial  direction.  Using  Figure  3.3,  the  three  equilibrium  equations  in  global 
coordinates  are 


— >  © 

1^=0.-. 

Ax  +  Bx-  Sw  -  0 

(3.1) 

T® 

2X =<>••• 

Ay  +  By  ~  Fw  =  0 

(3.2) 
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Figure  3.1 .  Instrumented  circular  arc  beam  (bumper) 


=0.-.  FWR  cos(97  -  0)  -  SH,/?  [1  -  sin  (97  -  6)] 


-  ByR  [cos  (54  +  0)  +  cos  (97  -  0)] 


(3.3) 


-BXR [sin  (97  -  0)  -  sin  (54  +  0)]  =  0 

In  addition  to  the  three  equilibrium  equations,  it  is  possible  to  introduce  the 
axis  transformation  from  local  to  global  coordinates.  These  four  new  equations 
are 


4c  =  Fm  cos[9-l)-SPA  sin(6>-7)  ■  (3.4) 

=  Fna  sin(0-7 )  +  SPA  cos(0-7)  (3.5) 

Bx  =  -Fnb  cos(36-  9)  +  SPB  sin(36  -  9)  (3,6) 

By  =Fm  sin(36-^)  +  5PB  cos(36-0)  (3,7) 


Note  that  in  these  four  equations,  the  use  of  angles  7  and  36  deg  establishes 
that  the  orientation  of  the  bumper  is  54  deg  relative  to  the  longitudinal  axis  of  the 
barges  (local  axis)  and  also  defines  the  orientation  of  the  pairs  of  reaction  forces 
at  each  pin,  which  can  be  observed  in  Figure  3.3. 
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Figure  3.2.  Structure,  idealization,  and  tree-body  diagram  of  the  bumper 
for  the  assumed  coefficient  of  friction  method 
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Figure  3.3.  (Concluded) 


One  additional  equation  is  introduced  to  produce  a  determinate  system.  As 
was  found  in  the  review  of  the  technical  literature,  the  value  for  dynamic 
coefficient  of  friction  between  steel  and  concrete  is  usually  within  the  range  of 
0.3  to  0.7.  An  additional  equation  being  used  is 


Sw=MkFw  (3.8) 

Solving  these  eight  equations,  the  problem  is  transformed  to  a  determinate 
system.  Now  the  known  forces  are  SPA  (load  cell  11)  and  SPB  (load  cell  10),  which 
are  the  radial  forces  measured  during  the  experiments  {F,,  and  Fw,  respectively). 
The  unknowns  are  the  reactions  in  global  coordinates  (Ax,  Ay,  Bx,  and  By),  the 
forces  in  the  wall  (Fw  and  Sw),  and  the  two  reactions  normal  to  the  cross  section 
of  the  beam  in  local  coordinates  (FNA  and  Fa®). 

Three  computer  software  programs  were  developed  to  calculate  these  forces. 
First,  a  spreadsheet  using  Excel  was  built.  Second,  a  Maple  worksheet  was 
created  to  calculate  numerically  and  symbolically  the  expressions  for  Fw  and  Sw. 
Finally,  a  FORTRAN  computer  program  was  developed  to  compute  the  load 
time-histories  for  the  experiments. 

To  determine  if  the  results  provided  by  the  proposed  formulation  are 
accurate,  two  finite  element  models  of  the  bumper  were  studied  using  the  com¬ 
puter  programs  SAP2000  (Computer  &  Structures,  Inc.  2000)  and 
VisualAnalysis  (Integrated  Engineering  Software  1998).  The  models  were 
constructed  fixing  the  approach  angle  0to  18  and  21.5  deg.  Thus,  we  have  the 
same  structure  but  with  different  positions  in  the  global  coordinate  system. 
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Figure  3.4  shows  the  idealization  of  the  bumper  structure,  and  Figures  3.5  and 
3.6  show  the  arc  layout  for  the  Visual  Analysis  and  SAP2000  model  for  6 
=  18  deg,  respectively.  The  properties  of  the  beam  elements  used  in  these  models 
are  as  follows: 

Cross-sectional  area  =  45  in.2  Moment  of  inertia  =  303,75  in.4 

Modulus  of  elasticity  =  29,000  ksi  Cord  length  =  35.5  in. 


Figure  3.4.  Idealization  of  the  bumper  structure 
for  the  finite  element  analysis 


Figure  3.5.  VisualAnalysis  arc  bumper  model  for  0  =  1 8  deg 


Figure  3.6.  SAP2000  arc  bumper  model  for  0  =  1 8  deg 
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The  sectional  properties  were  calculated  using  the  dimensions  5  by  9  in.  for 
the  arc  bumper.  The  external  applied  loads  were  1 .0  kips  for  Fw  and  0.50  for  Sw. 
These  are  assumed  values,  which  correspond  to  a  coefficient  of  friction  of  0.50. 
The  support  reactions  were  obtained  from  the  finite  element  analysis.  Using  the 
support  reactions  obtained  from  the  VisualAnalysis  and  S  AP2000  analyses,  the 
proposed  formulation  was  used  to  calculate  the  external  applied  loads  ( Fw  and 
Sw ).  The  results  produced  by  each  of  the  three  models  are  presented  in  Tables  3.1 
and  3.2. 


The  procedure  used  in  these  analyses  is  as  follows: 

a.  An  assumed  Fw  and  Sw  value  was  adopted  based  on  a  coefficient  of 
friction  of  0.5.  These  are  the  shaded  cells  (input  data  for  the  model)  in 
the  last  two  rows  in  Tables  3.1  and  3.2. 

b.  The  Gi  (load  cell  1 1)  and  02  (load  cell  10)  values  are  the  resulting 
approach  angles  used  in  the  finite  element  analysis  at  each  support. 

These  are  different  from  the  exact  values  0  =  18  deg  and  0  =  21 .5  deg, 
for  Table  3.1  and  3.2,  respectively,  because  the  arc  is  modeled  using 
straight-line  elements.  One  important  observation  is  that  the  average  of 
(0i  +  02)  produces  the  exact  approach  angle  value  used  in  the  equilibrium 
and  fixed  coefficient  of  friction  method.  [For  example:  (21 .030  deg 

+  21.970  deg)/2  =  21.5  deg.] 

c.  Using  these  values  as  input  data,  the  supports’  reactions  were  calculated 
in  the  global  direction  using  the  SAP2000  and  VisualAnalysis  programs, 
as  shown  in  Figures  3.5  and  3.6. 

d.  The  forces  that  are  needed  in  the  Equilibrium  and  Fixed  Coefficient  of 
Friction  method  are  the  radial  direction  forces  (F10  and  Fu) — the  same 
as  the  recorded  forces  during  the  experimental  procedure  in  1998.  Using 
Equations  3.4-3.7,  the  values  of  SPA  (load  cell  1 1)  and  SPP  (load  cell  10) 
can  be  calculated. 

e.  Using  these  two  forces  and  an  approximate  approach  angle  as  input  data 
(shaded  cells  in  the  Equivalent  and  Fixed  Coefficient  of  Friction 
method),  the  Fw  and  Sw  forces  are  calculated. 

f  Finally,  the  resulting  coefficient  of  friction  obtained  using  the  Equivalent 
and  Fixed  Coefficient  of  Friction  method  is  the  same,  |ik  =  0.5,  as  was 
assumed  initially  in  the  finite  element  analyses. 


From  Tables  3.1  and  3.2,  we  can  observe  good  agreement  between  the  first 
developed  model  (equilibrium  formulation  with  a  fixed  value  for  the  coefficient 
of  friction)  and  the  finite  element  models.  One  important  point  is  that  the  arc  was 
modeled  in  the  finite  element  model  using  straight  lines.  Because  of  this,  the 
initial  and  final  angles  for  the  arc  at  the  support  reaction  differ  for  the  equilibrium 
and  fixed  coefficient  of  friction  model  by  a  half-degree.  For  the  model  developed 
in  this  research,  the  angle  ©i,  and  02  changed  by  small  amounts,  as  shown  in 
Tables  3.1  and  3.2.  These  small  changes  did  not  affect  the  computed  results.  This 
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analysis  provides  an  independent  check  of  the  accuracy  of  the  Equilibrium  and 
Fixed  Coefficient  of  Friction  method  because,  using  the  results  from  the  finite 
element  analyses  as  input  data  in  the  proposed  method,  the  same  coefficient  of 
friction  was  obtained  (juk  =  0.5). 

This  numerical  model  solves  for  the  unknowns,  by  using  the  SPA  (at  pin  A 
where  load  cell  1 1  is  located)  and  Spb  (at  pin  B  where  load  cell  10  is  located)  as 
input  data.  In  any  traditional  arc  beam  the  known  values  are  the  applied  forces, 
and  the  unknown  values  are  the  support  reactions.  This  means  that,  in  this 
problem,  we  are  going  in  an  inverse  direction  as  in  the  traditional  problems.  One 
special  case  was  found  in  the  solution  process  of  the  numerical  model  presented 
in  this  chapter,  and  must  be  mentioned.  The  special  case  occurs  when  the 
approach  angle  0  has  a  value  that  produces  impact  at  the  midpoint  between  the 
supports,  and  the  same  magnitude  for  the  input  data  (SPA  at  pin  A  where  load 
cell  1 1  is  located  and  SPP  at  pin  B  where  load  cell  10  is  located)  forces  is  used.  In 
the  case  presented  in  this  chapter,  the  value  of  0  that  produces  this  situation  is 
18.5  deg.  If  the  same  input  forces  are  used  for  an  approach  angle  of  1 8.5  deg,  the 
model  will  produce  unreasonable  results  for  the  unknowns.  This  results  because 
it  is  impossible  to  have  the  same  support  reactions  if  the  two  applied  loads  are 
placed  at  the  center  of  the  arc.  The  behavior  of  the  arc  when  the  horizontal  load 
Sw  is  applied  is  to  produce  a  different  distribution  of  support  reactions  as  if  only 
the  normal  force  Fw  is  applied. 


3.2  Baltay  and  Gjelsvik  (1990)  Concrete-to-Steel 
Coefficient  of  Friction  Results 

In  order  for  the  model  presented  in  this  chapter,  which  is  based  on  an 
assumed  coefficient  of  friction,  to  be  successful,  it  is  necessary  to  adopt  a 
realistic  value  for  the  coefficient  of  friction.  Based  on  the  literature  research  that 
was  conducted,  the  range  of  coefficient  that  could  be  acceptable  ranges  in  value 
from  0.13  to  0.70.  Table  3.3  presents  different  values  for  the  static  coefficient  of 
friction  and  for  the  kinetic  coefficient  of  friction,  found  or  adopted  in  other 
research.  In  the  current  study,  the  kinetic  coefficient  of  friction  was  in  this 
magnitude  for  two  reasons.  First,  the  low-velocity  problem  moves  this  coefficient 
to  approach  the  static  value.  Second,  the  contact  area  (5  by  7.25  in.  =  36.25  in.2) 
produces  normal  stresses  in  a  range  found  in  other  research,  which  produces  a 
coefficient  of  friction  in  this  range.  The  normal  stress  in  all  experiments  was 
between  7.03  and  17.63  ksi,  as  shown  in  Chapter  5.  The  magnitudes  of  stress  in 
each  experiment  are  high  and  agree  with  the  stresses  presented  in  other  research, 
which  produced  a  coefficient  of  friction  in  the  range  mentioned  above. 

Baltay  and  Gjelsvik  (1990)  conducted  a  study  that  dealt  with  the  coefficient 
of  friction  for  steel  on  concrete  at  high  normal  stress.  These  researchers 
concluded  that  the  average  coefficient  of  friction  between  the  machined  mild 
steel  and  concrete  was  found  to  be  0.47  for  stress  levels  between  1  and 
68,000  psi,  as  shown  in  Figure  3.7.  Also,  the  coefficient  of  friction  for  a  mill 
scale  steel  surface  is  less  than  that  for  the  machined  surface  for  stress  levels 
below  10,000  psi,  as  shown  in  Figure  3.8.  It  was  observed  that,  at  medium  stress 
levels  of  100  to  1,000  psi,  the  steel  surfaces  were  heavily  scratched  while  the 
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Figure  3.7.  Coefficients  of  friction  for  machined 

mild  steel  surface  (after  Figure  3  Baltay 
and  Gjelsvik  1990) 
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Figure  3.8.  Coefficients  of  friction  for  mild  steel  surface 
with  mill  scale  (after  Figure  4  of  Baltay  and 
Gjelsvik  1 990) 


concrete  was  (or  appeared  to  be)  unscathed.  Steel  particles  were  left  behind  on 
the  concrete. 

In  addition,  Baltay  and  Gjelsvik  concluded  that,  at  high  stress  levels,  above 
10,000  psi,  the  steel  gouged  into  the  concrete.  When  the  steel  was  removed, 
small  sand  particles  from  the  concrete  were  found  between  the  steel  and  concrete. 
Finally,  the  orientation  of  the  machining  cuts  on  the  steel  surface  had  no  effect  on 
the  coefficient  of  friction,  and  the  coefficients  of  friction  obtained  are  very  close 
to  values  obtained  by  other  authors  at  low  stress  levels  (Rabbat  and  Russell  1 985) 
and  in  classic  experiments  of  steel  on  stone. 
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Friction  is  a  phenomenon  of  plastic  surface  yielding  (Bowden  and 
Tabor  1950).  Between  two  different  materials,  friction  is  due  in  part  to  the 
deformation  of  one  of  the  materials.  With  steel  and  concrete,  it  appears  that  the 
steel  yields  to  the  harder  particles  in  the  concrete  mixture,  resulting  in  gouges  in 
the  steel.  However,  at  higher  stress  levels  as  the  concrete  begins  to  crush,  small 
particles  are  locally  tom  free  from  the  concrete,  and  the  steel  begins  to  plow  a 
path  through  the  concrete.  It  seems  that  the  value  of  the  coefficient  of  friction  is 
governed  by  the  steel,  that  is,  by  the  local  plowing  of  hard  particles  through  the 
steel.  This  was  also  confirmed  by  the  tests  with  mill  scale.  The  mill  scale  is 
harder  than  the  steel  and  is  therefore  not  penetrated  by  the  concrete  particles  at 
stress  levels,  resulting  in  a  lower  coefficient  of  friction.  As  shown  in  the  second 
row  of  Table  3.4,  the  range  of  values  for  the  coefficient  of  friction  between 
concrete  and  steel  is  around  0.3  and  0.58  for  stress  levels  lower  than  10,000  psi. 
These  stress  levels  were  similar  to  the  stress  levels  reached  in  the  barge  impact 
experiments  in  1998. 
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4  Equilibrium  Formulation 
with  Energy  Method 
Application 


In  this  chapter,  the  energy  formulation  is  introduced  in  a  formulation  to 
calculate  the  forces  applied  to  the  wall  due  to  the  impact.  This  method  uses 
Castigliano’s  theorem,  which  states  that  the  partial  derivative  of  the  total  strain 
energy  of  the  arc  beam  with  respect  to  one  external  load  is  equal  to  the 
displacement  in  the  direction  of  the  applied  external  load.  Our  case  is  a  first- 
degree  indeterminate  system,  and  one  of  the  support  reactions  was  selected  as  the 
redundant.  All  the  internal  forces  (axial,  shear,  and  moment)  were  expressed  in 
terms  of  the  redundant  force  located  in  cell  10  (or  equivalently,  pin  B).  The  first 
configuration,  mentioned  in  Chapter  1,  is  used  here  to  compare  the  results  with 
the  values  obtained  in  Chapter  3.  Referring  to  Figures  4.1  and  4.2,  the 
equilibrium  equations  are 


£fx  =0.'.  Ax+Bx-Sw  =0 

t©  ^Fy=0:.  Ay+By-Fw=  0 

A  =  0  Fw/?sin(ft)- -sin(e)]-B^/?[sin(e)-sin(a)] 
-  Rti?[cos(dj)  +  sin  (ft)]  =  0 


where 


(4.1) 

(4.2) 

(4.3) 


a  =  54  +  0 
ft  =  0  -  7 
c  =  90  +  ft 
d  =  90  -a 
e  =  90  -  ft 
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Structure 


Free  Body  Diagram 


Figure  4.1 .  Structure,  idealization,  and  free-body  diagram  of  the 
bumper — energy  method 
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Y 


c.  Forces  acting  on  load  bumper 
Figure  4.2.  (Concluded) 


Note  that  the  values  assigned  to  angles  a  through  e  assume  that  the 
orientation  of  the  arc  is  the  same  as  the  first  configuration  used  to  analyze  the 
system,  and  the  forces  SPA  (load  cell  1 1)  and  SPB  (load  cell  1 0)  are  oriented  in  the 
radial  direction  (Figure  4.2c).  In  addition  to  the  three  equilibrium  equations,  we 
can  introduce  the  axis  transformation  equations,  which  transform  the  local  to 
global  coordinates.  These  four  new  equations  are 


4*  =  FNA  cos(0-7 )-Spa  sin(0-7) 

(4.4) 

Ay  =  Fm  sin(0  -7 )+SPA  cos{0  -  l) 

(4.5) 

Bx  =  -Fnb  cos(36  -  0)  +  S PB  sin(36  -  0) 

(4.6) 

By  =  Fnb  sin(36  -  9)  +  S PB  cos(36 -  6 ) 

(4.7) 

The  additional  equation  is  the  displacement  expression  obtained  from 
Castigliano’s  theorem,  which  is  equal  to  zero  because  the  redundant  was  taken  in 
the  support.  The  general  expressions  of  this  equation  are 

cU  f  N  cN  f  V  , 

Sb*  =  ffix  =  J  EA  ffix  Jj+  J  GA  oBx  ds  + 

r  M 

———^  =  0 

J  El  cBx 

(4.8) 

ds  =  Rd6 

(4.9) 
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where  N,  V,  and  M  are  the  internal  normal,  shear  forces,  and  internal  moment 
expressed  in  terms  of  the  reaction  Bx.  The  arc  beam  was  divided  into  two 
segments— one  from  support  B  to  the  vertical  global  axis  (point  of  application  of 
Fw  and  Sw),  and  the  second  from  support  £  to  a  section  after  the  application  of  Fw 
and  Sw.  In  other  words,  the  limits  of  the  integrals  in  the  problem  are  from  a 
(a  =  5  deg  +  9)  to  90  deg  and  from  90  deg  to  c  (c  =  90  deg  -  7  deg). 

This  new  formulation  did  not  introduce  any  assumption  regarding  the  value 
of  the  coefficient  of  friction  and  did  not  consider  the  arc  beam  as  a  rigid  body, 
but  considered  it  as  a  flexible  body.  Solving  these  eight  equations,  the  problem  is 
transformed  to  a  determinate  system.  Now,  our  knowns  are  SPA  (load  cell  1 1)  and 
SpB  (load  cell  10),  which  are  the  radial  forces  measured  during  the  experiments. 
The  unknowns  are  the  reactions  in  global  coordinates  (Ax,  Ay,  Bx,  and  By),  the 
forces  in  the  wall  ( Fw  and  Sw),  and  the  two  reactions  normal  to  the  cross  section 
of  the  beam  in  local  coordinates  (Fma  and  F^B).  In  this  case,  a  Maple  worksheet, 
presented  in  Appendix  B,  was  created  to  calculate  numerically  the  values  for  Fw 
and  Sw. 

To  determine  if  the  results  provided  by  the  energy  model  formulations  are 
accurate,  two  finite  element  models  of  the  bumper  were  studied  using  the  com¬ 
puter  programs  SAP2000  and  VisualAnalysis.  The  two  models  were  constructed 
by  fixing  the  approach  angle  to  18  and  21.5  deg.  This  results  in  the  same  model 
structure  but  with  different  positions  in  the  global  coordinate  system.  Figure  4.3 
shows  the  arc  layout  used  for  the  finite  element  analyses  for  0  =  18  deg. 

Figure  4.4  shows  the  arc  layout  for  the  VisualAnalysis  model  for  0  =  18  deg,  and 
Figure  4.5  shows  the  model  for  0  =  18  deg  built  using  SAP2000.  The  properties 
of  the  beam  elements  used  in  these  models  are  as  follows: 

Cross-sectional  area  =  45  in.2 

Moment  of  inertia  =  303.75  in.4 

Modulus  of  elasticity  (£)  =  29,000  ksi 

E 

Shear  modulus  (G)  - - 

2(1  +  v) 

Poisson  ratio  (  V  )  =  0.33 

The  sectional  properties  were  calculated  using  the  dimensions  of  5  by  9  in.  In 
this  case,  the  same  approach  as  in  Chapter  3  was  used.  A  unit  load  Fw  and  a  load 
of  0.50  kip  for  Sw  was  applied  to  the  finite  element  models,  corresponding  to  a 
coefficient  of  friction  of  0.50  (shaded  cells  in  second  and  third  rows  of  Tables  4.1 
and  4.2).  Then,  the  values  of  the  supports’  reactions  were  computed  from  the 
finite  elements  models  using  both  computer  programs.  The  proposed  equilibrium 
formulation  with  energy  method  was  then  used  with  these  support  reactions  as 
input  data,  to  calculate  the  resulting  normal  and  shear  force  applied  to  the  arc 
beam.  Remember  that  the  problem  of  the  barge  impact  is  an  inverse  problem. 
That  is,  the  known  forces  are  two  reactions,  and  the  unknown  forces  are  the  other 
two  reactions  at  the  pins  and  the  forces  normal  and  parallel  to  the  wall  at  the 
point  of  contact.  It  is  a  common  practice  to  apply  a  traditional  procedure,  which 
states  that,  with  a  known  external  applied  load,  the  reactions  and  internal  forces 
are  calculated.  The  case  with  which  we  are  dealing  in  the  barge  impact  problem 
is  the  opposite  of  this  traditional  procedure. 
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Figure  4.3.  Idealization  of  the  bumper  structure 
for  the  finite  element  analysis 


Figure  4.4.  VisualAnalysis  model  for  0  =  1 8  deg — energy  method 


Figure  4.5.  SAP2000  model  for  0  =  1 8  deg— energy  method 
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Table  4.1 

Comparisons  of  Results,  Energy  Method— Approach  Angle  21.5  deg  (shaded  areas  are  input  data  for  each  model) 
Without  Shear  Term 
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The  comparisons  between  the  results  produced  by  each  model  are  presented 
in  Tables  4.1  and  4.2.  For  the  model  presented  in  this  chapter,  the  coefficient  of 
friction  is  calculated  from  the  ratio  of  Sw  divided  by  Fw.  The  results  obtained  with 
this  model  produce  a  coefficient  of  friction  in  the  range  that  was  found  in  the 
technical  literature  review  (e.g.,  Baltay  and  Gjelsvik  1990).  The  procedure  to 
obtain  these  results  was  as  follows  (data  presented  are  obtained  from  Table  4.1 , 
without  shear  term): 

a.  The  arc  was  fixed  using  a  specific  0  value.  Due  to  the  use  of  straight 
elements  to  model  the  arc,  the  angles  that  make  the  initial  element  at 
each  support  with  the  horizontal  (global  axis)  have  a  small  difference 
from  the  exact  values  (0-7  deg)  and  (36  deg  -  0)  for  pin  A  (load  cell  11) 
and  pin  B  (load  cell  10),  respectively.  These  values  are  specified  by  ©i 
and  02  in  the  last  two  rows  of  Tables  4.1  and  4.2.  Although  they  are 
different  angles,  the  average  value  results  in  the  exact  approach  angle,  0, 
used  in  the  proposed  method,  as  shown  in  the  first  row  in  Tables  4.1 
and  4.2.  For  example,  the  finite  element  models  (in  Table  4.1)  have 

(0!  +  02)/2  =  (21 .03  deg  +  21 .97  deg)/2  =  21 .5  deg. 

b.  Using  SAP2000  and  VisualAnalysis,  the  model  was  constructed. 

c.  The  Fw  =  1.0  and  a  coefficient  of  friction  (0.5)  were  applied  to  the 
structure  ( Fw  =1.0  kip,  Sw  =  0.5  kip,  and  //  =  Fw  /  Sw  =  0.5). 

d.  The  reactions  in  the  global  coordinates  were  calculated  {Ax  =  1.108  kip, 

Ay  =  0.532  kip,  Bx  =  -0.581  kip,  and  By  =  0.465  kip)  (Pin  A  =  load 

cell  11;  Pin  B  =  load  cell  10). 

e.  The  components  in  the  local  coordinates  were  obtained  {SPA  =  0.238  kip; 
Spb  =  0.305  kip). 

f  The  radial  forces  (SPA  and  SPB)  were  used  as  input  data  for  the 
equilibrium  formulation  with  energy  model. 

g.  The  F„  and  Sw  were  determinate  using  the  equilibrium  formulation  with 
energy  model  ( Fw  =  0.998  kip;  Sw  =  0.526  kip). 

h.  The  coefficient  of  friction  was  calculated  as  Sw/ Fw  (/i  =  FJ Sw  =  0.527). 

The  coefficient  of  friction  obtained  with  the  equilibrium  formulation  with 
energy  model,  after  following  this  procedure,  was  very  close  to  the  initial 
assumption  of  n  =  05 .  From  this  point  of  view,  Tables  4.1  and  4.2  show  the  good 
agreement  between  the  energy  model  and  the  finite  element  models. 

Recall  that  the  proposed  method  includes  the  shear,  axial,  and  bending 
effects  in  the  arc  beam,  as  shown  in  Equation  4.8.  The  finite  element  models  used 
in  these  analyses  consider  the  axial  and  bending  effects  only.  It  can  be  observed 
from  Tables  4.1  and  4.2  that  a  difference  occurs  in  the  results  when  the  shear 
term  is  considered  via  the  energy  method.  If  the  shear  term  is  not  considered,  as 
in  the  finite  elements  models  used  here,  the  Castigliano’s  theorem  takes  the 
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following  form,  which  is  the  formulation  used  to  obtain  the  results  presented  in 
Tables  4.1  and  4.2  without  shear  term. 


rNcN  rM  cM 
- ds+ - 

J  ea  mx  J  ei  mx 


ds-  0 


(4.10) 


The  procedure  to  obtain  the  Fw  and  Sw  when  the  shear  term  is  considered  is 
the  same  as  the  approach  presented  earlier  in  this  chapter.  For  example,  the 
results  using  the  energy  method  presented  in  Table  4.1  with  shear  term  are  as 
follows: 

a.  The  arc  was  fixed  using  a  specific  0  value.  Due  to  the  use  of  straight 
elements  to  model  the  arc,  the  angles  that  make  the  initial  element  at 
each  support  with  the  horizontal  (global  axis)  have  a  small  difference 
from  the  exact  values  (0  -  7  deg)  and  (36  deg  -  0)  for  pin  A  (load  cell  1 1 ) 
and  pin  B  (load  cell  10),  respectively.  These  values  are  specified  by  0i 
and  02  in  the  last  two  rows  of  Tables  4.1  and  4.2.  Although  they  are 
different  angles,  the  average  value  results  in  the  exact  approach  angle,  0, 
used  in  the  proposed  method,  as  shown  in  the  first  row  in  Tables  4.1  and 
4.2.  For  example,  the  finite  element  models  in  Table  4.1  have  (0,  +02)/2 
=  (2 1 .03  deg  +  2 1 .97  deg)/2  =  21.5  deg. 

b.  Using  SAP2000  and  VisualAnalysis,  the  model  was  constructed. 

c.  The  Fw  =  1 .0  and  a  coefficient  of  friction  of  0.5  were  applied  to  the 
structure  (Fw  =  1.0  kip,  Sw  =  0.5  kip,  and  //=  Fw/Sw  =  0.5). 

d.  The  reactions  in  the  global  coordinates  were  calculated  (Ax  =  1 .264  kip. 
Ay  —  0.573  kip,  Bx  =  -0.737  kip,  and  By  =  0.506  kip)  (Pin  A  =  load  cell 

1 1 ;  Pin  B  =  load  cell  1 0). 

e.  The  components  in  the  local  coordinates  were  obtained  ( SPA  =  0  238  kip- 
SpB  =  0.305  kip). 

f.  The  radial  forces  (Spa  and  Spg)  were  used  as  input  data  for  the 
equilibrium  formulation  with  energy  model. 

g.  The  Fw  and  Sw  were  determinate  using  the  equilibrium  formulation  with 
energy  model  ( Fw  =  1 .07  kip,  and  Sw  =  0.526  kip). 

h.  The  coefficient  of  friction  was  calculated  as  Sw/Fw  (//  =  FJSW  =  0.492). 


After  following  this  procedure,  the  coefficient  of  friction  obtained  with  the 
equilibrium  formulation  and  the  energy  method  including  the  shear  term  was 
very  close  to  the  initial  assumption  of  fi  =  05  .  From  this  point  of  view. 

Tables  4.1  and  4.2  show  the  good  agreement  between  the  energy  model  and  the 
finite  element  models  if  the  shear  term  is  or  is  not  included  in  the  formulation. 
Almost  the  same  coefficient  of  friction  is  obtained  if  the  shear  term  is  or  is  not 
included.  However,  the  force  normal  to  the  wall  is  not  the  same.  For  example,  the 
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value  of  Fw  is  7  percent  greater  when  the  shear  term  is  considered  than  the  result 
obtained  without  the  shear  effect,  as  shown  in  Table  4.1 .  This  validation  is 
presented  to  confirm  that  the  energy  formulation  agrees  with  the  results  obtained 
with  the  finite  elements  programs.  From  now  on,  all  the  analyses  done  with  the 
energy  formulation  will  include  the  shear  term  because  it  has  a  significant 
contribution  to  the  internal  energy  of  the  arc  beam.  That  is,  in  Chapter  5,  when 
the  energy  formulation  is  used,  the  three  internal  effects  (axial,  shear  and 
moment)  are  considered  in  the  formulation. 

In  the  solution  process  of  the  numerical  model  presented  in  this  chapter,  we 
found  two  special  cases  of  the  model  that  should  be  mentioned.  These  cases 
occur  if  the  load  is  applied  directly  over  the  supports.  That  is,  if  0  =  7  deg  or  0 
=  36  deg,  the  equilibrium  formulation  with  energy  model  does  not  work.  The 
reason  for  this  problem  is  that  one  of  the  equations  in  the  system  (energy 
equation)  has  to  have  the  loads  over  the  beam  to  produce  internal  forces  resulting 
from  internal  deflection.  If  the  external  load  is  applied  over  the  supports,  no 
internal  energy  is  developed  inside  the  arc  beam.  This  leads  us  to  a  system  of 
equations  with  more  unknowns  than  equations.  To  use  the  energy  model,  the  Fw 
and  Sw  loads  should  be  placed  between  the  supports  to  avoid  an  indeterminate 
system.  Suppose  that  the  Fw  and  Sw  loads  are  applied  over  the  support  A  (load 
cell  11),  which  means  that  0  =  7  deg,  taking  the  moment  at  support  A  (load 
cell  11)  produces  zero  reactions  at  B  (load  cell  10).  Then,  making  the  summation 
of  forces  in  the  vertical  direction,  Fw  is  equal  to  Ay.  Finally,  making  the 
summation  of  forces  in  the  horizontal  direction,  we  obtain  one  equation  and  two 
unknowns  ( Ax  and  SH).  The  energy  equation  cannot  be  used  because  the  internal 
forces  are  zero,  resulting  in  an  indeterminate  system. 
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5  Comparison  Between  the 
Assumed  Coefficient  of 
Friction  Method  and  the 
Energy  Method 


5.1  Interpretation  of  Results  from  the  Initial 
Configuration  Assumed  for  F10  and 

Several  cases  will  be  presented  in  this  chapter  that  compare  the  results 
between  the  assumed  coefficient  of  friction  model  and  the  energy  model.  The 
initial  configuration  (Figure  1 .8)  was  used  with  the  energy  method  to  determine 
the  forces  normal  and  parallel  to  the  wall.  This  configuration  assumed  the  arc 
beam  to  be  at  54  deg  with  the  longitudinal  axis  of  the  barges,  and  the  measured 
loads  were  oriented  in  the  radial  direction,  as  shown  in  Figure  5.1. 

With  the  exception  of  the  results  for  Experiment  43,  this  configuration 
produces  unreasonable  values  of  coefficient  of  friction  when  compared  with  the 
values  reported  in  the  technical  literature  (e.g..  Table  3.4).  The  values  for  the 
coefficient  of  friction  for  this  initial  model  range  from  a  low  of  0.59  to  a  high  of 
5.69,  with  a  coefficient  of  friction  well  above  unity  for  nearly  all  impact 
experiments.  These  are  clearly  well  above  the  0.3  to  0.58  values  reported  by 
Baltay  and  Gjelsvik  (1990).  Also,  in  some  experiments,  this  configuration 
produced  tension  in  the  arc.  Tension  is  observed  in  the  values  of  the  horizontal 
reactions.  As  shown  in  Figure  3.3c,  the  reaction  Bx  (load  cell  10)  has  to  be 
negative  to  produce  compression  in  the  arc.  For  example,  the  reaction  Bx  (load 
cell  10)  in  Experiments  29,  30,  31,  37,  39,  40,  and  41  is  positive,  indicating 
tension  in  the  arc.  Then,  the  assumed  coefficient  of  friction  method  was  used 
with  the  unreasonable  values  of  coefficient  of  friction  as  input  data  obtained  with 
the  energy  method  to  check  the  energy  methods  results.  The  comparisons  of 
these  results  are  presented  in  Table  5.1.  Good  agreement  between  the  methods 
can  be  observed;  however,  the  unreasonable  results  obtained  (high  values  of 
coefficient  of  friction  or  tension  in  the  arc)  must  be  evaluated. 
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b.  Arc  geometry  used  for  force  analysis 

Figure  5.1 .  Initial  configuration  and  recorded  load  direction  (Continued) 
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c.  Forces  acting  on  load  bumper 
Figure  5.1.  (Concluded) 


After  a  detailed  analysis,  three  possible  sources  of  error  in  the  recorded  data 
were  found.  The  proposed  two  methods  (Chapter  3  and  4)  use,  as  input  data,  the 
(1)  approach  angle  of  the  barges,  (2)  arc  geometry,  and  (3)  location  of  the  arc  and 
the  orientation  of  the  recorded  forces.  Each  one  of  these  possible  sources  of  error 
will  be  discussed. 

After  analyzing  the  approach  angle,  0,  it  was  concluded  that  it  has  the  lower 
probability  of  error  since  it  was  measured  using  GPS  technology.  From 
Table  A. 3,  we  observe  that  for  the  last  5  sec  prior  to  impact  in  each  experiment, 
the  values  for  the  approach  angles  are  about  the  same.  This  indicates  that  the 
approach  angle  was  almost  constant  before  impact  occurs.  In  addition,  the 
registered  approach  angle  obtained  during  the  experiments  was  corrected  to  align 
the  line  joining  the  two  global  positioning  systems  used  with  the  barge’s 
longitudinal  axis.  This  correction  was  made,  and  has  a  value  of  0.4468  deg  in  all 
experiments.  The  nearly  constant  approach  angle  values  obtained  during  the  last 
5  sec  before  impact  are  supported  by  the  small  values  computed  for  the 
coefficient  of  variation  (COV)  in  the  GPS  approach  angle  data  for  each 
experiment.  The  largest  COV  is  computed  equal  to  2  percent  for  Experiment  30. 
In  the  other  experiments,  this  coefficient  had  a  value  of  less  than  1  percent. 

The  second  possible  source  of  error  was  the  wrong  determination  of  the  arc 
geometry.  This  kind  of  error  has  also  a  low  probability  of  being  responsible  for 
the  unreasonable  results  obtained  in  the  original  configuration  studied.  The 
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Table  5.1 

Comparison  of  Models  Using  Radial  Forces  and  Arc  Layout  of  7  and  36  deg  -  Initial  Configuration 


■s 

E  c 

S  o 

1  E 

o 

o 

i 

X 

X 

CD 

B 

in 

CD 

cd 

B 

00 

CN 

cd 

CN 

cd 

CO 

in 

cd 

CO 

cd 

o 

in 

CN 

o 

X 

CD 

in 

X 

If) 

ID 

CD 

X 

in 

o 

CN 

CN 

2 

I! 

i 

a 

CD 

3 

IN 

CN 

CD 

i 

1 

00 

in 

3 

CN 

1 

o 

r-. 

8 

00 

in 

CO 

CN 

CD 

00 

m 

CN 

3 

X 

CD 

3 

CD 

X 

3 

2 

3 

o 

CD 

CN 

CN 

X 

CN 

CN 

8 

co 

CN 

X 

CD 

X 

CN 

to 

X 

o 

CN 

CO 

in 

8 

CO 

CN 

1 

B 

1 

i 

i 

8 

8 

CO 

CD 

CO 

CD 

CN 

m 

a 

CD 

CD 

CO 

o 

NT 

CN 

CO 

8 

NT 

CN 

O 

| 

3 

8 

NT 

X 

CD 

00 

CD 

CO 

CD 

8 

co 

in 

co 

in 

CO 

a 

a 

CO 

3 

a 

In 

a 

in 

1 

i 

CD 

CN 

o> 

B 

1 

CD 

CN 

cd 

<o 

8 

ed 

co 

m 

o 

CD 

8 

in 

00 

3 

3 

8 

co 

d 

m 

a 

CD 

CO 

X 

in 

co 

8 

X 

CO 

CN 

I 

X 

CO 

3 

CN 

a 

oo 

00 

X 

8 

CO 

8 

8 

i 

<6 

s 

in 

CO 

r- 

x 

CO 

<N 

CD 

8 

CO 

CN 

i 

m 

3 

CN 

1 

CO 

in 

a 

o> 

CD 

O 

3 

CN 

NT 

co 

3 

a 

NT 

CN 

CO 

CN 

3 

O 

in 

& 

co 

X 

CO 

X 

co 

X 

<0 

m 

t'- 

CN 

CO 

in 

CN 

CN 

a 

X 

f- 

o 

X 

8 

co 

r- 

X 

r-~ 

§ 

3 

CO 

X 

1 

$ 

s 

CD 

00 

IN 

CO 

CD 

CD 

CO 

05 

in 

r) 

3 

ON 

in 

3 

8 

8 

NT 

i 

i 

| 

1 

Sj 

o 

in 

05 

in 

m 

CN 

d 

CN 

o 

3 

05 

00 

in 

CD 

X 

_ 

m  . 

o 

id 

CN 

CD 

CN 

3 

NT 

CD 

in 

NT 

X 

a 

CN 

in 

in 

CN 

in 

o 

8 

X 

a 

i 

1 

CO 

© 

CN 

r~- 

8 

8 

CN 

8 

CD 

CD 

§ 

3 

CD 

CO 

i 

: 

3 

3 

X 

in 

1 

in 

1 

S 

£ 

g 

3 

8 

8 

in 

8 

a 

1 

i 

8 

§ 

o 

§ 

a 

cd 

NT 

CO 

cd 

CO 

CO 

i 

a 

3 

a 

1 

8 

CO 

X 

X 

X 

3 

X 

in 

S 

co 

in 

s 

2 

8 

8 

3 

CN 

CN 

cn 

a 

& 

8 

B 

B 

hr 

co 

8 

CO 

8 

2 

1 

3 

1 

CD 

NT 

in 

X 

in 

CN 

X. 

8 

8 

CN 

g 

CO 

2 

| 

*5 

i 

x 

s i 

s 

a 

CO 

x 

8 

8 

3 

CN 

ed 

CN 

a 

a 

8 

X 

CO 

V 

in 

r~ 

cd 

NT 

<D 

cd 

NT 

1 

a 

CD 

in 

a 

m 

cd 

CO 

3 

3 

a 

| 

3 

a 

3 

8 

3 

8 

3 

8 

1 

a 

(N 

O 

CO 

s 

co 

CD 

CD 

o 

CD 

x 

I 

CO 

■n 

CO 

a 

CO 

a 

$ 

a 

05 

8 

CO 

i 

1 

CO 

3 

cO 

i 

O 

m> 

in 

co 

cO 

co 

in 

in 

co 

co 

I 

i 

i 

£ 

X 

1 

i 

£ 

8 

co 

X 

8 

co 

8 

00 

CD 

N 

in 

<o 

1 

3 

CD 

8 

O 

CN 

00 

o 

a 

r-~ 

x 

8 

8 

X 

CO 

B 

CN 

1 

i 

in 

X 

3 

3 

co 

co 

1 

I 

1 

3 

i 

1 

co 

0 

x 

CD 

O 

a 

d 

3 

B 

B 

B 

B 

1 

8 

8 

s 

CN 

e 

n 

3 

00 

CD 

£ 

1 

1 

D 

8 

B 

B 

B 

X 

in 

8 

s 

E 

o 

E 

cn 

CD 

x 

in 

x 

CO 

3 

CO 

00 

IN 

CO 

i 

1 

N 

n 

1 

8 

cd 

CN 

B 

1 

8 

8 

8 

o 

2 

£ 

U 

8 

8 

CO 

CO 

00 

co 

X 

co 

o 

NT 

NT 

CN 

X 

co 

X 

3 

■o 

• 

H 

E 

+ 

E 

3 

X 

i 

& 

Ul 

1 

2 

£ 

i 

* 

X 

• 

X 

E 

♦ 

E 

3 

X 

.a 

£ 

3 

X 

Ul 

X 

o 

§ 

Z 

| 

c 

Ul 

* 

X 

• 

X 

u. 

+ 

E 

X 

JD 

3 

X 

Ul 

X 

o 

i 

Z 

£ 

€> 

c 

Ul 

n 

X 

• 

X 

E 

♦ 

E 

€ 

£ 

3 

X 

Ul 

X 

o 

f 

Z 

| 

c 

Ul 

X 

• 

E 

♦ 

E 

3 

X 

£ 

3 

X 

Ul 

X 

o 

£ 

z 

1 

c 

Ul 

* 

X 

• 

£ 

♦ 

E 

3 

X 

D 

□ 

X 

Ul 

X 

O 

£ 

z 

£ 

• 

c 

Ul 

m. 

X 

• 

£ 

♦ 

X 

.O 

3 

X 

Ul 

X 

o 

i 

Z 

£ 

• 

c 

Ul 

«. 

X 

• 

M 

E 

♦ 

E 

3 

X 

O 

3 

UJ 

X 

o 

i 

Z 

£ 

• 

c 

Ul 

X 

• 

M 

E 

+ 

x 

a 

□ 

X 

Ul 

X 

o 

1 

1 

«. 

X 

• 

£ 

+ 

X 

§ 

3 

& 

X 

o 

i 

Z 

£ 

l5 

« 

X 

• 

H 

E 

+ 

| 

X 

Z 

£ 

1 

54 
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reason  for  assigning  a  low  probability  is  because  the  arc  bumper  was  made  in  a 
machine  shop  with  specific  requirements  set  forth  by  shop  drawings.  These 
requirements  can  be  measured  after  the  beam  was  constructed.  The  two  intrinsic 
values  of  the  arc  are  the  central  angle  made  by  the  radial  directions,  which  is 
29  deg,  and  the  chord  between  the  cell  pins,  which  is  35.5  inches.  Both  items 
were  measured  before  and  after  the  experiments,  and  were  unchanged.  For  that 
reason,  this  possibility  was  determined  to  be  an  unlikely  source  of  error  in  the 
models. 

The  third  source  of  error  is  the  most  probable,  because  the  exact  position  of 
the  arc  on  the  barge  and  the  exact  load  cell  orientation  were  not  precisely 
measured  after  mounting,  nor  checked  after  each  experiment.  With  this  in  mind, 
two  additional  configurations  were  studied.  These  two  configurations  will  be 
discussed  in  detail  in  Sections  5.2  and  5.3. 

The  second  configuration  was  the  same  as  the  first  one,  but  with  a  change  in 
the  load  cell  orientation.  This  configuration  was  obtained  from  a  combination  of 
design  drawings  and  documentary  photos  taken  during  the  experiments.  The 
angle  obtained  from  the  drawings  and  photos  were  3  deg  from  the  radial  direction 
at  support  A  (load  cell  11)  and  6  deg  from  the  radial  direction  at  support  B  (load 
cell  10).  After  a  detailed  geometrical  analysis,  this  model  was  rejected  because 
unreasonable  values  of  coefficient  of  friction  were  obtained. 

The  third  configuration  was  obtained  after  performing  a  parametric  analysis. 
In  this  case,  the  values  of  several  geometrical  parameters  were  changed  until 
reasonable  values  of  coefficient  of  friction  and  compression  in  the  arc  were 
obtained.  A  final  geometrical  configuration  was  obtained,  and  it  was  very  similar 
to  the  first  and  second  configurations  studied.  The  resulting  value  for  the  location 
of  the  arc  was  57.5  deg  from  the  longitudinal  axis  of  the  barge  instead  of  54  deg 
as  used  in  the  previous  configurations.  In  this  third  configuration,  the  recorded 
forces  at  the  supports  are  oriented  5.5  and  1.5  deg  outside  the  radial  direction  at 
supports  A  (load  cell  11)  and  B  (load  cell  10),  respectively. 


5.2  Configuration  with  Recorded  F10  and  Fu 
Outside  the  Radial  Direction 

After  the  initial  configuration  was  analyzed,  and  unreasonable  index  values 
of  coefficient  of  friction  were  computed  for  all  the  experiments,  a  second 
configuration  was  developed  from  a  combination  of  design  drawings  and  docu¬ 
mentary  photos  taken  during  the  experiments.  This  second  assessment 
established  that  the  arc  bumper  was  placed  54  deg  from  the  longitudinal  axis  of 
the  barge.  In  addition,  the  recorded  forces  FI0  and  Fn  were  oriented  6  and  3  deg 
outside  the  radial  direction,  respectively.  This  configuration  is  presented  in 
Figure  5,2.  If  the  recorded  forces  are  outside  the  radial  direction,  a  new  set  of 
coordinate  transformation  has  to  be  included.  This  set  of  equations  has  the 
following  form: 
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where  F,0n  and  FUn  are  the  new  unknown  support  reactions  acting  perpendicular 
to  the  recorded  forces  F 10  and  Fu.  Angles  8  and  C,  are  those  that  the  recorded 
forces  Fu  and  F10  make  with  the  radial  direction.  With  this  modification  of  the 
formulation,  the  energy  method  was  applied  to  determine  the  normal  and  shear 
forces  in  the  wall  generated  during  impact.  Again,  the  model  produced 
unreasonable  results,  as  can  be  observed  in  Table  5.2. 

The  values  for  the  coefficient  of  friction  for  this  second  model  range  from  a 
low  of  -0.82  to  a  high  of  5.53,  with  a  coefficient  of  friction  above  unity  for  a 
majority  of  the  impact  experiments.  This  value  is  clearly  well  above  the  0.3  to 
0.58  values  reported  by  Baltay  and  Gjelsvik  (1990).  However,  in  this  second 
configuration,  the  arc  is  in  compression  in  all  the  experiments.  Compression  in 
the  arc  is  observed  in  the  values  of  the  horizontal  global  reactions.  As  shown  in 
Figure  3.3c,  the  reaction  Bx  (load  cell  10)  has  to  be  negative  to  produce 
compression  in  the  arc,  as  is  the  case  with  this  second  configuration,  as  shown  in 
Table  5.2. 

To  check  the  results  obtained  with  the  energy  method,  the  fixed  coefficient  of 
friction  method  was  employed  using  the  coefficient  of  friction  obtained  from  the 
energy  method.  The  same  results  were  obtained,  indicating  that  two  models 
produce  the  same  behavior.  That  is,  the  high  values  of  coefficient  of  friction  and, 
even  worse,  negative  Sw  force,  as  in  Experiments  43  and  44,  indicate  that  it  acts 
in  the  direction  of  the  motion,  which  is  impossible.  That  is,  this  second 
configuration  produces  a  negative  shear  force  S„,  or  equivalently  a  shear  force  in 
the  direction  of  the  motion,  according  to  Figure  3.3c. 

The  explanation  for  this  odd  behavior  was  found  after  a  detailed  analysis  of 
the  geometry  of  the  system.  First,  a  detailed  drawing  of  the  supports  was  made  to 
see  if  the  angles  adopted  from  the  combination  of  design  drawings  and 
documentary  photos  taken  during  the  experiments  follow  all  geometrical  rules. 
Figure  5.3  shows  the  triangle  that  forms  the  radial  axis,  the  load  direction,  and 
the  local  axis  of  the  barges.  For  support  B  where  F10  was  recorded,  no  problem 
was  found  because  the  summation  of  the  internal  angles  of  the  triangle  formed  by 
the  action  line  of  the  force  Fw,  the  longitudinal  axis  of  the  barge  (local  axis),  and 
the  radial  direction  is  equal  to  1 80  deg.  A  different  situation  occurs  at  support  A 
where  Fu  was  recorded. 

As  shown  in  Figure  5.3,  it  is  geometrically  impossible  that  the  summation  of 
the  internal  angles  of  a  triangle  be  greater  than  1 80  deg.  In  other  words,  if  the 
recorded  force  Fn  makes  an  angle  of  3  deg  with  the  perpendicular  to  the 
longitudinal  axis  of  the  barge,  then  the  location  of  the  arc  bumper  could  not  be 
54  deg  with  the  longitudinal  axis  of  the  barges.  That  is,  the  54  deg  with  the 
longitudinal  axis  produces  an  angle  of  7  deg  with  the  perpendicular  axis  because 
the  central  angle  of  the  arc  is  29  deg.  Then,  as  shown  in  Figure  5.3,  the  7  deg  plus 
the  8  angle  must  be  equal  to  the  external  3  deg,  which  is  impossible.  This  is  why 
the  second  configuration  produces  unreasonable  results. 
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y  Local  A»s 
V  NA 


Figure  5.3.  Support  reaction  geometry  layout:  (left)  support  A,  F1h  and  (right) 
support  B,  F10 


5.3  Configuration  with  New  Location  of  the  Arc 
Beam  and  Recorded  F10  and  Fu  Outside  the 
Radial  Direction 

After  the  analysis  of  the  second  proposed  configuration  was  concluded,  a 
new  configuration  for  the  arc  beam  was  proposed.  Recall  that  it  is  necessary  to 
change  the  arc  position  because  the  angle  that  is  made  with  the  radial  direction 
must  be  less  than  the  external  angle,  as  was  demonstrated  in  the  previous  section 
(Figure  5.3).  After  several  attempts,  a  new  configuration  was  derived  (shown  in 
Figure  5.4)  and  was  adopted  because  it  was  the  only  configuration  (of  the  three) 
that  produced  reasonable  results  in  terms  of  the  coefficient  of  friction  and 
compression  in  the  arc.  Also,  the  normal  stresses  developed  at  the  instant  of 
maximum  normal  force  ( Fw)max  are  in  the  range  established  in  other  research 
(Table  3.4).  Now  the  arc  bumper  is  located  at  57.5  deg  from  the  longitudinal  axis 
of  the  barges  and  the  support  reactions’  geometry  agrees  with  the  basic  geometric 
law  that  the  internal  angle  of  a  triangle  sum  1 80  deg,  as  shown  in  Figure  5.5.  In 
this  configuration,  the  coordinates’  transformations  (presented  in  Equation  5.1) 
were  also  used. 

To  check  the  energy  model,  the  fixed  coefficient  of  friction  model  was  used 
with  the  coefficient  of  friction  obtained  from  the  energy  model  calculation.  The 
values  presented  in  Table  5.3  indicate  the  good  agreement  between  both  models, 
and  the  reasonable  coefficient  of  friction  values  obtained  with  both  models.  For 
example,  in  Experiment  29,  the  resulting  Fw  and  Sw  are  practically  the  same  in  the 
two  models.  We  can  compare  the  286.78  and  286.63  kips  for  the  fixed  coefficient 
of  friction  and  the  energy  method,  respectively. 
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a.  Arc  geometry  used  for  force  analysis 


Fv 

V 

\  Global  Axis 


b.  Forces  acting  on  load  bumper 

Figure  5.4.  Proposed  (third)  load  cell  and  arc  location  layout 
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Table  5.3 

Comparison  of  Models  Using  Recorded  Forces  Outside  the  Radial  Direction  and  Arc  Layout  of  3.5  and  32.5  deg  -  Third 
Configuration 
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Figure  5.5.  Support  reaction  geometry  layout  of  proposed  model:  (left) 
support  A,  F1U  and  (right)  support  B,  F10 


In  terms  of  the  resulting  coefficient  of  friction,  with  the  exception  of  one 
impact  experiment  (No.  42),  the  obtained  values  for  this  third  model  are  within 
the  0.43  to  0.6  range.  These  values  are  consistent  within  the  Baltay  and  Gjelsvik 
(1990)  data.  Using  this  third  configuration,  the  arc  is  in  compression  for  all  the 
impact  experiments.  It  can  be  observed  in  column  Bx  in  Table  5.3.  As  mentioned, 
the  horizontal  reaction  Bx  has  to  be  negative  (as  shown  in  Figure  3.3c)  for  the  arc 
to  be  in  compression.  This  condition  is  also  satisfied  with  this  third  configuration. 

One  additional  column  is  introduced  in  Table  5.3,  which  is  the  normal  stress 
developed  at  the  instant  of  maximum  normal  force.  These  stresses  are  calculated 
by  dividing  the  normal  force  to  the  wall  Fw  by  the  contact  area  (5  x  7.25  in. 

=  36.25  sq  in.).  As  shown  in  Table  5.3,  the  range  of  maximum  normal  stresses 
developed  for  all  the  experiments  is  6.35  to  15.93  ksi.  For  example,  for 
Experiment  38,  the  value  of  Fw  is  230.31  kips,  resulting  in  a  maximum  normal 
stress  of  a  =  F^/Area  =  230.31/36.25  =  6.35  ksi.  These  stresses  are  consistent 
with  the  range  presented  by  the  Baltay  and  Gjelsvik  (1990)  data  shown  in 
Table  3.4.  As  a  result  of  these  careful  evaluations  of  the  data,  the  third 
configuration  is  deemed  to  be  an  accurate  model  for  this  indeterminate  bumper 
system.  Key  factors  leading  to  this  conclusion  include  the  fact  that  this  third 
configuration  results  in  the  arc  bumper  being  in  compression  for  all  impact 
experiments  as  well  as  the  consistency  between  the  resulting  coefficients  of 
friction  and  the  stresses  developed  with  those  values  measured  in  the  Baltay  and 
Gjelsvik  (1990)  experiments. 

The  maximum  normal  force  (Fw)max  is  equal  to  the  reaction  force  provided  by 
the  lock  wall  on  the  barge  flotilla  during  the  impact.  Neither  the  forces  F row,  F drag, 
and  F heiper  acting  on  the  barge  flotilla  (and  depicted  in  Figure  2.3)  nor  the 
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hydrodynamic  forces  (corresponding  to  the  hydrodynamic  added  masses)  were 
measured  during  the  impact  tests.  However,  the  Table  5.3  computed  values  of 
(Fv,)mo.x  reflect  the  effects  of  all  external  forces  identified  in  Figure  2.3  acting  on 
the  barge  flotilla,  as  well  as  the  effects  of  the  hydrodynamic  masses  and  the 
inertia  of  the  barge  flotilla  mass  decelerating  during  impact. 


5.4  ETL  1110-2-338  Engineering  Procedure 


ETL  1 1 1 0-2-338  provides  an  engineering  procedure  for  the  collision  of  a 
barge  flotilla  with  a  rigid  structure.  According  to  the  ETL  1 1 10-2-338 
engineering  procedure,  the  magnitude  of  the  impact  forces  generated  by  a 
particular  collision  event  is  dependent  on  the  mass  (including  hydrodynamic 
added  mass  of  the  barge  flotilla),  the  approach  velocity,  the  approach  angle,  the 
barge  flotilla  moment  of  inertia,  damage  sustained  by  the  barge  structure,  and 
friction  between  the  barge  and  the  wall.  A  major  distinction  between  this 
procedure  and  the  traditional  Navy  method  for  determining  berthing  forces  is  the 
estimation  of  collision  energy  dissipated  in  deformation  of  the  barge  structure 
and  transferred  to  the  rotation  of  the  barge  flotilla.  The  analytical  method  uses  the 
structural  interaction  mechanism  of  Minorsky,  which  provides  an  empirical 
relationship  between  the  (nonrecoverable)  hull  deformation  and  the  energy 
absorbed  in  a  collision.  The  relationship  between  kinetic  energy  lost  in  a  collision 
and  the  volume  of  in-plane  (barge)  material  damaged  is  used  to  determine  impact 
force  as  a  relationship  to  instantaneous  contact  area  of  damaged  structure. 

Minorsky  used  the  conservation  laws  of  momentum  and  energy  and  the 
principles  of  rigid  body  mechanics  to  estimate  the  kinetic  energy  lost  during  a 
collision  between  two  vessels.  He  then  calculated  a  resistance  factor  that  is 
essentially  the  volume  of  material  damaged  in  the  bow  of  the  striking  ship  and  in 
the  side  of  the  struck  ship.  Minorsky  reasoned  that  the  principal  resistance  to 
collision  penetration  is  provided  by  deep  structure  that  suffers  in-plane  damage. 
For  the  case  of  a  barge  striking  a  fixed  wall,  the  main  deck,  the  bottom  plate,  the 
head  log,  and  the  transverse  frames  would  offer  resistance  to  damage. 

Minorsky  selected  and  analyzed  26  actual  ship  collisions  and  correlated  the 
energy  absorbed  in  the  collision  with  the  Minorsky  resistance  factor.  Using  the 
equivalency  between  energy  absorbed  and  the  work  performed  in  deforming  the 
structure,  a  constant  described  as  the  force  per  unit  of  damaged  surface  area  was 
defined  (=  13.7  ksi).  The  Minorsky  structural  interaction  mechanism  is  a  constant 
pressure  process  operating  with  a  pressure  of  13.7  ksi  over  the  instantaneous  face 
area  of  the  damaged  element.  This  allows  for  the  definition  of  an  equivalent, 
linear  spring  constant  representing  the  crushing  of  the  barge  structure  in  the 
ETL  1 110-2-338  analytical  formulation.  It  is  important  to  note  that  the  entire 
structural  interaction  mechanism  is  modeled  as  a  linear  spring  in  the  direction  of 
collision  corresponding  to  the  energy  absorption  in  the  crushed  barge  structure. 
The  formulation  becomes  one  of  an  initial  value  problem  for  a  barge  flotilla 
collision  with  a  rigid  wall,  representing  a  lock  wall  in  this  case,  and  leads  to  the 
solution  given  in  ETL  1110-2-338. 
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When  computing  values  of  maximum  impact  force  normal  to  the  wall, 
(Fy^max,  following  the  engineering  procedure  outlined  in  ETL-1 1 10-2-338,  the 
mass,  M,  and  moment  of  inertia,  I,  of  the  barge  flotilla  are  increased  to  account 
for  hydrodynamic  effects  by  means  of  hydrodynamic  added-mass  terms. 

Different  values  for  these  hydrodynamic  added  masses  are  used  in  different  barge 
flotilla  directions.  In  the  longitudinal  direction  of  the  barge  flotilla,  Mx  is  set 
equal  to  1 .05  M;  in  the  transverse  direction  My  is  set  equal  to  1  A  M,  and  7e  is  set 
equal  to  1 .47. 

The  ETL  1110-2-338  engineering  procedure  was  used  to  compute  values  of 
maximum  impact  force  normal  to  the  wall,  {Fw)max,  for  the  eight  full-scale, 
controlled  barge  impact  experiments  listed  in  Table  5.3.  The  coefficients  of 
friction  used  in  these  computations  were  the  same  as  those  derived  from  field-test 
data  summarized  in  Table  5.3.  The  results  of  these  ETL-based  computations  are 
given  in  Table  5.4.  Values  of  (Fy,)^  are  between  263.8  and  423.8  kips  for  these 
computations. 


Table  5.4 

Values  of  ( Fw)max ,  Time  of  End  of  Collision,  and  Deformation  of  Barge  at  Contact  Point  at 
End  of  Collision  by  the  ETL  1110-2-338  Engineering  Procedure 

Experiment 

Number 

Approach 

Velocity 

fps 

») »]  n 

9 

■ 

Coefficient  of 
Friction  - 
Table  5.3 

Time  of  End 
of  Collision 
sec 

Deformation  of 
Barge  at  Contact 
Point  at  End  of 
Collision 
ft 

CHS i 

29 

2.2 

12.63 

0.48 

0.6 

1.464 

0.619 

409.5 

30 

2.35 

12.19 

0.50 

0.48 

1.412 

0.619 

420.9 

31 

1.61 

10.6 

0.30 

0.43 

1.299 

0.342 

263.8 

37 

1.95 

10.29 

0.35 

0.52 

1.292 

0.399 

317.4 

38 

1.83 

11.94 

0.38 

0.57 

1.412 

0.472 

327.6 

39 

1.61 

14.12 

0.39 

0.51 

1.54 

0.532 

317.2 

41 

2.86 

8.76 

0.44 

0.51 

1.182 

0.458 

423.8 

42 

1.83 

17.48 

0.55 

0.09 

1.625 

0.793 

386.5 

*  Note:  (Fw)max computed  following  the  engineering  procedure  outlined  in  ETL-1 11 0-2-338,  which  uses  Mx  =  1.05  /Win  the 
longitudinal  direction,  My  =  1 .4  M  in  the  transverse  direction,  and  4  =  1 .47,  where  M  is  the  mass  and  /  is  the  moment  of  inertia  of 
the  flotilla. 

The  values  of  (Fw)max  given  in  Table  5.4  for  the  eight  controlled-impact 
experiments  occur  at  the  times  denoted  as  “time  of  end  of  collision.”  Note  that 
computed  values  for  the  time  of  end  of  collision  are  between  1.182  and  1.625  sec 
by  the  ETL  1 110-2-338  engineering  procedure,  while  the  actual  times  to 
maximum  force  recorded  for  these  eight  controlled-impact  experiments  were 
between  0.1282  and  0.2352  sec  (Table  A.2).  The  ETL  1 1 10-2-338  procedure 
overestimates  the  time  to  maximum  force  by  about  a  factor  of  8. 

Table  5.4  also  lists  the  computed  deformation  of  barge  at  contact  point  (i.e., 
impact  point)  at  end  of  collision  by  the  ETL-1 110-2-338  engineering  procedure. 
Values  range  from  0.342  to  0.793  ft  in  these  eight  computations.  This  computed 
deformation  is  the  deformation  of  the  linear  spring  in  the  direction  of  collision 
corresponding  to  the  energy  absorption  in  the  crushed  barge  structure.  Because 
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the  ETL  1 110-2-338  analytical  method  uses  the  structural  interaction  mechanism 
of  Minorsky,  which  provides  an  empirical  relationship  between  the 
(nonrecoverable)  hull  deformation  and  the  energy  absorbed  in  a  collision,  this 
deformation  may  be  viewed  as  an  indicator  of  the  nonrecoverable  damage  to  the 
barge  comer  as  a  result  of  the  barge  flotilla  impacting  the  lock  wall.  Note  that  no 
damage  was  observed  to  the  barge  comer  during  any  of  these  low-velocity, 
controlled-impact  experiments  at  Robert  C.  Byrd  Lock. 

The  results  of  these  eight  ETL  1 1 10-2-338  computations  are  repeated  in 
Table  5.5,  along  with  the  Table  5.3  values  for  ( Fw)max .  For  example,  for 
Experiment  29,  the  ETL  1 110-2-338  engineering  procedure  produces  an  (Fw)nmx 
value  of  409.5  kips,  and  the  Table  5.3  energy  method  or  the  equilibrium  with  a 
fixed  coefficient  of  friction  value  produces  an  (/vW  value  of  287.53  kips.  This 
corresponds  to  a  42-percent  overprediction  by  the  ETL  1 110-2-338  approach. 
The  maximum  force  normal  to  the  wall,  ( Fw)niax ,  computed  using  the  ETL  1110- 
2-338  engineering  procedure  differs  from  the  maximum  force  {Fw)nmx  values 
presented  in  Table  5.3  by  a  33-percent  underprediction  to  a  42-percent 
overestimate  for  these  eight  impact  tests. 


Table  5.5 

Comparist 

Engineerii 

an  Between  Table  5.3  ( Fw)max  Values  and  (Fw)max  by  the  ETL  1110-2-338 
ig  Procedure 

Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal 
to  Wall 
fps 

Coefficient 
of  Friction  - 
Table  5.3 

( Fw)  max 

Table  5.3 
kips 

(Fw)mtlx  ETL 
1110-2-338* 
kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.6 

287.53 

409.5 

42 

30 

2.35 

12.19 

0.50 

0.48 

371.16 

420.9 

13 

31 

1.61 

10.6 

0.30 

0.43 

235.30 

263.8 

12 

37 

1.95 

10.29 

0.35 

0.52 

327.49 

317.4 

-3 

38 

1.83 

11.94 

0.38 

0.57 

230.56 

327.6 

42 

39 

1.61 

14.12 

0.39 

0.51 

272.06 

317.2 

17 

41 

2.86 

8.76 

0.44 

0.51 

419.35 

423.8 

1 

42 

1.83 

17.48 

0.55 

0.09 

577.46 

386,5 

-33 

*  Note:  (Fw)max  computed  following  the  engineering  procedure  outlined  in  ETL-1 11 0-2-338,  which  uses  Mx 
the  longitudinal  direction,  My=  1.4  Min  the  transverse  direction,  and  l0  =  1.47,  where  Mis  the  mass  and 
moment  of  inertia  of  the  flotilla. 

=  1.05  Min 
is  the 

5.5  Fw  Time-Histories 

Using  the  third  model  configuration  that  was  described  in  Section  5.3,  a 
complete  time-history  of  force  normal  to  the  wall,  Fw,  during  the  approximately 
9  sec  of  impact  was  developed  for  eight  of  the  impact  experiments.  The  value  of 
the  coefficient  of  friction  used  in  this  calculation  was  obtained  in  Chapter  5  with 
the  energy  method  for  the  last  configuration  studied.  Using  a  computer  program 
presented  in  Appendix  B,  the  time-histories  for  eight  experiments  were 
calculated.  It  is  important  to  mention  that  three  of  the  experiments  had  some 
discrepancies  in  the  recorded  F10  and  Fu. 
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These  experiments  are  Nos.  40, 43,  and  44.  The  reason  for  neglecting  these 
experiments  is  that,  for  some  unknown  reason,  the  time-histories  for  Fjo  and  Fu 
cross  over  several  times,  or  at  some  instant  of  time  they  have  the  same  value  and 
in  another  instant  they  have  very  different  values  during  the  9  sec  of  the  impact 
process.  For  example,  Figure  A.15  presents  the  force  time-history  for 
Experiment  40.  From  this,  it  easy  to  see  that  the  maximum  force  was  recorded  at 
the  front  cell  (cell  10)  and  then,  from  the  third  peak,  the  maximum  load  was 
recorded  at  the  rear  cell  (cell  11). 

For  Experiments  43  and  44,  the  story  is  different.  We  can  easily  notice  from 
Figures  A.18  and  A.19  that,  at  the  beginning  of  the  time-history,  the  recorded 
force  at  both  load  cells  is  the  same  (that  is,  they  have  the  same  slope  for 
approximately  the  first  0.04  sec),  and  then  they  have  a  great  difference  in 
magnitude.  This  is  a  situation  that  can  be  explained  if  different  barge-to-wall 
angles  are  adopted  during  the  9  sec  of  impact,  which  was  checked  with  the  GPS 
data  recorded  during  the  entire  9  sec  of  the  impact  process.  We  concluded  that 
these  angles  remain  constant  during  the  impact  process  after  a  detailed  analysis 
of  the  GPS  data  of  these  experiments.  Experiments  40, 43,  and  44  were  the  last 
experiments  performed  during  the  full-scale,  low-velocity,  controlled  barge 
impact  experiments  at  Robert  C.  Byrd  Lock  in  December  1998.  Of  all  the 
experiments,  they  are  associated  with  the  highest  impact  angles  (17.52,  21.16, 
and  20.92,  respectively). 

We  can  observe  from  Table  5.3  that  the  only  experiment  for  which  the 
coefficient  of  friction  did  not  agree  with  the  values  reported  by  Baltay  and 
Gjelsvik  (1990)  is  Experiment  42,  which  produced  a  coefficient  of  friction 
of  0.09.  The  approach  angle  of  Experiment  42  was  17.48  deg,  which  is  very  close 
to  the  approach  angle  of  Experiment  40.  In  addition,  it  is  possible  that  the 
concrete  wall  surface  had  more  irregularities  (damage)  than  in  the  previous 
experiments  due  to  the  repeated  impacts  at  the  same  unarmored  concrete  wall 
location  by  the  flotilla.  It  is  important  to  notice  that  the  experiments  with  high 
approach  angles  produced,  in  some  way,  unreasonable  values  (that  is, 
unreasonable  coefficient  of  friction,  crossover  of  the  recorded  forces,  or  a 
variable  value  (not  constant)  of  the  difference  between  the  F/0  and  Fu  forces 
during  the  time-history).  The  high  approach  angles  associated  with  these 
experiments  and  the  fact  that  they  were  the  last  experiments  conducted  (and 
probably  were  impacting  a  damage  zone  in  the  concrete  wall  resulting  from  the 
previous  impact  experiments)  could  explain  the  unreasonable  coefficient  of 
friction  or  the  crossover  of  the  recorded  forces  in  the  time-history. 

As  shown  in  Table  5.6,  it  can  be  observed  that,  from  the  calculation  of  the 
time-histories  for  the  acceptable  eight  experiments  (Figures  5.6-5.13),  the 
resulting  maximum  force  normal  to  the  wall  Fw  is  not  equal  to  the  sum  of  the 
maximum  recorded  forces  F10  and  Fu.  It  is  not  true  that  the  Fw  has  to  be  equal  to 
F io  +  Fu,  because  the  forces  are  not  oriented  in  the  same  direction.  Force  Fw  is 
oriented  in  the  vertical  global  direction,  and  Fw  and  Fn  are  oriented  with  some 
angle  from  the  vertical  global  direction.  In  addition,  the  shear  term  affects  the 
distribution  of  forces  in  the  arc. 


66 


Chapter  5  Comparison  Between  the  Friction  Method  and  the  Energy  Method 


Table  5.6 

Comparison  of  the  Maximum  Recorded  Load  Cells  and  Maximum 
Force  Normal  to  the  Wall 


Komm  u 

FIOmax  +  Fllmax 
(kips) 

Fw 

(kips) 

29 

59,62 

172.68 

232.30 

286.78 

30 

94,39 

245.95 

340.34 

369.82 

31 

32.86 

161.68 

194.54 

236.32 

37 

31.66 

220.67 

252.33 

327.32 

38 

32.59 

138.24 

170.83 

230.31 

30 

84.91 

158.77 

243.69 

271,98 

41 

28.92 

315.49 

344.41 

419.38 

42 

153.14 

211.17 

17.48 

364.31 

577.44 
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Figure  5.7.  F10,  F1U  and  Fw  time-histories  for  Experiment  30 
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Figure  5.8.  F10,  F1h  and  Fw  time-histories  for  Experiment  31 
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| - F10  - F11  - Fw  j 


Figure  5.9.  F10,  F1U  and  Fw  time-histories  for  Experiment  37 
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Figure  5.10,  F10,  F1h  and  Fw  time-histories  for  Experiment  38 
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Load  (kips) 


Figure  5.1 1 .  F10,  F11t  and  Fw  time-histories  for  Experiment  39 
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6  Scaled  Impulsive  Force 
Formulation  and  Empirical 
Correlation 


6.1  Scaled  Impulsive  Force  Formulation 


The  impulse  and  linear  momentum  equation  normal  to  the  wall  can  be 
expressed  as 


1 2 

where  \Fwdt  is  the  total  impulse  due  to  the  force  normal  to  the  wall  or,  in  other 

words,  the  total  area  under  the  Fw  versus  time  of  impact  curve.  Recall  that  the 
velocity  normal  to  the  wall,  vsinO ,  is  assumed  to  be  zero  during  impact.  The 
maximum  normal  force  occurs  at  a  tt  from  the  beginning  of  the  impact  time- 
history.  The  total  impulse  of  the  force  normal  to  the  wall  can  be  divided  into  two 
impulses,  one  up  to  the  maximum  Fw,  and  a  second  calculated  to  the  end  of  the 
time-history  of  impact  at  time  f2,  as  shown  in  Figure  6.1 .  These  two  impulses  are 
designated,  as  shown  in  Figure  6.1,  as  areas  A /  and  A2,  respectively. 

The  experimental  total  impulse  equals  ( A\  +  A2).  Now,  Equation  6.1  becomes 

h  .  h 

j Fwdt  =  A,  +  A2  =  [m flgfuia  - mwaler  jv,  sin G  +  \Ft0W sin Gdt 

*'  '■  (6.2) 

h  h 

+  \F helper  cos  Gdt  -  \Fdrag  sin  Gdt 

h  h 


where 


h+&i 

A\  =  J  Fwalldt  =  f[Fwa„  Aq  (6.3) 

h 

and  (3  is  a  constant.  By  definition,  A2  is  equal  to  04/  +  A2 )  minus  A]t  or 
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(6.4) 


j^wall^  P inwall  Lax  A^i 

<1 


A  constant,  y,  is  introduced,  which  is  defined  as  the  ratio  of  the  impulse  up  to 
the  maximum  normal  force  Fw  (A;),  divided  by  the  experimental  total  impulse 
(A/  +  A2),  written  as 


y  =  a-T4 

Then,  Equation  6.3  can  be  written  in  the  form 


Ay  =/*(A,  +  A2) 
and  Equation  6.4  results  in 

4  =  0-7)*  (A +  4) 

Substituting  Equation  6.6  into  Equation  6.3,  we  obtain 

=y(a+a) 


(6.5) 


(6.6) 


(6.7) 


(6.8) 


Finally,  the  maximum  force  normal  to  the  wall  can  be  determined  from  the 
following  expression: 


L 


y(a+a) 

PAt, 


(6.9) 


Note  that  the  last  three  terms  on  the  right-hand  side  of  Equation  6.2  were  not 
measured  during  the  experiments.  However,  we  can  associate  the  linear 
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momentum  with  the  total  impulse  at  the  left  side  of  the  same  equation.  The  linear 
momentum  can  be  calculated  for  each  of  the  experiments  and  also  the  total 
impulse  normal  to  the  wall,  which  is  the  summation  of  A,  and  A2.  A  new  factor, 
which  is  the  ratio  of  the  total  impulse  due  to  the  normal  force  divided  by  the 
linear  momentum,  is  now  introduced  as 

FRAC  =  ^1Al  (6.10) 

mv  sin  0 

then,  Equation  6.9  transforms  to 

(p  )  -  T*  FRAC  *  mv  sin  6  (611) 

V  wJmax  pAtj 

This  approach  was  used  to  calculate  the  maximum  normal  force  in  the  wall 
using  the  relation  between  the  Fw  and  the  linear  momentum.  For  all  the 
experiments  presented  in  Table  6.1,  a  statistical  evaluation  of  all  these  factors  (J3, 
y,  FRAC,  etc.)  was  done.  These  factors,  p,  y,  FRAC,  and  At,,  are  dependent 
variables  of  the  independent  variable  (FP),^.  Then  using  the  average  values  of 
these  factors,  an  estimate  of  the  force  normal  to  the  wall  was  calculated.  These 
average  values  are  presented  in  Table  6.1a. 

As  an  example,  Experiment  31  was  selected.  Recall  that  the  total  mass  of  the 
flotilla  equals  1,865.59  k-sec2/ft  (see  Appendix  A).  Note  that  the  masses  used  to 
develop  the  linear  momentum  normal  to  the  wall  correlation  with  values  of 
( Fw)max  in  this  chapter  do  not  include  the  calculation  of  hydrodynamic  added 
masses.  For  these  experiments,  the  following  values  were  obtained,  as  shown  in 
Table  6.1a. 

Average  values: 

At,  =0.16  sec 

y  =  0.06279 
p  =  0.72224 
FRAC  =  0.8 


For  Experiment  31:  [The  area  under  the  force  time-history  A,  +  A2  was 
calculated  using  the  numerical  integration  scheme  presented  in  Appendix  C.] 

v  =  1.61  ft  /  sec 
0  =  10.6° 

A1  +  A2  =  466.83  k  -sec  (see  center  shaded  column  in  Table  6.1a) 
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Table  6.1 

(Fw)max  Scaled  Method 


Chapters  Sealed 


Computation: 


mv  sin  9  = 


^1,865 .59  k  SCC  ^ 
V  ft  j 


*  1.61- 


sec  J 


*  sin  (10.60  deg  )=  552.5162  k-s 


The  total  impulse  normal  to  the  wall,  the  (A,  +  A2)  shaded  column  in 
Table  6.1a,  due  to  the  force  normal  to  the  wall  is  less  than  the  linear  momentum 
normal  to  the  wall. 

For  Experiment  31,  the  total  impulse  normal  to  the  wall  is  computed  to  be 
84  percent  of  the  linear  momentum  normal  to  the  wall  (=  466.83  /  552.5162 
in  decimal  fraction).  This  ratio  is  identified  as  FRAC  in  Table  6.1a.  The 
remaining  16  percent  of  the  linear  momentum  is  provided  to  the  system  by  the 
drag  force,  the  helper  boat,  and  the  tow  boat,  as  shown  in  Figure  2.3.  For  all  eight 
experiments,  the  Table  6.1a  ratio  FRAC  establishes  that  the  total  impulse  normal 
to  the  wall  is  between  65  and  85  percent  of  the  linear  momentum  normal  to  the 
wall. 


Computation: 


y  *  FRAC  *  mv  sin  6  (0.06279)  (0.8)  (552.5 162) 
pAt{  ~  (0.72224)  (0.16) 


240.17 


kips 


Table  6.1a  reports  the  average  values  of  the  At],  y,  p,  FRAC,  and  impact 
experiment-specific  value  (A/  +A2)  used  to  calculate  the  maximum  Fw  value  by 
the  scaled  method,  using  Equation  6.1 1.  Applying  Equation  6.1 1  using  the 
average  values  for  At,,  y,  p,  and  FRAC  given  in  Table  6.1a,  the  estimate  of 
(Fw)max  for  Experiment  31  is  computed  equal  to  240.17  kips.  This  value  can  be 
compared  with  the  value  found  using  the  assumed  coefficient  of  friction  model  or 
the  energy  model  presented  in  Table  5.3  and  repeated  in  the  eighth  column  in 
Table  6.1a,  equal  to  236.20  kips.  The  approximation  of  the  maximum  force 
normal  to  the  wall  for  this  experiment  is  within  2  percent  of  the  Table  5.3  value. 

It  is  important  to  mention  that  if  the  actual  values  of  At,,  y,  p,  and  FRAC  are  used 
to  calculate  the  maximum  force  normal  to  the  wall,  the  values  obtained  for  the 
maximum  Fw  are  equal  to  the  maximum  Fw  reported  in  Table  5.3.  This  is  because 
the  maximum  Fw  reported  in  Table  5.3  is  used  as  the  independent  variable  to 
determine  the  factors  At,,y  /?,  apd  FRAC,  making  these  factors  dependent 
variables. 


The  resulting  values  of  the  maximum  force  normal  to  the  wall  using  the 
average  factors  ( At,,y ;  /?,  FRAC),  appear  in  Table  6.1b,  and  the  trend  of  the  scaled 
method  is  presented  in  Figure  6.2.  It  can  be  observed  that  the  trend  of  the  data 
points  is  to  an  increase  in  the  magnitude  of  the  force  normal  to  the  wall  with  an 
increase  in  the  linear  momentum  normal  to  the  wall.  This  trend  can  be  observed 
from  Equation  6.1 1  (middle  line  in  Figure  6.2),  which  was  used  to  relate  the 
linear  momentum  of  the  flotilla  and  the  maximum  force  normal  to  the  wall.  In 
addition,  a  least  squares  regression  procedure  applied  to  the  four  variables  {At,,  y, 
p,  FRAC)  was  used  to  develop  the  best-fit  straight  line  through  the  eight  data 
points  (for  the  eight  impact  experiments).  The  line  was  assumed  to  start  at  the 
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origin  (i.e.,  no  intercept  term  was  used  for  the  linear  regression).  The  average 
minus  one  standard  error  and  average  plus  one  standard  error  lines  are  also 
shown  in  Figure  6.2.  The  eight  data  points  shown  in  Figure  6.2  are  the  maximum 
force  normal  to  the  wall  obtained  from  Table  5.3  for  each  experiment.  These  data 
points  are  presented  only  for  comparison  purposes,  and  because  they  are  the 
independent  variables  from  which  the  factors  At,,  y,  /},  and  FRAC  were  obtained 
and  then  used  in  Equation  6.1 1.  The  statistical  evaluation  applied  to  the  four 
variables  (At,,  y,  j3,  FRAC)  was  done  using  the  computer  program  presented  in 
Appendix  B  and  verified  using  the  standard  Excel  statistical  package. 


Figure  6.2  is  a  correlation  based  on  low-velocity,  shallow-impact  (up  to 
21.1  deg)  experiments  that,  by  definition,  do  not  account  for  factors  that  manifest 
themselves  at  higher  velocities.  Additionally,  no  damage  occurred  to  the  flotilla 
of  barges,  and  no  lashings  broke  during  these  eight  impact  experiments. 
Therefore,  at  this  time  and  pending  the  results  of  additional  research.  Figure  6.2 
is  deemed  applicable  to  a  barge  flotilla  that  has  a  velocity  normal  to  the  wall  up 
to  and  not  exceeding  0.57  fps  (0.39  mph)  with  no  damage  occurring  during 
impact  events,  and  for  impact  angles  up  to  21.1  deg,  for  a  barge  flotilla  with  a 
linear  momentum  normal  to  the  wall  between  649.84  and  1,025.48  k-sec.  Finally, 
Equation  6.11  transforms  to 

(Fw  )max  =  0.43469  *  (mv  sin  6 ) 

using  the  average  At,,  y,  /3,  and  FRAC  (given  in  Table  6.1a  and  shown  in  Figure 
6.2),  where  the  units  for  the  mass,  velocity,  and  maximum  Fw  are  k-see2/ft,  ft/sec, 
and  kips,  respectively. 

The  maximum  normal  force  (Fw)nmx  by  the  scaled  method  is  equal  to  the 
reaction  force  provided  by  the  lock  wall  on  the  barge  flotilla  during  the  impact. 
Note  in  Table  6.1b  that  the  mass  used  to  develop  Figure  6.2  correlation  of  linear 
momentum  normal  to  the  wall  with  values  of  (F*,)^  is  the  mass  of  the  barge 
flotilla  (given  in  Appendix  A)  and  does  not  include  the  computation  of  any 
hydrodynamic  added  masses.  However,  the  Table  6.1a  values  of  (Fw)nKa 
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(computed  using  Equation  6.1 1)  that  are  used  in  the  Figure  6.2  correlation  reflect 
the  effects  of  the  Figure  2.3  external  forces  of  Ftom  Fdrag,  and  F helper  acting  on  the 
barge  flotilla,  as  well  as  the  effects  of  the  hydrodynamic  masses  and  the  inertia  of 
the  barge  flotilla  mass  decelerating  during  impact. 

A  summary  of  the  values  of  the  Table  5.3  maximum  impact  force  normal  to 
the  wall,  (Fw)max,  for  the  energy  method  or  the  equilibrium  with  a  fixed 
coefficient  of  friction  values  for  the  eight  full-scale,  controlled  barge  impact 
experiments  and  of  the  Figure  6.2  best-fit  straight-line  (JFM)max  values  is  given  in 
Table  6.2.  For  example,  the  energy  method  or  the  equilibrium  with  a  fixed 
coefficient  of  friction  value  (=  0.6)  for  the  full-scale,  controlled  barge  impact 
experiment  (No.  29)  produces  an  (FK.)mM  value  of  286.63  kips,  and  the  best-fit 
straight  line  of  the  scaled  impulsive  relationship  produces  an  ( F 'w)mx  value  of 
390.09  kips.  This  corresponds  to  a  35-percent  overprediction  by  the  best-fit 
straight  line  of  the  scaled  impulsive  relationship  approach.  The  maximum  force 
normal  to  the  wall  (Fw)max  computed  using  the  best-fit  straight  line  of  the  scaled 
impulsive  relationship  differs  from  the  maximum  force  (Fw)max  values  presented 
in  Table  5.3  by  a  23-percent  underprediction  to  a  36-percent  overestimate  for 
these  eight  impact  tests. 


Table  6.2 

Comparison  Between  Table  5.3  (Fw)max  Values  and  the  Figure  6.2  Best-Fit 
Straight-Line  (Fw)max  Values  of  the  Scaled  Impulsive  Relationship 

Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal 
to  Wall 
fps 

Coefficient 
of  Friction 
Table  5.3 

(F  w)max 

Table  5.3 
kips 

(F  w)max 

Figure  6.2 

Average 

kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.6 

286.63 

390.09 

36 

30 

2.35 

12.19 

0.50 

0.48 

369.15 

402.40 

8 

31 

1.61 

10.6 

0.30 

0.43 

236.20 

240.17 

2 

37 

1.95 

10.29 

0.35 

0.52 

327.27 

282.47 

-14 

38 

1.83 

11.94 

0.38 

0.57 

230.29 

307.03 

33 

39 

1.61 

14.12 

0.39 

0.51 

271.07 

318.51 

17 

41 

2.86 

8.76 

0.44 

0.51 

419.37 

353.22 

-16 

42 

1.83 

17.48 

0.55 

0.09 

577.44 

445.76 

-23 

6.2  ETL  111  0-2-338  Engineering  Procedure 
Comparison  with  Scaled  Impulsive  Force 
Formulation  Results 

ETL  1 1 10-2-338  provides  an  engineering  procedure  for  the  collision  of  a 
barge  flotilla  with  a  rigid  structure.  According  to  the  ETL  1 1 10-2-338 
engineering  procedure,  the  magnitude  of  the  impact  forces  generated  by  a 
particular  collision  event  is  dependent  on  the  mass  (including  hydrodynamic 
added  mass  of  the  barge  flotilla),  the  approach  velocity,  the  approach  angle,  the 
barge  flotilla  moment  of  inertia,  damage  sustained  by  the  barge  structure,  and 
friction  between  the  barge  and  the  wall.  A  major  distinction  between  this 
procedure  and  the  traditional  Navy  method  for  determining  berthing  forces  is  the 
estimation  of  collision  energy  dissipated  in  deformation  of  the  barge  structure 
and  transferred  to  the  rotation  of  the  barge  flotilla.  The  analytical  method  uses  the 


78 


Chapter  6  Scaled  Impulsive  Force  Formulation  and  Empirical  Correlation 


structural  interaction  mechanism  of  Minorsky,  which  provides  an  empirical 
relationship  between  the  (nonrecoverable)  hull  deformation  and  the  energy 
absorbed  in  a  collision.  The  relationship  between  kinetic  energy  lost  in  a  collision 
and  the  volume  of  in-plane  (barge)  material  damaged  is  used  to  determine  impact 
force  as  a  relationship  to  instantaneous  contact  area  of  damaged  structure. 
Additional  details  regarding  this  approach  are  discussed  in  Section  5.4  of 
Chapter  5. 

The  ETL  1 110-2-338  engineering  procedure  was  used  to  compute  values  of 
maximum  impact  force  normal  to  the  wall,  (Fw)max,  for  the  eight  full-scale, 
controlled  barge  impact  experiments  listed  in  Table  5.3.  The  coefficients  of 
friction  used  in  these  computations  were  the  same  as  those  derived  from  field-test 
data  summarized  in  Table  5.3.  The  results  of  these  ETL-based  computations  are 
given  in  Table  6.3  along  with  the  results  from  the  Figure  6.2  scaled  impulsive 
relationship  best-fit  straight  line.  For  example,  the  best-fit  scaled  impulsive 
relationship  for  barge  impact  Experiment  29  produces  an  (/vW*  value  of 
390.09  kips,  and  the  ETL  1 110-2-338  engineering  procedure  produces  an  (Fn)max 
value  of  409.5  kips.  This  corresponds  to  a  5-percent  overprediction  by  the  ETL 
1 1 10-2-338  approach.  The  maximum  force  normal  to  the  wall  (Fw)ntax  computed 
using  the  ETL  1 110-2-338  engineering  procedure  differs  from  the  best-fit  straight 
line  of  the  scaled  impulsive  Figure  6.2  relationship  by  a  1 3-percent 
underprediction  to  a  20-percent  overestimate  for  these  eight  impact  tests. 


Table  6.3 

Comparison  Between  the  Figure  6.2  Best-Fit  Straight-Line  ( Fw)max  Values  of  the 
Scaled  Impulsive  Relationship  and  ( Fw)max  by  the  ETL  1110-2-338  Engineering 
Procedure 


Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal 
to  Wall 
fps 

Coefficient 
of  Friction 
Table  5.3 

(F  w)max 

Figure  6.2 

Average 

kips 

(Fw)max 

ETL  1110- 

2-338* 

kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.6 

390.09 

409.5 

5 

30 

2.35 

12.19 

0.50 

0.48 

402.40 

420.9 

5 

31 

1.61 

10.6 

0.30 

0.43 

240.17 

263.8 

10 

37 

1.95 

10.29 

0.35 

0.52 

282.47 

317.4 

12 

38 

1.83 

11.94 

0.38 

0.57 

307.03 

327.6 

7 

39 

1.61 

14.12 

0.39 

0.51 

318.51 

317.2 

0 

41 

2.86 

8.76 

0.44 

0.51 

353.22 

423,8 

20 

42 

1.83 

17.48 

0.55 

0.09 

445.76 

386.5 

-13 

*  Note:  {Fw)max computed  following  the  engineering  procedure  outlined  in  ETL  1110-2-338  that  uses  Mx=  1,05'Min  the 
longitudinal  direction,  My=  1.4  Min  the  transverse  direction,  and  k  =  1.41,  where  Mis  the  mass  and  /  is  the  moment  of 
inertia  of  the  flotilla. 


Table  6.4  compares  the  values  of  maximum  impact  force  normal  to  the  wall, 
(F\v)max,  for  the  eight  full-scale,  controlled  barge  impact  experiments  computed 
using  the  ETL  1 110-2-338  engineering  procedure  with  the  range  in  values  from 
the  Figure  6.2  scaled  impulsive  relationship.  The  results  given  in  Table  6.4  show 
that  the  maximum  impact  force  normal  to  the  wall,  {F^^,  for  the  eight  full-scale 
experiments  computed  using  the  ETL  1 1 10-2-338  approach  are  within  one 
standard  error  of  the  Figure  6.2  best-fit  straight  line  of  the  scaled  impulsive 
relationship. 
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Table  6.4 

Comparison  Between  the  Figure  6.2  Scaled  Impulsive  Relationship  for  ( Fw)max 
and  {Fw)max  by  the  ETL 1110-2-338  Engineering  Procedure 

Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal 
to  Wall 
fps 

(F w)msx 

Figure  6.2 

Ave-SE 

kips 

(F w)max 

Figure  6.2 

Average 

kips 

(F wjmBx 

Figure  6.2 

Ave+SE 

kips 

(F„)m„  ETL 
1110-2-338* 
kips 

29 

2.2 

12.63 

0.48 

304.77 

390.09 

475.43 

409.5 

30 

2.35 

12.19 

0.50 

317.07 

402.40 

487.73 

420.9 

31 

1.61 

10.6 

0.30 

154.84 

240.17 

325.50 

263.8 

37 

1.95 

10.29 

0.35 

197.15 

282.47 

367.81 

317.4 

38 

1.83 

11.94 

0.38 

221.70 

307.03 

392.36 

327.6 

39 

1.61 

14.12 

0.39 

233.18 

318.51 

403.84 

317.2 

41 

2.86 

8.76 

0.44 

267.89 

353.22 

438.55 

423.8 

42 

1.83 

17.48 

0.55 

360.44 

445.76 

531.10 

386.5 

*  Note:  (Fw)max computed  following  the  engineering  procedure  outlined  in  ETL  1110-2-338  that  uses  Mx  =  1.05  /Win  the 
longitudinal  direction,  My=  1.4 /Win  the  transverse  direction,  and  l0  =  1.41,  where  M is  the  mass  and  /  is  the  moment  of 
inertia  of  the  flotilla. 

6.3  Empirical  Correlation  Between  (Fw)max  and 
Linear  Momentum  Normal  to  the  Wall 

In  this  section,  an  empirical  correlation  between  the  maximum  force  normal 
to  the  wall  and  the  linear  momentum  normal  to  the  wall  is  presented.  The 
maximum  force,  Fw,  used  in  this  correlation  was  obtained  from  the  models 
presented  in  Chapter  5.  Using  values  for  the  maximum  normal  force,  Fw,  and  the 
linear  momentum  normal  to  the  wall  given  in  Table  6.5  (values  from  Chapter  5), 
a  best-fit  straight  line  was  calculated  for  the  eight  experiments.  This  approach 
relates  the  Fw  obtained  from  the  energy  method  directly  to  the  linear  momentum. 
These  results  can  be  compared  to  the  results  from  the  scaled  impulsive 
relationship.  It  is  important  to  mention  that  only  one  data  point  of  the  whole  Fw 
time-history  for  each  of  the  eight  experiments  was  used  to  develop  this  empirical 
correlation.  However,  in  the  scaled  method,  the  whole  Fw  time-history  was  used 
because  the  total  impulse  (A/  +  A2 )  for  each  experiment  was  used  to  determine 
Equation  6.11. 

The  least  squares  regression  procedure  was  used  to  develop  the  best-fit 
straight  line  through  the  eight  data  points  (for  the  eight  impact  experiments)  for 
the  empirical  correlation.  The  line  was  assumed  to  start  at  the  origin  (i.e.,  no 
intercept  term  was  used  for  the  linear  equation).  The  resulting  best-fit  straight 
line  and  the  average  minus  one  standard  error  and  average  plus  one  standard 
error  lines  are  shown  in  Figure  6.3. 

The  statistical  evaluation  was  done  using  the  computer  program  presented  in 
Appendix  B  and  verified  using  the  standard  Excel  statistical  package.  The 
concepts  that  apply  to  the  statistical  theory  are  presented  in  Appendix  D. 
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Table  6.5 

( Fw)max  Empirical  Correlation 

Exp# 

Approach 
Angle,  deg 

Vel 

fps 

Mass 

k-sec!/ft 

mvsin  (8) 
k-sec 

29 

12,63 

2.20 

1865.59 

286.63 

897.42 

30 

12.19 

1865.59 

369.15 

925.73 

31 

10.60 

1865.59 

236.2 

37 

10.29 

1865.59 

327.27 

649.84 

38 

11.94 

1.83 

1865.59 

230.29 

706.32 

39 

14.12 

1.61 

1865.59 

271.07 

732.74 

41 

8.76 

1 

1865.59 

419.37 

812.59 

42 

17.48 

1865.59 

577.44 

1025.48 

The  resulting  best-fit  equation  for  this  set  of  eight  data  values 
is  (Fw  J  =  0.435  mv  sin  8  ■  That  is,  a  coefficient  times  the  linear  momentum 

normal  to  the  wall  determines  the  maximum  force  normal  to  the  wall.  It  is 
interesting  to  note  that  the  Equation  6.1 1  scaled  impulsive  relationship  has  the 
same  form  as  this  new  empirical  equation.  Using  the  average  of  the  y,  (3,  Ath  and 
FRAC  coefficients  presented  in  Table  6.1a,  Equation  6.11  transforms  to 


(Fw) 


max 


yFRAC 


mv  sin  8  = 


0.06279 

0.72224 


*  0.8 
*  0.L6 


mv  sin  8  =  0.43469 


mv  sin  9 


Note  that  the  coefficients  are  nearly  the  same  value  for  the  two  approaches. 
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We  can  observe  that  the  greater  the  magnitude  for  the  linear  momentum,  the 
larger  will  be  the  maximum  value  for  the  impact  force  normal  to  the  wall,  as 
shown  in  Figure  6.3.  Figure  6.3  is  a  correlation  based  on  low-velocity,  shallow- 
impact  (up  to  21.1  deg)  experiments  that,  by  definition,  do  not  account  for  factors 
that  manifest  themselves  at  higher  velocities.  Additionally,  no  damage  occurred 
to  the  flotilla  of  barges,  and  no  lashings  broke  during  these  eight  impact 
experiments.  Therefore,  at  this  time  and  pending  the  results  of  additional 
research,  Figure  6.3  is  deemed  applicable  to  a  barge  flotilla  that  has  a  velocity 
normal  to  the  wall  up  to  and  not  exceeding  0.57  fps  (0.39  mph)  with  no  damage 
occurring  during  impact  events,  for  impact  angles  up  to  21 .1  deg,  for  a  barge 
flotilla  with  a  linear  momentum  normal  to  the  wall  between  649.84  and 
1,025.48  k-sec.  The  equation  to  determine  the  maximum  force  normal  to  the  wall 
using  the  empirical  correlation  is 

(fwL  =0.435(mvsin6) 

where  the  units  for  the  mass,  velocity  and  maximum  Fw  are  k-sec2/ft,  ft/sec,  and 
kips,  respectively. 

The  maximum  normal  force  (Fh,)^  by  the  empirical  correlation  is  equal  to 
the  reaction  force  provided  by  the  lock  wall  on  the  barge  flotilla  during  the 
impact.  Note  that,  in  Table  6.5,  the  mass  used  to  develop  the  Figure  6.3 
correlation  of  linear  momentum  normal  to  the  wall  with  values  of  ( Fw)max  is  the 
mass  of  the  barge  flotilla  (given  in  Appendix  A)  and  does  not  include  the 
computation  of  any  hydrodynamic  added  masses.  However,  the  values  of  ( Fw)max 
in  Table  6.5  (which  are  derived  from  the  field-test  data  and  originally  reported  in 
Table  5.3)  that  are  used  in  the  Figure  6.3  correlation  reflect  the  effects  of  the 
Figure  2.3  external  forces  of  Ftom  Fdras,  and  Fhetper  acting  on  the  barge  flotilla,  as 
well  as  the  effects  of  the  hydrodynamic  masses  and  the  inertia  of  the  barge 
flotilla  mass  decelerating  during  impact. 

A  summary  of  the  values  of  the  Table  5.3  maximum  impact  force  normal  to 
the  wall,  ( Fw)max ,  for  the  energy  method  or  the  equilibrium  with  a  fixed 
coefficient  of  friction  values  for  the  eight  full-scale,  controlled  barge  impact 
experiments  and  the  values  from  the  Figure  6.3  best-fit  straight-line  (Fw)max 
values  is  given  in  Table  6.6.  For  example,  the  energy  method  or  the  equilibrium 
with  a  fixed  coefficient  of  friction  value  (=  0.6)  for  the  full-scale,  controlled 
barge  impact  experiment  (No.  29)  produces  an  ( Fw)max  value  of  286.63  kips,  and 
the  best-fit  straight  line  of  the  empirical  correlation  produces  an  (FK)max  value  of 
390.38  kips.  This  corresponds  to  a  36-percent  overprediction  by  the  best-fit 
straight  line  of  the  empirical  correlation  approach.  The  maximum  force  normal  to 
the  wall  ( Fw)max  computed  using  the  best-fit  straight  line  of  the  empirical 
correlation  differs  from  the  maximum  force  (Fw)max  values  presented  in  Table  5.3 
by  a  23-percent  underprediction  to  a  36-percent  overestimate  for  these  eight 
impact  tests. 
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Table  6.6 

Comparison  Between  Table  5.3  (F„)max  Values  and  the  Figure  6.3  Best-Fit 
Straight-Line  (Fw)max  Values  of  the  Empirical  Correlation 


Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal 
to  Wall 
fps 

Coefficient 
of  Friction 
Table  5.3 

(  F w)max 

Table  5.3 
kips 

(F w)max 

Figure  6.3 

Average 

kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.6 

286,63 

390.38 

36 

30 

2,35 

12.19 

0.50 

0.48 

369.15 

402,69 

9 

31 

1.61 

10.6 

0.30 

0.43 

236.20 

240,35 

2 

37 

1.95 

10.29 

0.35 

0.52 

327.27 

282,68 

-14 

38 

1.83 

11.94 

0.38 

0.57 

230.29 

307.25 

33 

39 

1.61 

14.12 

0.39 

0.51 

271.07 

318.74 

18 

41 

2.86 

8.76 

0.44 

0.51 

419.37 

353,48 

-16 

42 

1.83 

17.48 

0.55 

0.09 

577.44 

446.08 

-23 

6.4  ETL  1110-2-338  Engineering  Procedure 
Comparison  With  Empirical  Correlation 
Results 

The  ETL  1 1 1 0-2-338  engineering  procedure  was  used  to  compute  values  of 
maximum  impact  force  normal  to  the  wall,  (Fw)max,  for  the  eight  full-scale, 
controlled  barge  impact  experiments  listed  in  Table  5.3.  The  coefficients  of 
friction  used  in  these  computations  were  the  same  as  those  derived  from  field-test 
data  summarized  in  Table  5.3.  The  results  of  these  ETL-based  computations  are 
gi  ven  in  Table  6.7  along  with  the  results  from  the  Figure  6.3  empirical 
correlation  best-fit  straight  line.  For  example,  the  best-fit  empirical  correlation 
for  barge  impact  Experiment  29  produces  an  {Fw)ntax  value  of  390.38  kips,  and  the 
ETL  1 110-2-338  engineering  procedure  produces  an  (F^^  value  of  409.5  kips. 
This  corresponds  to  a  5-percent  overprediction  by  the  ETL  1110-2-338  approach. 
The  maximum  force  normal  to  the  wall  computed  using  the  ETL  1 1 10-2- 

338  engineering  procedure  differs  from  the  best-fit  straight  line  of  the  empirical 
correlation  Figure  6.3  relationship  by  a  13-percent  underprediction  to  a  20- 
percent  overestimate  for  these  eight  impact  tests. 

Table  6.8  compares  the  values  of  maximum  impact  force  normal  to  the  wall, 

C Fw)max ,  for  the  eight  full-scale,  controlled  barge  impact  experiments  computed 
using  the  ETL  1 1 10-2-338  engineering  procedure  with  the  range  in  values  from 
the  Figure  6.3  empirical  correlation.  The  results  given  in  Table  6.8  show  that  the 
maximum  impact  force  normal  to  the  wall,  (Fw)mcLK,  for  the  eight  full-scale, 
computed  using  the  ETL  1 1 10-2-338  approach  are  within  one  standard  error  of 
the  Figure  6.3  best-fit  straight  line  of  the  empirical  correlation. 
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Table  6.7 

Comparison  Between  the  Figure  6.3  Best-Fit  Straight-Line  ( Fw)max  Values  of  the 
Empirical  Correlation  and  {Fw)max  by  the  ETL  1110-2-338  Engineering  Procedure 

Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal  to 
the  Wall 
fps 

Coefficient 
of  Friction 
Table  5.3 

( F „)  — 

Figure  6.3 

Average 

Kips 

(F„)mar  ETL 
1110-2-338* 
kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.60 

390.38 

409.5 

5 

30 

2.35 

12.19 

0.50 

0.48 

402.69 

420.9 

5 

31 

1.61 

10.6 

0.30 

0.43 

240.35 

263.8 

10 

37 

1.95 

10.29 

0.35 

0.52 

282.68 

317.4 

12 

38 

1.83 

11.94 

0.38 

0.57 

307.25 

327.6 

7 

39 

1.61 

14.12 

0.39 

0.51 

318.74 

317.2 

-1 

41 

2.86 

8.76 

0.44 

0.51 

353.48 

423.8 

20 

42 

1.83 

17.48 

0.55 

0.09 

446.08 

386.5 

-13 

*  Note:  (Fw)mBX computed  following  the  engineering  procedure  outlined  in  ETL  1110-2-338  that  uses  M,=  1.05  Min  the 
longitudinal  direction,  My=  1.4'Min  the  transverse  direction,  and  l0  =  1.41,  where  Mis  the  mass  and  /  Is  the  moment  of 
inertia  of  the  flotilla. 

Table  6.8 

Comparison  Between  the  Figure  6.3  Empirical  Correlation  for  (Fw)max  and  ( Fw)max 
by  the  ETL 1110-2-338  Engineering  Procedure 


Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

deg 

Velocity 
Normal  to 
the  Wall 
fps 

(Fv«)m»x 

Figure  6.3 

Ave-SE 

kips 

(Fw)m»x 

Figure  6.3 

Average 

kips 

(F  w)mtix 

Figure  6.3 
Ave+SE 

Kips 

(Fw)mo,  ETL 
1110-2-338* 
kips 

29 

2.2 

12.63 

0.48 

305.05 

390.38 

475.71 

409.5 

30 

2.35 

12.19 

0.50 

317.36 

402.69 

488.02 

420.9 

31 

1.61 

10.6 

0.30 

155.05 

240.35 

325.67 

263.8 

37 

1.95 

10.29 

0.35 

197.35 

282.68 

368.01 

317.4 

38 

1.83 

11.94 

0.38 

221.92 

307.25 

392.58 

327.6 

39 

1.61 

14.12 

0.39 

233.41 

318.74 

404.07 

317.2 

41 

2.86 

8.76 

0.44 

268.15 

353.48 

438.80 

423.8 

42 

1.83 

17.48 

0.55 

360.76 

446.08 

531.41 

386.5 

*  Note:  (F„)max computed  following  the  engineering  procedure  outlined  in  ETL  1110-2-338  that  uses  Mx  =  1.05  /Win  the 
longitudinal  direction,  My  =  1.4 /Win  the  transverse  direction,  and  /0  =  1.41,  where  Mis  the  mass  and  /  is  the  moment  of 
inertia  of  the  flotilla.  _  ...  -  ••  y  .  .  ■ _ 


6.5  Empirical  Correlation  Between  (Fw)max  and 
Kinetic  Energy  Normal  to  the  Wall 

In  this  section  an  empirical  correlation  between  the  maximum  force  normal 
to  the  wall  and  the  kinetic  energy  normal  to  the  wall  of  the  barge  flotilla  is 
presented. 

1  2 

The  kinetic  energy  is  defined  as  T  =  -mv  ,  where  m  is  the  mass  of  the 

object  and  v  is  the  speed  the  object.  This  is  a  positive  scalar  quantity  since  it  does 
not  depend  on  the  direction  of  the  particle’s  velocity.  Furthermore,  the  kinetic 
energy  has  the  same  units  as  work,  e.g.,  joules  (J)  or  foot-pounds,  force  (ft-lb). 
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The  force  Fw  used  in  this  correlation  was  obtained  from  the  models  presented  in 
Section  5.3  of  Chapter  5.  Using  the  values  of  maximum  normal  force  (Fw)max 
from  eight  of  the  impact  experiments  and  the  kinetic  energy  concept  shown  in 
Table  6.9  (values  from  Chapter  5),  a  best-fit  straight  line  was  calculated.  This 
approach  relates  the  (Fw)nua  values  obtained  from  the  energy  method  (see  Section 
5.3)  interpretation  of  the  bumper  instrumentation  data  directly  with  the  kinetic 
energy.  A  least  squares  regression  procedure  was  used  to  develop  the  best-fit 
straight  line  through  the  eight  data  points  (for  the  eight  impact  experiments).  The 
line  was  assumed  to  start  at  the  origin  (i.e.,  no  intercept  term  was  used  for  the 
linear  equation).  The  resulting  best-fit  straight  line  and  the  average  minus  one 
standard  error  and  average  plus  one  standard  error  lines  are  shown  in  Figure  6.4, 


Table  6.9 

Comparison  Between  Table  5.3  (Fw)max  Values  and  the  Figure  6.4  Best-Fit 
Straight-Line  {Fw)max  Values  of  the  Empirical  Correlation  Using  the  Kinetic 

Energy  Concept 

Experiment 

Number 

Approach 

Velocity 

fps 

Approach 

Angle 

Deg 

Velocity 
Normal 
to  Wall 
fps 

Coefficient 
of  Friction 
Table  5.3 

(F  w)max 

Table  5.3 
kips 

(F w)max 

Figure  6.4 

Average 

kips 

Percent 

Difference 

29 

2.2 

12.63 

0.48 

0.60 

286.63 

406.62 

42 

30 

2.35 

12.19 

0.50 

0.48 

369.15 

441.21 

31 

1.61 

10.6 

0.30 

0.43 

236.20 

158.84 

-33 

37 

1.95 

10.29 

0.35 

0.52 

327.27 

216.19 

-34 

38 

1.83 

11.94 

0.38 

0.57 

230.29 

254.84 

11 

39 

1.61 

14.12 

0.39 

0.51 

271.07 

-1 

41 

2.86 

0.44 

0.51 

419.37 

-19 

42 

1.83 

0.55 

577.44 

E11SI 

-8 

Figure  6.4.  Empirical  correlation  using  the  kinetic  energy  concept 
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The  statistical  evaluation  was  done  using  the  computer  program  presented  in 
Appendix  B  and  verified  using  the  standard  Excel  statistical  package.  The 
concepts  that  apply  to  the  statistical  theory  are  presented  in  Appendix  D.  It  can 
be  observed  that,  the  greater  the  magnitude  for  the  kinetic  energy,  the  larger  will 
be  the  maximum  value  for  the  impact  force  normal  to  the  wall,  as  shown  in 
Figure  6.4.  Figure  6.4  is  a  correlation  based  on  low-velocity,  shallow-impact  (up 
to  21.1  deg)  experiments  that,  by  definition,  do  not  account  for  factors  that 
manifest  themselves  at  higher  velocities.  Additionally,  no  damage  occurred  to  the 
flotilla  of  barges,  and  no  lashings  broke  during  these  eight  impact  experiments. 

At  this  time  and  pending  the  results  of  additional  research,  Figure  6.4  is  deemed 
applicable  to  a  barge  flotilla  that  has  a  velocity  normal  to  the  wall  up  to  and  not 
exceeding  0.57  fps  (0.39  mph)  with  no  damage  occurring  during  impact  events, 
for  impact  angles  up  to  21.1  degrees,  and  for  a  barge  flotilla  with  kinetic  energy 
normal  to  the  wall  between  83.95  and  282.17  k-ft. 

Note  that  v  is  the  velocity  normal  to  the  wall  in  the  equations  given  in 
Figure  6.4.  We  can  observe  that,  the  greater  the  magnitude  for  the  kinetic  energy 
normal  to  the  wall,  the  larger  will  be  the  maximum  value  for  the  impact  force 
normal  to  the  wall,  as  shown  in  Figure  6.4.  In  these  expressions,  the  units  for  the 
mass,  velocity  and  maximum  Fw  are  k-sec2/ft,  ft/sec,  and  kips,  respectively. 

The  maximum  normal  force  (Fw)max  by  the  empirical  correlation  is  equal  to 
the  reaction  force  provided  by  the  lock  wall  on  the  barge  flotilla  during  the 
impact.  Note  that,  in  Table  6.9,  the  mass  used  to  develop  the  Figure  6.4 
correlation  of  kinetic  energy  normal  to  the  wall  with  values  of  ( F„.)max  is  the  mass 
of  the  barge  flotilla  (given  in  Appendix  A)  and  does  not  include  the  computation 
of  any  hydrodynamic  added  masses.  However,  in  Table  6.9,  the  values  of  (Fw)mnr 
(which  are  derived  from  the  field-test  data  and  originally  reported  in  Table  5.3) 
that  are  used  in  the  Figure  6.4  correlation  reflect  the  effects  of  the  Figure  2.3 
external  forces  of  Ftow,  Fdrag,  and  Fhelper  acting  on  the  barge  flotilla  acting,  as  well 
as  the  effects  of  the  hydrodynamic  masses  and  the  inertia  of  the  barge  flotilla 
mass  decelerating  during  impact. 

A  summary  of  the  values  of  the  Table  5.3  maximum  impact  force  normal  to 
the  wall,  ( Fw)max ,  for  the  energy  method  or  the  equilibrium  with  a  fixed 
coefficient  of  friction  values  for  the  eight  full-scale,  controlled  barge  impact 
experiments  and  of  the  Figure  6.4  best-fit  straight-line  ( Fw)max  values  is  given  in 
Table  6.9.  For  example,  the  energy  method  or  the  equilibrium  with  a  fixed 
coefficient  of  friction  value  (=  0.6)  for  the  full-scale,  controlled  barge  impact 
Experiment  29  produces  an  (Fw)max  value  of  286.63  kips,  and  the  best-fit  straight 
line  of  the  empirical  correlation  using  the  kinetic  energy  concept  produces  an 
(Fw)max  value  of  406.62  kips.  This  corresponds  to  a  42-percent  overprediction  by 
the  best-fit  straight  line  of  the  empirical  correlation  approach.  The  maximum 
force  normal  to  the  wall  (Fh,)^  computed  using  the  best-fit  straight  line  of  the 
empirical  correlation  differs  from  the  maximum  force  (Fw)max  values  presented  in 
Table  5.3  by  a  34-percent  underprediction  to  a  42-percent  overestimate  for  these 
eight  impact  tests. 
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6.6  Summary,  Conclusions  and 
Recommendation 

This  chapter  summarizes  the  research  and  development  of  three  engineering 
procedures  to  predict  the  value  for  the  maximum  impact  force,  acting 

normal  to  the  lock  wall  during  impact;  the  scaled  method  (Section  6.1);  the 
empirical  correlation  using  linear  momentum  normal  to  the  wall  (Section  6.3); 
and  the  empirical  correlation  using  kinetic  energy  normal  to  the  wall  (Section 
6.5).  After  a  careful  evaluation  of  all  three  formulations  to  calculate  (Fw)niax,  the 
authors  recommend  that  the  empirical  correlation  using  linear  momentum  normal 
to  the  wall  (given  in  Figure  6.3)  be  used.  Recall  that  Figure  6.3  is  a  correlation 
based  on  low -velocity,  shallow-impact  (up  to  21 .1  deg)  experiments  that,  by 
definition,  do  not  account  for  factors  that  manifest  themselves  at  higher 
velocities.  Additionally,  no  damage  occurred  to  the  flotilla  of  barges  during  these 
impact  experiments.  Therefore,  at  this  time  and  pending  the  results  of  additional 
research,  Figure  6.3  is  deemed  applicable  to  a  barge  flotilla  that  has  a  velocity 
normal  to  the  wall  up  to  and  not  exceeding  0.57  fps  (0.39  mph)  with  no  damage 
occurring  during  impact  events,  for  impact  angles  up  to  21.1  deg,  and  for  a  barge 
flotilla  with  a  linear  momentum  normal  to  the  wall  between  649.84  and 
1,025.48  k-sec. 

The  ETL  1 110-2-338  engineering  procedure  was  used  to  compute  values  of 
maximum  impact  force  normal  to  the  wall,  (Fw)max,  for  the  eight  full-scale, 
controlled  barge  impact  experiments.  A  key  aspect  of  this  engineering 
formulation  is  computation  of  collision  energy  dissipated  in  nonrecoverable, 
plastic  hull  deformation  of  (i.e.,  damage  to)  the  comer  of  the  barge  where  impact 
with  the  wall  occurs.  Note  that  no  damage  was  observed  to  the  barge  comer 
during  any  of  these  low-velocity,  controlled-impact  experiments  at  Robert  C. 
Byrd  Lock.  After  a  careful  evaluation  of  the  ETL  1 1 10-2-338  formulation  and 
results  of  the  computations  given  in  this  chapter  and  in  Chapter  5,  the  authors 
recommend  that  this  engineering  procedure  not  be  used  when  damage  to  the 
barge  will  not  occur  during  impact. 
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7  Possible  Layout  for  Future 
Experiments 


The  experiments  that  form  the  basis  of  this  research  contain  many  unknowns, 
producing  a  sensitive  indeterminate  system  that  complicates  the  interpretation  of 
the  measured  results.  The  force  exerted  by  the  towboat  and  helper  boat,  the  drag 
force,  and  the  orientation  of  the  measured  forces  at  the  load  cells  in  the  bumper 
are  some  examples  of  the  variables  affecting  the  problem  formulation.  Sugges¬ 
tions  are  made  in  the  following  list  to  avoid  some  of  the  problems  encountered 
with  using  the  results  from  the  December  1998  experiments  in  the  creation  of  the 
numerical  models  in  this  barge  impact  research  project. 

a.  Measure  four  reactions  and  the  exact  orientation  to  avoid  the 
approximation  of  the  orientation  of  the  forces  F10  and  Fu.  This  involves 
the  use  of  a  pair  of  two-axis  load  cell  pins. 

b.  Measure  layout  geometry  of  bumper  on  barge. 

c.  Take  pictures  directly  vertical  over  the  impact  zone. 

d.  Do  not  use  a  helper  boat. 

e.  Align  a  pair  of  GPS  monitoring  points  parallel  to  the  side  of  the  lead 
impact  barge. 

/.  Analyze  the  first  experiments  to  adjust  the  system  in  the  field  prior  to 
finishing  the  test. 

Table  7.1  presents  some  suggested  values  for  future  experiments.  These 
values  consider  a  range  for  the  number  of  barges  in  the  flotilla,  the  approach 
velocity,  and  the  approach  angle  to  produce  a  linear  momentum  below  the  values 
reported  in  the  experiments  with  a  15-barge  flotilla.  The  mass  per  barge  was 
calculated  based  on  the  mass  of  the  15-barge  flotilla  divided  by  15  to  obtain  a 
unit  mass  value. 
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Table  7.1 

Testing  Range  Values  for  Future  Experiments 

Number  of 

Barges 

Mass  per  Barge 
k-s2/ft 

Approach 
Velocity,  ft/s 

Approach 

Angle 

deg 

Linear 

Momentum,  k-s 

1 

130 

6 

30 

390 

2 

130 

5 

30 

650 

4 

130 

4 

20 

711 

8 

130 

3 

15 

808 
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8  Conclusions  and 
Recommendations 


In  this  research  project,  a  complex  dynamic  problem  was  modeled  using  a 
single  mass  model  to  determine  the  force  applied  to  the  lock  wall  by  a  barge 
flotilla.  The  input  data  for  the  research  were  obtained  from  the  experiments 
conducted  in  1998  at  the  Robert  C.  Byrd  Lock  in  West  Virginia.  The  input  data 
were  the  approach  angle,  the  approach  velocity,  and  the  force  time-histories 
recorded  in  the  experiment.  The  approach  velocity  and  angle  were  obtained  using 
a  global  positioning  system. 

The  load  time-histories  were  obtained  using  an  instrumented  arc  beam.  This 
beam  was  made  of  A-36  structural  steel  and  placed  at  the  comer  of  the  barge 
using  a  pin-pin  connection.  The  loads  were  recorded  in  a  direction  other  than  the 
radial  direction  of  the  arc  bumper.  To  determine  the  maximum  force  applied  to 
the  wall,  the  impulse  and  linear  momentum  principle  was  applied.  However,  the 
system  of  equations  that  this  principle  uses  was  not  enough  to  solve  the  complete 
time-history  of  impact  from  time  t\  to  time  t2.  There  were  more  unknown  forces 
in  the  system  than  equations.  Two  sets  of  methodologies  were  used  to  interpret 
the  barge  impact  data. 

One  of  the  methodologies  used  to  determine  the  normal  load  applied  to  the 
lock  wall  by  the  barge  flotilla  involved  the  use  of  the  equations  of  equilibrium 
combined  with  an  assumed  value  for  the  coefficient  of  friction.  Another 
methodology  used  to  determine  the  normal  load  applied  to  the  lock  wall  by  the 
barge  flotilla  involved  the  relation  between  the  linear  momentum  and  the  impulse 
normal  to  the  wall  produced  by  the  force  computed  using  the  model  referred  to  as 
the  assumed  coefficient  of  friction  method,  or  the  energy  method.  An  empirical 
correlation  was  developed  relating  the  maximum  force  on  the  bumper  (normal  to 
the  wall)  to  the  linear  momentum  m  •  v  •  sin  0  at  time  of  initial  impact. 

Three  possible  bumper  configurations  were  analyzed  based  on  the  location  of 
the  bumper  in  relation  with  the  longitudinal  axis  of  the  barges  and  the  load  cell 
orientations.  The  first  configuration  was  based  on  the  configuration  reported  at 
the  time  of  the  experiments.  The  bumper  was  placed  54  deg  from  the  longitudinal 
axis  of  the  barges,  and  the  load  cells  were  oriented  in  the  radial  direction.  This 
configuration  produced  some  unreasonable  values  of  coefficient  of  friction  and 
can  produce  tension  in  the  arc,  which  is  physically  impossible.  After  a  detailed 
study  of  drawings  and  photographs,  a  new  configuration  was  estimated.  In  this 
second  case,  the  arc  was  considered  to  be  placed  in  the  same  position  as  the  first 
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configuration,  but  the  load  cells  were  oriented  in  a  direction  outside  the  radial 
direction.  An  unreasonable  solution  was  again  obtained  in  terms  of  the 
coefficient  of  friction  (greater  than  1 ,  or  negative  values)  and  tension  in  the  arc. 
Then,  a  final  (third)  configuration  was  obtained  based  on  the  values  of  coefficient 
of  friction  (between  0  and  1)  and  compression  in  the  arc.  With  this  final 
configuration,  the  force  normal  to  the  wall  time-histories  was  calculated  for  eight 
experiments,  and  an  empirical  relationship  using  the  impulse  and  linear 
momentum  principle  was  obtained  (Figure  6.3)  to  calculate  the  maximum  normal 
force  at  the  wall  due  to  the  flotilla  impact. 

Chapter  6  summarized  the  development  of  three  engineering  procedures  to 
predict  the  value  for  the  maximum  impact  force,  (FJ^,  acting  normal  to  the 
lock  wall  during  impact:  the  scaled  method  (Section  6.1);  the  empirical 
correlation  using  linear  momentum  normal  to  the  wall  (Section  6.3);  and  the 
empirical  correlation  using  kinetic  energy  normal  to  the  wall  (Section  6.5).  After 
a  careful  evaluation  of  all  three  formulations  to  calculate  (FJ^,  the  authors 
recommend  that  the  empirical  correlation  using  linear  momentum  normal  to  the 
wall  given  in  Figure  6.3  be  used.  Recall  that  Figure  6.3  is  a  correlation  based  on 
low-velocity,  shallow-impact  (<21.1  deg)  experiments  that,  by  definition,  do  not 
account  for  factors  that  manifest  themselves  at  higher  velocities.  Additionally,  no 
damage  occurred  to  the  flotilla  of  barges  during  these  impact  experiments. 
Therefore,  at  this  time  and  pending  the  results  of  additional  research.  Figure  6.3 
is  deemed  applicable  to  a  barge  flotilla  that  has  a  velocity  normal  to  the  wall  up 
to  and  not  exceeding  0.57  fps  (0.39  mph)  with  no  damage  occurring  during 
impact  events,  for  impact  angles  up  to  21 .1  deg,  and  for  a  barge  flotilla  with  a 
linear  momentum  normal  to  the  wall  between  649.84  and  1,025.48  k-sec. 

The  ETL  1 1 10-2-338  engineering  procedure  was  also  used  to  compute  values 
of  maximum  impact  force  normal  to  the  wall,  (FJ^,  for  the  eight  full-scale, 
controlled  barge  impact  experiments.  A  key  aspect  of  this  engineering 
formulation  is  computation  of  collision  energy  dissipated  in  nonrecoverable, 
plastic  hull  deformation  of  (i.e.,  damage  to)  the  comer  of  the  barge  where  impact 
with  the  wall  occurs.  No  damage  was  observed  to  the  barge  comer  during  any  of 
these  low-velocity,  controlled-impact  experiments  at  Robert  C.  Byrd  Lock.  After 
a  careful  evaluation  of  the  ETL  1 110-2-338  formulation  and  the  results  of  the 
computations  given  in  Chapters  5  and  6,  the  authors  recommend  that  this 
engineering  procedure  not  be  used  when  damage  to  the  barge  will  not  occur 
during  impact. 

The  interpretation  described  in  this  report  did  not  include  an  evaluation  of 
how  the  presence  of  the  bumper  alters  the  barge -to-lock  wall  impacts  or 
consideration  of  any  flexibility  actions  of  the  lock  wall  (e.g.,  rocking  on  its 
foundation)  during  the  impact  event.  The  effect  of  wall/system  flexibility  during 
impact  events  as  well  as  the  response  of  a  multi-degree-of-freedom  system  (to 
account  for  the  lashings  among  the  barges  that  comprise  the  flotilla  during  an 
impact  event)  are  important  issues  that  will  be  the  subject  of  future  research.  This 
series  of  events  will  be  interpreted  using  acceleration  time-histories  that  were 
recorded  during  impact  events  using  Newton’s  second  law  (F  =  ma ). 

Finally,  the  suggestions  made  in  Chapter  7  should  be  considered  if  future 
instrumented  barge  impact  experiments  are  to  be  conducted.  The  use  of  simple 
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geometric  bumper  configurations  and  the  use  of  four  load  cells  is  recommended 
to  address  the  problem  of  the  measured  nonradial  forces  at  the  pin-pin 
connections. 
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Appendix  A 

Summary  of  Experimental 
Results 


The  data  and  analysis  presented  in  this  appendix  were  generated  from  field 
data  acquired  during  the  full-scale  barge  impact  experiments  conducted  at 
Robert  C.  Byrd  Lock  and  Dam  in  December  1998.  The  data  presented  here 
represent  only  the  portion  relevant  to  this  report,  mainly  the  load  bumper  and 
global  positioning  system  (GPS)  data. 

Figure  A.l  shows  the  tow  configuration  and  orientation  for  the  barge  impact 
experiments.  The  flotilla  consisted  of  15  jumbo  open-hopper  barges  (35 
by  195  ft),  ballasted  with  anthracite  coal  to  a  draft  of  9  ft,  a  typical  configuration 
for  commercial  traffic  on  the  Ohio  River  main  stem.  The  front  and  rear  rows 
were  single-raked  barges  (approximately  23-deg  rake  at  the  front,  flat  at  the 
back).  The  middle  nine  barges  were  double  raked  (rakes  at  both  ends).  The 
drawing  below  shows  the  two  coordinate  systems  used  in  this  analysis:  the  global 
axis  represents  the  lock  wall,  and  the  local  axis  is  on  the  tow.  The  impact  angle 
(0)  is  the  angle  of  rotation  between  these  axes.  The  impact  angle,  0,  referenced 
throughout  this  report  is  the  angle  between  the  port  side  of  the  right-front  barge 
and  the  lock  wall,  as  measured  with  the  GPS  systems  (see  Table  A.1). 

Force  time-histories  for  the  load-bumper  gages  F10  (front)  and  FI  1  (rear)  for 
Experiments  22-31  and  37-44  are  presented  in  this  appendix.  Full-duration  plots 
for  all  experiments  are  included,  as  well  as  expanded  plots  of  most  experiments 
showing  the  first  load  cycle.  Table  A.l  summarizes  the  significant  aspects  of 
these  plots.  In  Table  A.2,  rise  time  (defined  here  as  AT  between  contact  and  peak 
load)  is  listed  for  each  measurement  along  with  the  time  of  contact. 

Figure  A.2  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  22  during  the  9  sec  of  impact.  This  experiment  was 
conducted  at  a  shallow  approach  angle  of  10  deg  and  prior  to  removal  of  the  flat- 
pack  skid  plate  along  the  front  face  of  the  bumper.  For  this  shallow  approach 
angle,  the  rear  load  cell  absorbed  nearly  all  the  load  during  the  impact,  compared 
with  the  near-zero  loads  recorded  by  the  front  load  cell.  Note  that  the  initial 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  xiii  of  the  main  text. 
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Figure  A.1 .  Experiment  configuration  diagram 
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r  AT  =  Ti  at  time  of  maximum  force  (load  cell  10  for  experiment  26,  40,  43  and  44,  load  cell  11  otherwise)  minus  time  of  minimum 


For  the  statistical  parameters  experiment  40  not  used. 


impact  force  pulse  contains  the  maximum  magnitude  force  that  is  recorded 
during  the  experiment  for  both  load  cells.  This  observation  is  valid  for  all 
experiments. 

Figure  A.3  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  27.  This  experiment  was  conducted  at  a  shallow 
approach  angle  of  8  deg  at  a  location  farther  up  the  wall  from  the  primary  target 
site,  for  the  purpose  of  “shearing  off’  the  damaged  flat-pack  bracket  hardware. 
This  experiment  should  not  be  used  for  evaluation  purposes.  Only  those 
experiments  designated  with  numbers  above  27  were  used  in  the  data  reduction 
that  is  discussed  in  this  appendix. 
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Figure  A.2.  Experiment  22-Time  variation  of  recorded  bumper  forces  FI  0  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  0.88  fps  and 
a  10-deg  impact  angle 
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Figure  A.3.  Experiment  27-Time  variation  of  recorded  bumper  forces  FI  0  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  2.2  fps  and 
an  8-deg  impact  angle 


Figures  A.4-A.6  show  the  force  time-histories  of  the  front  load  cell  F10  and 
rear  load  cell  FI  1  for  Experiments  23-25,  respectively.  These  experiments  were 
conducted  at  shallow  approach  angles  ranging  from  15.5  to  18  deg.  Note  that  the 
maximum  forces  were  recorded  by  the  rear  load  cell  (FI  1)  during  these 
experiments. 

Figure  A.7  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  26.  This  experiment  was  conducted  at  a  steep 
approach  angle  of  23.75  deg.  Note  that  the  magnitude  of  the  forces  recorded  by 
the  front  load  cell  (F10)  during  impact  exceeds  the  forces  recorded  by  the  rear 
load  cell  (FI  1).  The  authors  believe  that  the  front  load  cell  recorded  the  largest 
force,  as  compared  to  that  recorded  by  the  rear  load  cell,  due  to  the  proximity  of 
the  approach  angle  of  the  flotilla  in  comparison  with  the  midpoint  of  the  arc  of 
the  bumper.  Midpoint  of  the  arc  geometry  of  the  bumper  is  equal  to  21.5  deg  (the 
average  of  7  and  36  deg).  The  point  of  impact  is  a  function  of  the  approach  angle 
and  the  mounting  coordinates  (in  terms  of  angles  for  the  bumper 
connection/reaction  points  on  the  barge),  and  the  approach  angle  of  the  flotilla 
and  barge  on  which  the  bumper  is  mounted.  The  approach  angle  of  23.75  deg  is 
in  excess  of  the  midpoint  angle  of  21 .5  deg  along  the  arc  of  the  bumper;  thus,  the 
front  bumper  load  cell  F10  will  record  the  maximum  force. 
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Figure  A.4.  Experiment  23-Time  variation  of  recorded  bumper  forces  FI  0  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  0.81  fps  and 
a  15.5-deg  impact  angle 
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Figure  A.5.  Experiment  24-Time  variation  of  recorded  bumper  forces  FI  0  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  1.1  fps  and 
an  18-deg  impact  angle 
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Figures  A.8-A.13  show  the  force  time-histories  of  the  front  load  cell  F10  and 
rear  load  cell  FI  1  for  Experiments  28-31  and  Experiments  37  and  38, 
respectively.  These  experiments  were  conducted  at  shallow  approach  angles 
ranging  from  12.5  to  15  deg.  Note  that  the  maximum  forces  were  recorded  by  the 
rear  load  cell  (FI  1)  during  these  experiments. 

Figure  A.  14  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  Fll  for  Experiment  39.  This  experiment  was  conducted  at  a  moderate 
approach  angle  of  17.25  deg.  Note  that  the  maximum  forces  were  recorded  by  the 
rear  load  cell  (Fll)  for  this  experiment.  During  the  first  force  pulse,  the  rear  load 
cell  FI  1  records  a  maximum  force  that  is  nearly  twice  the  magnitude  of  force 
recorded  by  the  front  load  cell  F10. 

Figure  A.  15  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  40.  This  experiment  was  conducted  at  a  moderate 
approach  angle  of  20.25  deg.  Note  that  during  the  first  two  impact  force  pulses, 
the  maximum  forces  were  nearly  the  same  for  both  load  cells,  with  the  front  load 
cell  (F10)  recording  slightly  larger  forces.  This  trend  was  reversed  during  the 
subsequent  impact  pulses.  Recall  that  the  midpoint  along  the  arc  of  the  bumper  is 
21.5  deg  and  that  the  approach  angle  in  this  experiment  is  close  to  this  value.  The 
trend  of  the  rear  load  cell  recording  slightly  larger  values  than  the  front  load  cell 
a  few  seconds  after  impact  commences  implies  that  there  was  a  change  in  angle 
of  the  flotilla  (relative  to  the  wall)  during  the  9  sec  of  impact. 

Figure  A.16  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  41.  This  experiment  was  conducted  at  a  shallow 
approach  angle  of  1 1 .5  deg.  Note  that  the  maximum  forces  were  recorded  by  the 
rear  load  cell  (FI  1)  during  this  experiment. 

Figure  A.  17  shows  the  force  time-histories  of  the  front  load  cell  F10  and  rear 
load  cell  FI  1  for  Experiment  42.  This  experiment  was  conducted  at  a  moderate 
approach  angle  of  18.5  deg.  Note  that  the  maximum  forces  were  recorded  by  the 
rear  load  cell  (Fll)  for  this  experiment.  During  the  first  force  pulse,  the  rear  load 
cell  Fll  records  a  maximum  force  but,  because  of  the  moderate  approach  angle, 
the  front  load  cell  F10  records  a  maximum  force  that  is  only  about  25  percent 
less  than  the  maximum  value  recorded  by  the  rear  load  cell. 

Figures  A.18  and  A.19  show  the  force  time-histories  of  the  front  load  cell 
F10  and  rear  load  cell  FI  1  for  Experiments  43  and  44,  respectively.  These 
experiments  were  conducted  at  steep  approach  angles  (25  and  23  deg, 
respectively).  Note  that  the  maximum  forces  were  recorded  by  the  front  load  cell 
(F10)  for  these  two  experiments.  During  the  first  force  pulse,  the  rear  load  cell 
FI  1  records  a  maximum  force  that  is  only  about  10  percent  less  than  the 
maximum  value  recorded  by  the  front  load  cell  in  Experiment  43  and  about 
40  percent  less  for  Experiment  44. 
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Figure  A. 10.  Experiment  30-Time  variation  of  recorded  bumper  forces  F10  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  2.35  fps  and 
a  15-deg  impact  angle 
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Figure  A.1 5.  Experiment  40-Time  variation  of  recorded  bumper  forces  FI  0  and 
F1 1  during  barge  impact  with  an  approach  velocity  of  1 ,91  fps  and 
a  20.25-deg  impact  angle 
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Figure  A.1 9.  Experiment  44-Time  variation  of  recorded  bumper  forces  FI  0 
and  F1 1  during  barge  impact  with  an  approach  velocity  of 
1 .22  fps  and  a  23-deg  impact  angle 
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The  bumper  force  data  reduction  and  data  assessment,  summarized  in 
Table  A.l,  demonstrate  that,  for  Experiments  29-31  and  37-44,  the  maximum 
force  values  recorded  by  the  front  and  rear  load  cells  F10  and  FI  1 ,  respectively, 
occurred  at  approximately  the  same  time. 

This  data  reduction  and  data  assessment  also  shows  that  the  duration  of 
impact  between  the  flotilla  and  the  wall  was  about  9  sec  (Table  A.2).  The  time  to 
maximum  load  was  always  observed  to  be  recorded  during  the  first  impact  force 
pulse  for  these  experiments  and  occurred,  on  average,  0.186  sec  after  initial 
contact  was  made  with  the  wall,  with  a  standard  deviation  of  0.042  sec  and  a 
coefficient  of  variation  of  0.25  (Table  A.2). 

Some  comments  regarding  the  development  of  Table  A.3  should  be  noted. 
The  exact  angle  of  impact  is  critical  in  determining  the  reaction  forces  acting  on 
the  barge  comer.  To  this  end,  the  differential  GPS  data  (DGPS)  documenting  the 
position  of  the  comer  barge  have  been  carefully  analyzed  to  identify  the  exact 
time  that  impact  occurred,  as  well  as  the  precise  orientation  of  the  barge  at 
impact.  This  table  presents  the  angle  for  the  last  5  sec  (points)  prior  to  impact. 
The  column  “Theta  GPS  ”  gives  the  angle  formed  by  the  line  between  GPS 
stations  1  and  2  and  the  lock  wall.  The  correction  value  listed  in  the  “Theta  Corr” 
column  is  that  needed  to  bring  this  line  parallel  to  the  barge’s  port  side.  “Theta 
AVE”  is  then  found  for  the  five  points.  Velocity  is  taken  from  Table  A.l . 

The  data  in  Table  A.4  summarize  the  individual  barge  weights  of  the  15 
barges,  towboat,  and  helper  boat.  The  total  mass  is  1,865.59  k-sec2/ft,  equal  to  the 
total  weight  divided  by  the  gravitational  constant,  g. 

The  precise  tow  position  during  the  44  impact  experiments  was  determined 
using  kinematic  DGPS  instruments.  These  instruments  were  located  on  the  barge 
flotilla,  as  shown  in  Figure  A.20.  The  two  most  significant  units  (GPS-1  and 
GPS-2)  were  deployed  along  the  port  side  of  the  comer  barge,  with  a  third  unit 
(GPS-3)  on  the  stem  area  of  the  center  barge.  These  units  recorded  tow  position 
each  second  to  an  accuracy  of  2  cm.  The  time  and  point  of  impact  were  obtained 
through  the  positioning  plots  of  the  forwardmost  DGPS  station  (Figures  A.21- 
A.30).  This  station  followed  the  general  course  of  the  tow  up  to  the  time  of 
impact.  After  impact,  the  front  of  the  barge  tracked  parallel  to  the  lock  wall. 

Once  the  time  of  impact  was  identified,  the  impact  angle  (the  angle  formed  by 
the  port  side  of  the  comer  barge  with  the  lock  wall)  was  determined  from  the 
GPS  data  and  corrected  for  the  relative  positions  of  GPS  stations  1  and  2  on  the 
barge  (see  Table  A.3).  This  angle  is  critical  to  the  bumper  geometry  and  the 
resulting  force  system.  Velocity  (actually,  speed)  was  calculated  from  the 
displacement  of  the  front  comer  GPS  unit  per  unit  time  (1  sec). 
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Figure  A.20.  Differential  GPS  instrument  locations  for  Experiments  22-44 


Figure  A. 21 .  Route  of  Experiment  29  and  wall  alignment 
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Figure  A.24.  Route  of  Experiment  37  GPS-2  and  wall  alignment  (Note:  GPS-1 
data  were  not  available  for  this  impact) 


Figure  A.25.  Route  of  Experiment  38  GPS-2  and  wall  alignment  (Note:  GPS- 
data  were  not  available  for  this  impact) 
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An  additional  assessment  was  performed  to  determine  the  variation,  if  any,  of 
the  impact  angle  during  the  9  sec  of  impact.  In  Figures  A.31-A.36,  the  route 
recorded  for  GPS-1  and  GPS-2  is  presented.  Three  additional  lines  are  also 
presented.  The  red  line  (on  the  right)  is  the  reference  line  that  joins  the  initial 
point  of  contact  of  both  trajectories.  The  green  line  (on  the  left)is  the  reference 
line  that  joins  the  final  point  of  contact  of  both  trajectories.  It  is  also  indicated  as 
the  line  that  makes  reference  to  the  alignment  of  the  concrete  wall. 

Figures  A.3 1  -A.36  present  (in  the  upper  portion)  the  trajectory  of  both  GPS 
units  for  the  corresponding  experiment.  At  the  lower  portion  of  the  figure,  each 
trajectory  is  considered  separately  to  observe  that  both  angles  (initial  and  final 
angle  with  respect  to  the  concrete  wall)  are  the  same.  The  presented  angles  in 
each  figure  are  very  close  to  the  angle  used  in  the  calculation  in  this  technical 
report.  It  is  important  to  mention  that  the  route  of  GPS-2  was  not  provided  for 
Experiments  37  and  38  and,  for  that  reason,  these  routes  were  not  included  in 
these  plots. 
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Figure  A.32.  GPS-1  and  GPS-2  route  of  Experiment  30  and  wall  alignment 
during  9  sec  of  impact 
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Figure  A.33.  GPS-1  and  GPS-2  route  of  Experiment  31  and  wall  alignment 
during  9  sec  of  impact 
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Figure  A.34.  GPS-1  and  GPS-2  route  of  Experiment  39  and  wall  alignment 
during  9  sec  of  impact 
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Figure  A.35.  GPS-1  and  GPS-2  route  of  Experiment  41  and  wall 
alignment  during  9  sec  of  impact 
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Appendix  B 

FORTRAN  Source  Programs 
and  Maple  Worksheets 


B.1  Program  to  Calculate  the  Forces  in  the  Arc  Beam 
Using  the  Assumed  Coefficient  of  Friction  Model 


Program  COMPONENTS 


*  Program  to  calculate  the  reactions  components  in  the  global  and 

*  local  coordinates  systems .  The  known  values  are  the  radial 

*  components  measured  in  experimental  test.  The  wall  coincides  with 
the  global  "x"  axis  and  the  flotilla  is  inclined  "theta"  degrees 

*  with  the  horizontal  "x"  global  axis.  The  "x"  local  axis  of  the 

*  flotilla  coincides  with  the  longitudinal  axes  of  it,  which  is 

*  inclined  "theta"  degrees  with  the  wall. 


*  Identifiers  used  are: 

*  F10 , Fll  :  known  radial  forces 

NF  :  number  of  forces  in  the  force  time  history 

*  I,J  :  indeces 


integer *4  NLOAD,  I,  J,  NEXP,  NFRIC,  III 
real*8  THETA,  F10,  Fll,  FRIC(IO),  FNEW(IO), 

1  SW,  AX,  AY,  BX,  BY,  RESULA,  RESULB,  ANGA,  ANGB,  SOL(20,20), 

2  SOLM  (20,20)  ,  FG,  FGG,  F10N,  FUN,  Cl ,  C2  ,  C3  ,  C4 ,  C5  ,  C6  ,  a,  b,  c ,  d,  e , 

3  aa,  bb,  cc,  dd,  ee,  f  f ,  gg ,  hh,  ii ,  j  j  ,  kk,  11 ,  mm,  nn,  oo,  pp,  gq,  rr ,  ss ,  tt , 

4  uu,  w,  ww,  xx, yy,  aaa,  bbb,  ccc,  ddd,  eee,  f  f  f 
character*20  NAME1,  NAME 2 ,  NAME 3 ,  NAME 4 

write (*,*)  'Enter  force  input  filename  with  extension' 
read ( * , * )  NAME1 

write {*,*)  'Enter  reactions  output  filename  with  extension' 
read {  *  ,  * )  NAME2 

write(*,*)  'Enter  Fw  output  filename  with  extension' 
read ( * , * )  NAME 3 

write {*,*}  'Enter  Maximum  output  filename  with  extension' 

read ( * , * )  NAME 4 

open (UNIT=1 , FILE=NAME1 ) 

open (UNIT=2 , FILE=NAME2 ) 
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open (UNIT=3 , FILE=NAME3 ) 

open (UNIT=4 , FILE= ' FRICCOEFF . DAT ' ) 

open (UNIT=15 , FILE=NAME4 ) 

read(l,*)  NEXP 

read(l,*)  THETA 

read (4,*)  NFRIC 

read (4, *)  (FRIC(I),  1=1,  NFRIC) 
read (1 , * )  NLOAD 

write {2, 106)  NEXP,  THETA,  NFRIC,  NLOAD 

write (3, 106)  NEXP,  THETA,  NFRIC,  NLOAD 

THETA  =  THETA  *  3  d4159  /  180.0 

Cl  =  (3 .5/180.0) *3 .14159 

C2  =  (32. 5/180. 0)*3. 14159 

C3  =  (57.5/180.0) *3 .14159 

C4  =  (93. 5/180. 0)*3. 14159 

C5  =  (5. 5/180. 0)*3. 14159 

C6  =  (1. 5/180. 0)*3. 14159 

a  =  THETA-C1 

b  =  C2 -THETA 

C  =  C3+THETA 

d  =  C4-THETA 

e  =  C5 

f  =  C6 

aa=cos (a) *sin (e) 
bb=cos (a) *cos (e) 
cc=cos (b) *sin (f ) 
dd=cos (b) *cos (f ) 
ee=sin(a) *cos (e) 
f f =sin (a) *sin (e) 
gg=sin(b) *cos (f ) 
hh=sin (b) *sin(f) 
ii=sin(f ) *sin(b) *cos (c) 
j j=sin(f ) *sin(b) *cos (d) 
kk=sin (f)*cos(b)*sin(d) 
ll=sin(f)*sin(c) *cos(b) 
mm=cos (f)*sin(b) *cos(c) 
nn=cos (f)*sin(b)*cos(d) 
oo=cos (f)*cos(b)*sin(d) 
pp=cos (f ) *cos (b) *sin (c) 
rr=cos (f ) *cos (b) *cos (c) 
ss=cos ( f ) *cos (b) *cos (d) 
tt=cos ( f ) *sin (b) *sin (d) 
uu=cos (f) *sin(b) *sin(c) 
w=sin  ( f )  *cos  (b)  *cos  (c) 
ww=sin ( f ) *cos (b) *cos (d) 
xx=sin(f) *sin(b) *sin(d) 
yy=sin(f) *sin(b) *sin(c) 
aaa=aa-ee 
bbb=bb+f f 
ccc=-cc+gg 
ddd=-dd-hh 

eee=-ii- j j  +kk-ll-rr-ss-tt+uu 

f  f  f = -mm-nn+oo  -pp+ w+ww+xx-yy 

write (3, 107)  (FRIC(I),  1=1,  NFRIC) 

write (2 , 101) 

write (15, 108) 

write (2 , 102 ) 
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n  n 


write(15» 109) 
write (2 ,103) 
write (15, 110) 

FGG  =  -1000.00 
do  I  =  1  ,  NLOAD 

readfl,*)  TIME,F10,F11 
a  =  THETA- Cl 

b  =  C2 -THETA 

c  =  C3+THETA 

d  =  C 4 -THETA 

e  =  C5 

f  =  C6 

aa=cos (a) *sin (e) 
bb=cos (a) *cos (e) 
cc=cos (b) *sin (f ) 
dd=cos (b) *cos (f ) 
ee=sin(a) *cos (e) 
ff=sin(a) *sin(e) 
gg=sin (b) *cos { f ) 
hh=sin (b) *sin (f ) 
ii=sin ( f ) *sin (b) *cos (c) 
j j=sin(f ) *sin(b) *cos (d) 
kk=sin(f ) *cos (b) *sin(d) 
ll=sin ( f ) *sin (c) *cos (b) 
mm=cos ( f ) *sin (b) *cos (c) 
nn=cos ( f ) *sin (b) *cos (d) 
oo=cos (f) *cos (b) *sin (d) 
pp— cos (f ) *cos (b) *sin(c) 
rr=cos (f ) *cos (b) *cos (c) 
ss=cos (f ) *  cos (b) *cos (d) 
tt— cos (f) *sin(b) *sin(d) 
uu=cos (f) *sin(b) *sin(c) 
w=sin  (  f )  *cos  (b)  *cos  (c) 
ww=sin(f ) *cos (b) *cos (d) 
xx=sin { f ) *sin (b) *sin(d) 
yy=sin(f ) *sin (b) *sin(c) 
aaa=aa-ee 
bbb=bb+f f 
ccc=-cc+gg 
ddd=-dd-hh 

eee=-ii- j  j  +kk-ll-rr-ss-tt+uu 
f  f  f = -mm-nn+oo  -pp  + w+ww+xx-yy 
if  (F10  .gt.  Fll )  then 
FG  =  F10 
else 

FG  =  Fll 
endif 

do  J  =  1  ,  NFRIC 

qq  =  cos (d) -FRIC { J) +FRXC ( J) *sin (d) 


Begins  calculation  of  the  matrix  terms 
C 

Ff ON  =  - (qq*aaa**2*Fll+aaa*qq*ccc*F10+aaa*FRIC (J) *eee*F10 
l+qq*bbb**2*Fll-bbb*qq*ddd*F10+bbb*eee*F10) / (aaa*qq*ddd+aaa*FRIC ( J) 
2*f f f +bbb*f f f +bbb*qq*ccc) 

F11N  =  - (f f f *aaa*Fll+f f f *ccc*F10+qq*ccc*aaa*Fll+qq*ccc**2 
1 *F10+ccc*FRIC ( J) *eee*F10-qq*bbb*ddd*Fll-bbb*Fll *FRIC ( J) *f ff + 
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2gq*ddd**2*F10+ddd*F10*FRIC ( J) *f f f-eee*F10*ddd) / (aaa*qq*ddd 
3+aaa*FRIC ( J) *f f f +bbb*f f f +bbb*gq*ccc) 

FW  =  (-eee*F10*aaa*ddd-ccc*bbb*eee*F10+f f f*aaa**2*Fll 
1+f f f *aaa*ccc*F10  +  f f f *bbb**2  *Fll-f f f *bbb*ddd*F10) / 

2 (aaa*qq*ddd+aaa*FRIC ( J) *ff f +bbb*f f f +bbb*qq*ccc) 

FNEW(J)  =  FW 
SW  =  FW  *  FRIC(J) 

SPA  =  Fll *  cos (e) —FI IN*  sin ( e ) 

FNA  =  Fll*sin(e) +FllN*cos (e) 

SPB  =  F10*cos (f ) -FlON*sin(f ) 

FNB  =  F10*sin(f ) +F10N*COS (f ) 

AX  =  FNA*cos (a) -SPA*sin (a) 

BX  =  SPB*sin (b) -FNB*cos (b) 

AY  =  SPA*cos (a) +FNA*sin (a) 

BY  =  SPB*cos (b) +FNB*sin (b) 

RESULA  =  sqrt (AX*AX+AY*AY) 

RESULB  =  sqrt (BX*BX+BY*BY) 

ANGA  =  atan (AY/AX)  *  180  /  3.141 
ANGB  =  atan (BY/BX)  *  180  /  3.141 
SOL ( J, 1 )  =  FRIC(J) 

SOL ( J, 2 )  =  F10 
SOL ( J, 3 )  =  Fll 
SOL ( J, 4 )  =  FNA 
SOL ( J, 5 )  =  FNB 
SOL ( J, 6 )  =  SPA 
SOL ( J, 7 )  =  SPB 
SOL ( J, 8 )  =  FW 
SOL ( J, 9 )  =  SW 
SOL ( J, 10)  =  AX 
SOL (J, 11)  =  BX 
SOL (J, 12)  =  AY 
SOL (J, 13)  =  BY 
SOL (J, 14)  =  RESULA 
SOL (J, 15)  =  RESULB 
SOL (J, 16)  =  ANGA 
SOL (J, 17)  =  ANGB 
if  (FG  .gt.  FGG)  then 
SOLM ( J , 1 )  =  FRIC(J) 

SOLM ( J, 2 )  =  F10 
SOLM { J , 3 )  =  Fll 
SOLM ( J , 4 )  =  FNA 
SOLM (J, 5)  =  FNB 
SOLM (J, 6)  =  SPA 
SOLM (J, 7)  =  SPB 

SOLM (J, 8)  =  FW 
SOLM(J, 9)  =  SW 

SOLM{ J, 10)  =  AX 
SOLM { J , 11 )  =  BX 
SOLM( J, 12 )  =  AY 
SOLM (J, 13)  =  BY 
SOLM (J, 14)  =  RESULA 
SOLM (J, 15)  =  RESULB 

SOLM{ J, 16)  =  ANGA 
SOLM( J, 17)  =  ANGB 

else 

endif 

enddo 
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if  (FG  .gt.  FGG)  then 
FGG  =  FG 
else 

FGG  =  FGG 
endif 

write (2 , 104)  { (SOL ( III, JJJ) , JJJ=1 , 17) ,  III=1,NFRIC) 

write (3,105)  TIME, F10 , Fll , (FNEW (JJJ) , JJJ  =  1  ,  NFRIC) 
enddo 

do  III  =  1  ,  NFRIC 

write ( 15 , 104)  (SOLM (III , JJJ) , JJJ=1,17) 
enddo 

101  format (2X, ’FRIC.COEF. ' ,9X, 'F10' ,9X, 'Fll' ,9X, 1 FNA' ,9X, ' FNB 1 ,9X, 
l'SPAMOX,  1  SPB '  ,  10X,  •FW'  ,  10X,  'SW',  10X,  'AX'.IOX,  1  BX '  ,  10X,  'AY'  , 
21 OX,  ' BY  1 , 13X,  'RESULTANT'  , 9X,  'ANGLE  WITH  HORIZONTAL') 

102  format (158X, 'A  (11)  B  (10) 1 , 8X, ' A  B  ') 

103  format (185X, ' (+  counterclockwise) 1 ,/) 

104  format ( IX, 100F12 . 3 ) 

105  format ( IX, F12 . 4 , F12 . 4 , F12 . 4 , 10F12 . 4 ) 

106  format (IX, 14 , F12 . 4 , 14 , 18) 

107  format (37X, 10F12 .4) 

108  format ( 2X, 'FRIC.COEF. ' ,9X, 'F10' ,9X, 'Fll' ,9X, 'FNA' ,9X, 'FNB' ,9X, 
1 ' SPA ' ,10X, 'SPB' , 10X, ' FW' ,10X, 'SW' , 10X, 'AX' , 10X, 1 BX 1 , 10X, 'AY' , 
210X, 'BY' , 13X, 'RESULTANT' , 9X, 'ANGLE  WITH  HORIZONTAL') 

109  format (158X, 'A  (11)  B  (10)',8X,'A  B  ') 

110  f ormat ( 185X, ‘ ( +  counterclockwise)1,/) 
stop 

end 


B.2  Program  to  Calculate  Fwfor  Experimental  Data  Using 
the  Scale  Impulse  Force  Formulation 


Program  FWALLMAX 


* 

* 

* 

* 

* 


Program  to  calculate  the  maximum  force  perpendicular  to  the  wall 
due  to  the  impact  of  a  flotilla  of  harges  in  the  experiments . 

The  shear  developed  in  the  wall  due  to  the  impact  is  calculated 
using  the  coefficient  of  friction  between  steel  and  concrete. 


* 
★ 
k 
k 
k  k 


integer *4  ncases,  i,  nvel,  ii,  ntheta,  j,  jj 

real*8  fric,  mass,  vel(100),  theta (100),  dtl  ,  beta,  gama,  frac, 
1  fwm (100,100) 
character*20  namel,  name2 

write ( * , * )  'Enter  input  filename  with  extension' 
read ( * , * )  namel 

write ( * , * )  'Enter  output  filename  with  extension1 
read ( * , * )  name2 
open (UNIT=1 , FILE=namel ) 
open (UNIT=2 , FILE=name2 ) 
read { 1 , * )  ncases 
do  i  =  1  ,  ncases 
read(l,*)  fric 
read(l,*)  nvel 

read(l,*)  (vel(ii),  ii  =  1  ,  nvel) 
read(l,*)  mass 
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read(l,*)  ntheta 

read(l,*)  (theta(ii),  ii  =  1  ,  ntheta) 
readd,*)  dtl  ,  beta,  gama,  frac 
do  j  =  1  ,  nvel 

do  jj  =  1,  ntheta 

fwm(j,jj)  =  (gama  *  frac  *  mass  *  vel(j)  * 

1  sin  (theta  (jj  )  *3d4 159/180.0)  )  /  (dtl  *beta) 


enddo 

enddo 

write(2,100)  i 
write(2,101)  fric 
write (2 , 102 ) 

write(2,103)  (theta(ii),  ii  =  1  ,  ntheta) 
do  j  =  1,  nvel 

write(2,104)  vel ( j ) , ( fwm ( j , j j ) ,  jj  =  1,  ntheta) 
enddo 
enddo 

100  f ormat ( // ,  1  Case  Number  =  ',  16) 

101  f ormat (' Fric .  Coeff.  =',  F12.3) 

102  format (2 8x, 'Theta  Values') 

103  f ormat (llx, 30F12 . 3 ) 

104  format (30F12 . 3 ) 
stop 

end 


B3.  Program  to  Calculate  the  Numerical  Integration  of  the 
Force  Versus  Time  or  Acceleration-Velocity- 
Displacement 

program  numint 

integer*4  option,  np,  i,  nexp,  nnu,  j,  ncc 

real*8  t(50000),  force (50000) ,  rimp,  dt,  acc(20000), 

1  vel(20000),  disp(20000),  degree,  fric(10),  mat (50000, 15) 

character*20  namel,  name2 

*********************************************************************** 

*  Program  to  calculate  the  area  under  the  curve  using  the  * 

*  trapezoidal  rule.  It  has  three  options:  * 

*  * 

*  1.  calculate  the  integral  of  (F  *  dt)  =  I  * 

*  2.  calculate  the  integral  of  (a  *  dt)  =  v  and  v  time  history* 

*  3.  calculate  the  integral  of  (v  *  dt)  =  x  and  x  time  history* 

*********************************************************************** 

write (*,*)  '  Please  select  one  option  to  calculate:1 

write ( * , * ) 

write (*,*)  '  1.  integral  of  (F  *  dt)  =  I' 

write(*,*)  '  2.  integral  of  (a  *  dt)  and  v-t,  x-t  time  history' 
write ( * , * ) 

write (*,*)  '  Your  selection  is?' 

write ( * , * ) 
read(*,*)  option 
if  (option  .eq.  1)  then 

write (*,*)  'Enter  force  input  filename  with  extension' 

read(*,*)  namel 

open (UNIT=1 , FILE=namel ) 
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write (*,*)  ’Enter  force  output  filename  with  extension' 

read(*,*)  name2 

open (UNIT=2 , FILE=name2 ) 

read (I,*)  nexp,  degree,  nnu,  np 

read(l,*)  (fric(i),  i  =  1,  nnu) 

ncc  =  3  +  nnu 

do  i  =  1  ,  np 

readfl,*)  (mat(i,j),  j  =  1  ,  ncc) 
enddo 

do  i  =  1 ,  np 

t (i)  =  mat (i, 1) 
enddo 

dt  =  ( t (np)  -  t(l))  /  (np-1) 
do  j  =  2  ,  3 +nnu 
do  i  =  1  ,  np 

force (i)  =  mat (i, j ) 
enddo 

call  artra (np, force, dt, rimp) 
write ( * , * ) 
write (2, 101)  rimp 
write (* , *) 
enddo 
else 

write (*,*)  'Enter  acceleration  input  filename  with  extension’ 
read(*,*)  namel 

write (*,*)  'Enter  output  filename  with  extension' 

read(*,*)  name2 

open (UNIT=1 , FILE=namel ) 

open (UNIT=2 , FILE=name2 ) 

readd,  *)  np 

read(l,*)  veld) 

read (1 , * )  disp ( 1 ) 

do  i  =  1  ,  np 

readfl,*)  t(i),acc(i) 
enddo 

dt  =  ( t (np)  -  t(l))  /  (np-1) 
do  i  =  2 ,  np 

vel(i)  =  vel(i-l)  +  (acc(i)  +  acc(i-l))  *  dt  *  0.5 
enddo 

do  i  =  2 ,  np 

disp(i)  =disp ( i-1) +vel (i-1) *dt+ ( (2 *acc ( i-1 ) +acc ( i ) ) *dt**2) 
1  /6.0 
enddo 
write  (2 , *) 

write (2,*)  1  Time  Acc.  Vel.  Disp.1 

write (2,*) 
do  i  =  1  ,  np 

write (2, 100)  t(i) ,acc(i) ,vel(i) ,disp(i) 
enddo 
endif 

100  format (4F12 . 4) 

101  format ('The  impulse  due  to  this  force  is  =  ' ,F12.4) 
stop 

end 
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subroutine  artra (np , f , dx, area) 
integer *4  i,  np 
real*8  fx,  f (50000),  dx,  area 
fx  =  0.0 
do  i  =  2  ,  np-1 
fx  =  f (i)  +  fx 
enddo 

area  =dx*  (fx+  (f(l)  +f(np)  )  /2.0) 

return 

end 


B4.  Program  to  Calculate  the  Average,  Standard  Deviation, 
and  Coefficient  of  Variation  of  a  Set  of  Data  Points  and 
the  Slope  of  a  Set  “X”  and  “Y”  Data  Points 

Program  DATA 

*********************************************************************** 

*  Program  to  calculate  the  average,  standard  deviation,  coefficient  * 

*  of  variation  and  the  slope  for  the  linear  regresion  for  a  set  of 

* 

*  data  points.  The  intercept  of  the  resulting  line  is  zero. 
*********************************************************************** 


integer  NPOINTS,  I,  ND 

real*8  DATAX(IOOO),  DATAY(IOOO),  AVGX,  STDX,  COVX,  SLOPE,  COX, 
1COY,  AVGY,  STDY,  COVY,  COXX,  COXY,  COXXX 
character*20  NAMEl 

write (*,*)  'Enter  data  input  filename  with  extension' 

read ( * , * )  NAMEl 

open (UNIT=1 , FILE=NAME1 ) 

* 

*  ND  =  1  FOR  ONE  DATA  COLUMN  "X" 

*  ND  =  2  FOR  TWO  DATA  COLUMNS  "X"  AND  "Y" 

* 

read(l, *)  ND 
read(l,*)  NPOINTS 
do  I  =1,  NPOINTS 

if  (ND  .  eq.  1)  then 
read { 1 , * )  DATAX ( I ) 
else 

read(l,*)  DATAX ( I) ,  DATAY(I) 
endif 
enddo 
COX  =0.0 
do  I  =  1,  NPOINTS 

COX  =  COX  +  DATAX ( I ) 
enddo 

AVGX  =  COX  /  NPOINTS 

COXX  =  0.0 

do  I  =  1,  NPOINTS 

COXX  =  COXX  +  ( DATAX ( I ) -AVGX ) *  *  2 
enddo 

STDX  =  dsqrt (COXX  /  NPOINTS) 

COVX  =  STDX  /  AVGX 
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write ( * , * )  'Average  X  =',AVGX 
write (*,*)  'Standard  Deviation  X  =',  STDX 
write (*,*)  'Coefficient  of  Variation  X  =',  COVX 
write ( * , * ) 
if  (ND  . eq.  2)  then 
COY  =  0.0 
do  I  =  1,  NPOINTS 

COY  =  COY  +  DAT AY ( I ) 
enddo 

AVGY  =  COY  /  NPOINTS 

COY  =0.0 

do  I  =  1,  NPOINTS 

COY  =  COY  +  (DATAY (I) -AVGY) **2 
enddo 

STDY  =  dsqrt (COY  /  NPOINTS) 

COVY  =  STDY  /  AVGY 
write { * , * )  'Average  Y  =',AVGY 
write {*,*)  'Standard  Deviation  Y  =',  STDY 
write (*,*)  'Coefficient  of  Variation  Y  =',  COVY 
else 
endif 

if  (ND  .eq.  2)  then 
COXY  =0.0 
COXXX  =0.0 
do  I  =  1,  NPOINTS 

COXY  =  COXY  +  DATAX(I) * DATAY (I) 

COXXX  =  COXXX  +  DATAX ( I ) *  *  2 
enddo 

slope  =  COXY  /  COXXX 
write (*,*)  'Slope  =’, slope 
else 
endif 
stop 
end 


B.5  Worksheet  to  Calculate  the  Fw  Expression  for  the 
Assumed  Coefficient  of  Friction  Model 

Integrates : 

>  restart: 

>  with  (linalg) : 

Warning,  the  protected  names  norm  and  trace  have  been  redefined  and 
unprotected 

> 

>  Ay : =Spa*cos (b) +Fna*sin (b) ; 

>  Ax : =Fna*cos (b) -Spa* sin (b) ; 

>  By: =Spb*cos (d) +Fnb*sin (d) ; 

>  Bx: =Spb*sin (d) -Fnb*cos (d) ; 
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a= (54+theta) *pi/180 

b= (theta-7) *pi/180 

c=(90*pi/180)+b 

d=90*pi/180-a 

e=90*pi/180-b 

Sw=nu*Fw 

>  EQN1 : =Ay+By-Fw=0 ; 

>  EQN2 : =Bx+Ax-mu*Fw=0 ; 

>  EQN3 : =-Bx* ( sin (e) -sin (a) ) -By* (cos (a) +sin (b) ) +Fw* (sin (b) -mu* (1- 
sin(e) ) )=0; 

>  si : = solve ( {EQN1, EQN2 , EQN3 } , {Fna, Fnb, Fw} ) : 

>  Fw:=eval (Fw, si) ; 


B.6  Worksheet  to  Numerically  Determine  Fw  Using  the 
Energy  Model 

#  Integrales: 

# 

>  restart: 

>  with  (linalg) : 

Warning,  the  protected  names  norm  and  trace  have  been  redefined  and 
unprotected 

> 

VI : =Bx/ (tan (alpha) *sin (alpha) +cos (alpha) ) +By* tan (alpha) / ( tan (alpha) *si 

>  n (alpha) +cos (alpha) ) ; 

>  dVl : =1/ (tan (alpha) *s in (alpha) +cos (alpha) ) ; 

>  N1 : =By/cos (alpha) -Bx* tan (alpha) / ( tan (alpha) * sin (alpha) +cos (alpha) ) - 
By* 

>  (tan(alpha) ) A2/ (tan(alpha) *sin(alpha) +cos (alpha) ) ; 

>  dNl : =- tan (alpha) / (tan (alpha) * sin (alpha) +cos (alpha) ) ; 

>  Ml : =Cl*R*cos (a) *Fw-C2*R*cos (a) *Bx-C3*R*cos (a) *Sw- 
Cl*R*cos (alpha) *Fw+C2 

>  *R*cos (alpha) *Bx+C3*R*cos (alpha) *Sw+Bx*R*sin (alpha) -Bx*R*sin(a) ; 

>  dMl : =-C2 *R*cos (a) +C2 *R*cos (alpha) +R* sin (alpha) -R*sin (a) ; 

>  V2 : = (Fw*tan (beta) +Sw-Bx- 

By*tan(beta) ) / (tan(beta) *sin(beta) +cos (beta) ) ; 

>  dV2 : =-1/ (tan (beta) *sin (beta) +cos (beta) ) ; 

>  N2 : = By /cos (beta) -Fw/cos (beta) + (Fw* ( tan (beta) ) A2+Sw*tan (beta)- 
Bx*tan (be 

>  ta) -By* ( tan (beta) ) A2 ) / ( tan (beta) *sin (beta) +cos (beta) ) ; 
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>  dN2 :=- tan (beta) / (tan (beta) *sin(beta) +cos (beta) ) ; 

>  M2 : =Cl*R*cos (a) *Fw-C2*R*cos (a) *Bx-C3*R*cos (a) *Sw- 
Cl*R*eos (beta) *Fw+C2* 

>  R*cos (beta) *Bx+C3*R*cos (beta) *Sw+Bx*R*sin (beta) - 
Bx*R*sin (a) +Fw*R*cos (b 

>  eta) +R*Sw-R*sin(beta) *Sw; 

>  dM2 : =-C2*R*cos (a) +C2*R*cos (beta) +R* sin (beta) -R*sin{a) ; 

# 

>  INTI :=int (Ml*dMl, alpha=a, .pi/2) ; 

>  INT2 : =int (M2*dM2 , beta=pi/2 .  ,  z )  ; 

>  INT3:=int<Vl*dVl,alpha=a. .pi/2) ; 

>  INT4 : =int ( V2 *dV2 , beta=pi/2 . . z) ; 

>  INT5 : =int (Nl*dNl , alpha=a . .pi/2) ; 

>  INT6 : =int (N2*dN2 , beta=pi/2 . . z) ; 

> 

>  INTI : =subs (cos (pi/2 ) =0 , INTI) ; 

>  INT2 :=subs (cos (pi/2) =0, INT2) ; 

>  INT3 :=subs (cos (pi/2) =0, INT3) ; 

>  INT4:=subs (cos (pi/2) =0, INT4) ; 

>  INT5:=subs (cos (pi/2) =0, INT5) ; 

# 

>  INT6 : =subs (cos (pi/2 ) =0 , INT6) ; 

>  INTI : =subs (sin (pi /2 ) =1, INTI) ; 

>  INT2 ; =subs (sin (pi/2) =1, INT2 ) ; 

>  INT3 :=subs (sin (pi/2) =1, INT3 ) ; 

>  INT4 : =subs (sin (pi/2) =1, INT4) ; 

>  INT5 :=subs (sin (pi/2) =1 , INT5) ; 

>  INT6 : =subs (sin (pi/2) =1, INT6) ; 

>  INTI : =simplify (INTI) ; 

>  INT2 :=simplify (INT2) ; 

>  INT3 : =simplify ( INT3 ) ; 
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>  INT4 : =simplify (INT4) ; 

>  INT5 : =simplify (INT5) ; 

>  INT6 : =simplify (INT6 ) ; 

> 

INT : = ( (R/ (E*Ine) ) * (INT1+INT2) )+{ (R/ (E*Ar ) ) * (INT5+INT6) ) +( (R/ (G*Ar ) ) *(I 

>  NT3+INT4) ) ; 

>  INT: =simplify (INT) ; 

>  Ay:=Spa*cos (b) +Fna*sin(b) ; 

>  Ax:=Fna*cos (b) -Spa*sin(b) ; 

>  By : =Spb*cos (d) +Fnb*sin (d) ; 

>  Bx:=Spb*sin (d) -Fnb*cos (d) ; 

>  thetad: =21 . 5 ; 

>  pi:=3 .141592654; 

>  theta : =thetad*pi/180 ; 

>  Ine : =303 .  75 ; 

>  nu:=0.33; 

>  Ar : =45 ; 

>  E : =29000 ; 

>  G:=E/ (2* ( 1+nu) ) ; 

>  R: =70 . 89 ; 

>  a: = (54+thetad) *pi/180 ; 

>  b : = ( thetad-7 ) *pi/180 ; 

>  c :  =  ( 9 0*pi/180 ) +b; 

>  d:=90*pi/180-a; 

>  e : =90*pi/180-b; 

>  Cl : =sin (b) / (cos (a) +sin (b) ) ; 

>  C2 :  =  (sin (e) -sin (a) ) / (cos (a) +sin (b) ) ; 

>  C3 := (1-sin (e) ) / (cos (a) +sin (b) ) ; 

>  z:=29*pi/180+a; 
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>  Spa: =0.24; 


>  Spb: =0.304; 

>  EQN1:=INT=0; 

>  EQN2 : =Bx+Ax-Sw=0 ; 

>  EQN3 : =Ay+By-Fw=0 ; 

>  EQN4 : =Fw*R*sin (b) -Sw*R* ( 1-sin (e) ) -Bx*R* ( sin ( e) -sin (a) ) - 
By*R* (cos (a) +si 

>  n(b) ) =0 ; 

>  si : =solve ( { EQN1 , EQN2 , EQN3 , EQN4 } , {Fna, Fnb, Fw, Sw} ) ; 

>  si :=eval (si) ; 


B.7  Worksheet  to  Determine  Fw  in  Terms  of  Variable 
Using  the  Energy  Model 

Integrales : 

>  restart: 

>  with  (linalg) : 

Warning,  the  protected  names  norm  and  trace  have  been  redefined  and 
unprotected 

> 

Vl : =Bx/ (tan (alpha) *sin (alpha) +cos (alpha) ) +By* tan (alpha) / (tan (alpha) *sin 
(alpha) +cos (alpha) ) ; 

>  dVl : =1/ (tan (alpha) *sin (alpha) +cos (alpha) ) ; 

>  Nl:=By/cos (alpha) -Bx*tan (alpha) / (tan(alpha) *sin (alpha) +cos (alpha) )- 
By* (tan (alpha) ) A2/ (tan(alpha) *sin(alpha) +cos (alpha) ) ; 

>  dNl :=-tan( alpha) / (tan (alpha) *sin (alpha) +cos (alpha) ) ; 

>  Ml :=Cl*R*cos (a) *Fw-C2*R*cos (a) *Bx-C3*R*cos (a) *Sw- 

Cl*R*cos (alpha) *Fw+C2*R*cos (alpha) *Bx+C3*R*cos (alpha) *Sw+Bx*R*sin (alpha 
) -Bx*R*sin (a) ; 

>  dMl :=-C2*R*cos (a) +C2*R*cos (alpha) +R*sin (alpha) -R*sin(a) ; 

>  V2 : = (Fw*tan (beta) +Sw-Bx- 

By*tan (beta) ) / (tan (beta) *sin(beta) +cos (beta) ) ; 

>  dV2 :=-!/ (tan (beta) *sin (beta) +cos (beta) ) ; 

>  N2 : =By/ cos (beta) -Fw/cos (beta) + (Fw* (tan (beta) ) A2+Sw* tan (beta) - 
Bx*tan (beta) -By* (tan (beta) ) A2) / (tan (beta) *sin (beta) +cos (beta) ) ; 

>  dN2 :=- tan (beta ) / (tan(beta) *sin(beta) +cos (beta) ) ; 

>  M2 : =Gl*R*cos (a) *Fw-C2*R*cos (a) *Bx-C3*R*cos (a) *Sw- 

Cl*R*cos (beta) *Fw+C2*R*cos (beta) *Bx+C3*R*cos (beta) *Sw+Bx*R*sin(beta) - 
Bx*R*sin(a) +Fw*R*cos (beta) +R*Sw-R*sin (beta) *Sw; 

>  dM2 : =-C2 *R*cos (a) +C2*R*cos (beta) +R*sin(beta) -R*sin(a) ; 

>  INTI : =int (Ml*dMl , alpha=a. .pi/2) ; 

>  INT2 : =int (M2*dM2 , beta=pi/2 . . z) ; 

>  INT3 : =int (Vl*dVl , alpha=a . .pi/2) ; 

>  INT4 : =int (V2*dV2 , beta=pi/2 . . z)  ; 

>  INT5 : =int (Nl*dNl, alpha=a. .pi/2) ; 

>  INT6 : =int (N2*dN2 , beta=pi/2 . . z ) ; 


Appendix  B  FORTRAN  Source  Programs  and  Maple  Worksheets 


B13 


>  INTI : =subs (cos (pi/2 ) =0 , INTI) ; 

>  INT2 : =subs (cos (pi/2 ) =0 , INT2 ) ; 

>  INT3 : =subs (cos (pi/2 ) =0 , INT3 ) ; 

>  INT4 : =subs (cos (pi/2 ) =0 , INT4 ) ; 

>  INT5 : =subs (cos (pi/2 ) =0 , INT5) ; 

>  INT6 : =subs (cos (pi/2 ) =0 , INT6) ; 

>  INTI : =subs ( sin (pi/2 ) =1 , INTI) ; 

>  INT2 :=subs (sin (pi/2) =1, INT2) ; 

>  INT3 :=subs (sin (pi/2) =1, INT3) ; 

>  INT4 : =subs (sin (pi/2 ) =1 ,  INT4) ; 

>  INT5 : =subs (sin (pi/2 ) =1 , INT5) ; 

>  INT6 : =subs (sin (pi/2 ) =1 , INT6) ; 

>  INTI : =simplify (INTI) ; 

>  INT2 : =simplify (INT2 ) ; 

>  INT3 : =simplify (INT3 ) ; 

>  INT4 : =simplify (INT4 ) ; 

>  INT5 : =simplify (INT5) ; 

>  INT6 : =simplify (INT6) ; 

> 

INT : = ( (R/ (E*Ine) ) * (INT1+INT2) ) + ( (R/ (E*Ar ) ) * (INT5+INT6) ) +( (R/ (G*Ar) ) * ( IN 
T3+INT4) ) ; 

>  INT : =simplify ( INT) ; 

>  Ay :=Spa*cos (b) +Fna*sin(b) ; 

>  Ax:=Fna*cos (b) -Spa*sin(b) ; 

>  By : =Spb*cos (d) +Fnb*sin (d) ; 

>  Bx:=Spb*sin (d) -Fnb*cos (d) ; 

>  EQNl : =INT=0 ; 

>  EQN2 : =Bx+Ax-Sw=0 ; 

>  EQN3 : =Ay+By-Fw=0 ; 

>  EQN4 : =Fw*R*sin (b) -Sw*R* (1-sin (e) ) -Bx*R* (sin (e) -sin (a) ) - 
By*R* (cos (a) +sin (b) ) =0 ; 

>  si : =solve ( {EQNl , EQN2 , EQN3 , EQN4} , {Fna, Fnb, Fw, Sw) ) : 

>  Fw: =eval (Fw, si ) ; 

>  Sw:=eval (Sw, si) ; 

> 
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Appendix  C 
Numerical  Integration 


One  problem  in  which  numerical  methods  are  used  is  that  of  approximating 
the  area  under  the  graph  of  a  function  y  =  F(x )  from  x  =  atox-b,  thus  obtaining 
an  approximate  value  of  the  integral. 

b 

\f(x)dx 

a 


One  common  method  is  to  divide  the  internal  ( a,b )  into  n  subintervals,  each 
of  length  Ax  =  (b-  a)/n  and  then  form  trapezoids  having  a  base  equal  to  Ax  and 
with  altitudes  given  by  the  values  of  the  function  at  the  beginning  and  at  the  end 
subinterval.  This  is  illustrated  in  Figure  C.l. 


Numerical  Inteqration 


y 


Figure  C.l.  Numerical  integration  scheme 
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The  sum  of  the  areas  of  these  trapezoids  is  given  by 


A  =  'Elf(xi-i)  +  f(xi)]^- 

i=i  £ 

which  can  also  be  written  as 


A  - 


Jm±m 


+  £/(*  l) 


Example: 

Find  the  area  under  the  curve /(x)  =  x2  +  x  +  2  from  x  =  0  to  x  =  4. 


3  2 

A  =  |(x2  +x  +  2)dx  =  —  +  —  +  2x 


-  =  37.33 
0 


Using  num  int.  for  n  =  401  data  point  (PARABOL.DAT) 
A  =  37.2335  produces  a  0.26  percent  error 


C2 
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Example: 


Find  the  velocity  and  displacement  time-history  from  the  next  acceleration  time- 
history,  if 

V(0)  =  0  in./s  and  X  (0)  =  0  in. 


The  exact  solution  is 

a(t)  =  10sin(2^f) 


ft  t  in  1 

V(0-V(0)=  jadt=  jl  0sin(2rt)dt  =  1 0  jsin(2^)r/f  =  — [-cos(2/zf)] 

to= 0  0  0  0 

V(t)  =  [cos(2^f)-l]  +  V(0) 

7t 


x(t)  -  jc(0)  = 

n 


— sm{2rt)-t 
_2  n 


5  sin(2^f) 
n  2  n 
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Displacement  Velocity 


C4 
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Appendix  D 
Statistical  Theory 


D.1  Simple  Statistics 

The  following  information  on  simple  statistics  basically  follows  the 
information  given  in  Part  4  of  Chapra  and  Canale  (1988). 

Data  are  often  given  for  discrete  values  along  a  continuum.  However,  one 
may  require  estimates  at  points  between  the  discrete  values.  This  appendix 
describes  techniques  to  fit  curves  to  such  data  in  order  to  obtain  intermediate 
estimates.  In  addition,  there  may  be  a  need  for  a  simplified  version  of  a 
complicated  function.  One  way  to  do  this  is  to  compute  values  of  the  function  at 
a  number  of  discrete  values  along  the  range  of  interest.  Then,  a  simpler  function 
may  be  derived  to  fit  these  values.  Both  of  these  applications  are  known  as  curve 
fitting. 

There  are  two  general  approaches  for  curve  fitting  that  are  distinguished  from 
each  other  on  the  basis  of  the  amount  of  error  associated  with  the  data.  First, 
where  the  data  exhibit  a  significant  degree  of  error,  the  strategy  is  to  derive  a 
single  curve  that  represents  the  general  trend  of  the  data.  Because  any  individual 
data  point  may  be  incorrect,  one  makes  no  effort  to  intersect  every  point.  Rather, 
the  curve  is  designed  to  follow  the  pattern  of  the  points  taken  as  a  group.  One 
approach  of  this  nature  is  called  least-squares  regression. 

Second,  where  the  data  are  known  to  be  very  precise,  the  basic  approach  is  to 
fit  a  curve  or  a  series  of  curves  that  pass  directly  through  each  of  the  points.  Such 
data  usually  originate  from  tables.  Examples  of  such  data  are  values  for  the 
density  of  water  or  for  the  heat  capacity  of  gases  as  a  function  of  temperature. 

The  estimation  of  values  between  well-known  discrete  points  is  called 
interpolation. 

The  simplest  method  for  fitting  a  curve  to  data  is  to  plot  the  points  and  then 
sketch  a  line  that  visually  conforms  to  the  data.  Although  this  is  a  valid  option 
when  quick  estimates  are  required,  the  results  are  dependent  on  the  subjective 
viewpoint  of  the  person  sketching  the  curve.  This  is  a  very  common  practice  in 
engineering.  If  the  values  are  truly  close  to  being  linear  or  are  spaced  closely, 
such  an  approximation  provides  estimates  that  are  adequate  for  many  engineering 
calculations.  However,  where  the  underlying  relationship  is  highly  curvilinear  or 
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the  data  are  widely  spaced,  significant  errors  can  be  introduced  by  such  linear 
interpolation. 

The  most  common  statistics  term  is  the  arithmetic  mean.  The  arithmetic 
mean  (y)  of  a  sample  is  defined  as  the  sum  of  the  individual  data  points  (yf ) 
divided  by  the  number  of  points  («),  or 


y  = 


(D.l) 


where  the  summation  is  from  i  =  1  through  n.  The  most  common  measure  of 
spread  for  a  sample  is  the  standard  deviation  sy  about  the  mean 


s 


y 


(D.2) 


where  S,  is  the  total  sum  of  the  squares  of  the  residuals  between  the  data  points 

and  the  mean,  or  S  ,-xk-yt  .  Thus,  if  the  individual  measurements  are 

spread  out  widely  around  the  mean,  S,  (and  consequently,  sy)  will  be  large.  If  they 
are  grouped  tightly,  the  standard  deviation  will  be  small.  The  spread  can  also  be 
represented  by  the  square  of  the  standard  deviation,  which  is  called  the  variance. 


s 


2 

y 


S, 

n  —  1 


(D.3) 


Note  that  the  denominator  in  Equations  D.2  and  D.3  is  n  -  1.  This  quantity  is 
referred  to  as  the  degree  of  freedom.  This  nomenclature  derives  from  the  fact  that 
the  sum  of  the  quantities  upon  which  S,  is  based  is  zero.  Another  justification  for 
dividing  by  n  - 1  is  the  fact  that  there  is  no  such  thing  as  the  spread  of  a  single 
data  point. 

A  final  statistic  that  has  utility  in  quantifying  the  spread  of  data  is  the 
coefficient  of  variation  (COV).  This  statistic  is  the  ratio  of  the  standard  deviation 
to  the  mean.  As  such,  it  provides  a  normalized  measure  of  the  spread.  It  is  often 
multiplied  by  100  so  that  it  can  be  expressed  in  the  form  of  a  percent,  as 


5 

COV=-f- 100%  (D.4) 

y 

That  is,  it  is  the  ratio  of  a  measure  of  error  sy  to  an  estimate  of  the  true  value  y . 

Using  the  (Fw)„ulx  values  given  in  Table  5.3  of  the  main  text  (and  reported  in 
Table  6.1a,  eighth  column,  in  yellow),  statistical  values  are  obtained  in  the 
following  calculations  for  the  data  presented  in  Table  D.l . 


D2 
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Table  D.l. 

Linear  Momentum  Normal  to  the  Wall  and  Maximum 

Force  Normal  to  the  Wall  for  Eight  Barge  Impact 
Experiments 

Experiment 

Number 

Linear  Momentum  Normal 
to  the  Wall 

Xl 

{F w)max  Table  5,3 

y> 

29 

897.42 

286.63 

30 

925.73 

B.  ' 

31 

552.52 

37 

649.84 

38 

706.32 

39 

732.74 

41 

812.59 

42 

1025.48 

wsam  •^’r.'TTsn 

Calculation: 


y  = 


2717.42 

8 


=  339.68  kips 


St  - (286.63- 339.68)2  + (369.15 -339.68)2 

1=1 

+  (236.20 -339.68)2  +(327.27 -339.68)2  +  (230.29 -339.68)2 
+  (271.07-339.68)2  +(419.37-339.68)2  +  (577.44 -339.68)2  =94098.72 


5 

y 


94098.72 

V  8-1 


=  115.94 


s 

COV  =-^ 

y 


1 15  94 

- *100  =  34% 

339.68 


D.2  Linear  Regression  and  Quantification  of  Error 

The  Information  presented  in  this  section  (on  linear  regression  and 
quantification  of  error)  basically  follows  the  information  given  in  Chapter  12  of 
Chapra  and  Canale  (1988). 

Where  substantial  error  is  associated  with  data,  polynomial  interpolation  is 
inappropriate  and  may  yield  unsatisfactory  results  when  used  to  predict 
intermediate  values.  Experimental  data  are  often  of  this  type.  A  more  appropriate 
strategy  for  such  cases  is  to  derive  an  approximating  function  that  fits  the  shape 
or  general  trend  of  the  data  without  necessarily  matching  the  individual  points. 
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One  way  to  determine  this  line  is  to  visually  inspect  the  plotted  data  and  then 
sketch  a  “best”  line  through  the  points.  Although  such  “eyeball”  approaches  have 
common-sense  appeal  and  are  valid  for  “back-of-the-envelope”  calculations,  they 
are  deficient  because  they  are  arbitrary.  That  is,  unless  the  points  define  a  perfect 
straight  line  (in  which  case,  interpolation  would  be  appropriate),  different 
analysts  would  draw  different  lines. 

To  remove  this  subjectivity,  some  criterion  must  be  devised  to  establish  a 
basis  for  the  fit.  One  way  to  do  this  is  to  derive  a  curve  that  minimizes  the 
discrepancy  between  the  data  points  and  the  curve.  A  technique  for 
accomplishing  this  objective,  called  least-squares  regression,  will  be  discussed 
here. 

The  simplest  example  of  a  least-squares  approximation  is  fitting  a  straight 
line  to  a  set  of  paired  observations:  (xy ,yt),  feyz),...,  (x„,yn).  The  mathematical 
expression  for  the  straight  line  is 


y  =  a0  +a,x  +  e 


(D.5) 


where  aQ  and  at,  are  coefficients  representing  the  intercept  and  the  slope, 
respectively,  and  e  is  the  error,  or  residual,  between  the  model  and  the 
observations,  which  can  represented  by  rearranging  Equation  D.5  as 


e  =  y-a0  -atx  (D.6) 

Thus,  the  error,  or  residual,  is  the  discrepancy  between  the  true  value  of  y  and  the 
approximate  value,  a0  +  a,x,  predicted  by  the  linear  equation. 

One  strategy  for  fitting  a  “best”  line  through  the  data  would  be  to  minimize 
the  sum  of  the  residual  errors  for  all  the  available  data,  as  in 

=t(y-  ~ao  ~aixi)  (DJ) 

i=l  i=l 

where  n  is  the  total  number  of  points.  However,  this  is  an  inadequate  criterion, 
which  depicts  the  fit  of  a  straight  line  to  two  points.  Obviously,  the  best  fit  is  the 
line  connecting  the  points.  However,  any  straight  line  passing  through  the 
midpoint  of  the  connecting  line  (except  a  perfectly  vertical  line)  results  in  a 
minimum  value  of  Equation  D.7  equal  to  zero  because  the  errors  cancel. 

Another  criterion  would  be  to  minimize  the  sum  of  the  absolute  values  of  the 
discrepancies,  as  in 


D4 
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(D.8) 


n  n 

ZN  =  Zh  ~ao  ~aixi 

i =1  i=l 

This  criterion  also  does  not  yield  a  unique  best  fit.  A  third  strategy  for  fitting 
a  best  line  is  the  minimax  criterion.  In  this  technique,  the  line  is  chosen  that 
minimizes  the  maximum  distance  that  an  individual  point  falls  from  the  line.  This 
strategy  is  ill  suited  for  regression  because  it  gives  undue  influence  to  an  outlier, 
that  is,  a  single  point  with  a  large  error. 

A  strategy  that  overcomes  the  shortcomings  of  the  aforementioned 
approaches  is  to  minimize  the  sum  of  the  squares  of  the  residuals,  Sn  as  in 

n  n  2 

sr  =  Z^2  =  Z (y«  - ao  -  a,x, )  (D.9) 

i=l  i=l 

This  criterion  has  a  number  of  advantages,  including  the  fact  that  it  yields  a 
unique  line  for  a  given  set  of  data.  To  determine  values  for  a0  and 
Equation  D.9  is  differentiated  with  respect  to  each  coefficient.  Setting  these 
derivatives  equal  to  zero  will  result  in  a  minimum  Sr.  Now  the  equations  can  be 
expressed  as  a  set  of  two  simultaneous  linear  equations  with  two  unknowns,  a0 
and  at.  These  are  called  the  normal  equations,  and  they  can  be  solved 
simultaneously,  producing  the  following  equations: 


„  _"Zv,-I>,X>, 

"Xf  (D.10) 

a  =y -a  x 
o  J  1 

where  ybar  ( y  )  is  the  arithmetic  mean  of  data  points  yh  and  xbar  ( 3c )  is  the 
arithmetic  mean  of  data  points  jc,.  Substituting  the  computed  values  for  these 
coefficients  into  Equation  D.5  produces  the  smallest  error  that  can  be  obtained. 

Recalling  that  the  sum  of  the  squares  is  defined  as 

n 

Sr  =Z(>'<  (D.ll) 

»=/ 

where  the  square  of  the  residual  represents  the  square  of  the  vertical  distance 
between  the  data  and  another  measure  of  central  tendency  (the  straight  line),  the 
standard  deviation  for  the  regression  line  can  be  determined  as 


Sx/y  ~ 


(D.12) 
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where  sx/y  is  called  the  standard  error  of  the  estimate.  The  subscript  notation 

x/ y  designates  that  the  error  is  for  a  predicted  value  of  y  corresponding  to  a 
particular  value  of  x.  Also,  notice  that  we  now  divide  by  n-2  because  two  data- 
derived  estimates — a0  and  ai — were  used  to  compute  Sr ;  thus  we  have  lost  two 
degrees  of  freedom.  Another  justification  for  dividing  by  n  -  2  is  that  there  is  no 
such  thing  as  the  “spread  of  data”  around  a  straight  line  connecting  two  points. 
Thus,  for  the  case  where  n  =  2,  Equation  D.12  yields  a  meaningless  result  of 
infinity.  Just  as  was  the  case  with  the  standard  deviation,  the  standard  error  of  the 
estimate  quantifies  the  spread  of  the  data.  However,  sx/y  quantifies  the  spread 

around  the  regression  line  in  contrast  to  the  original  standard  deviation  sy  that 
quantified  the  spread  around  the  mean. 

The  above  concepts  can  be  used  to  quantify  the  “goodness”  of  our  fit.  This  is 
particularly  useful  for  comparison  of  several  regressions.  To  do  this,  we  return  to 
the  original  data  and  determine  the  total  sum  of  the  squares  around  the  mean  for 
the  dependent  variable  (in  our  case,  y).  As  was  presented  earlier,  this  quantity  is 
designated  S,.  This  is  the  uncertainty  associated  with  the  dependent  variable  prior 
to  regression.  After  performing  the  regression,  we  can  compute  Sr,  the  sum  of  the 
squares  of  the  residuals  around  the  regression  line.  This  represents  the 
uncertainty  that  remains  after  the  regression.  It  is,  therefore,  sometimes  called  the 
unexplained  sum  of  the  squares.  The  difference  between  the  two  quantities,  S, 

-  Sr,  quantifies  the  improvement  or  error  reduction  due  to  describing  the  data  in 
terms  of  a  straight  line  rather  than  as  an  average  value.  Because  the  magnitude  of 
this  quantity  is  scale-dependent,  the  difference  is  normalized  to  the  total  error  to 
yield 


r2  _  $ ,  Sr 

S' 


(D.13) 


where  r2  is  called  the  coefficient  of  determination  and  r  is  the  correlation 
coefficient.  For  a  perfect  fit,  Sn  =  0  and  r  -r2  =  1,  signifying  that  the  line 
explains  100  percent  of  the  variability  of  the  data.  For  r  =  r  =  0,  Sr  =  Sh  and  the 
fit  represents  no  improvement. 


There  exists  a  special  case  for  a  straight  line  when  the  resulting  equation 
passes  through  the  origin.  That  is,  the  intercept  of  the  resulting  equation  is  set 
equal  to  zero  prior  to  regression,  the  coefficient  a0  is  set  equal  to  zero,  and  the 
equation  of  the  line  becomes  y  =  ajx.  This  case  is  important  for  the  barge  impact 
data,  as  will  be  discussed  subsequently.  If  the  best-fit  equation  starts  at  the  origin, 
Equation  D.10  transforms  to 


a, 


2X 


for  a0  =0 


(D.14) 
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D.3  Example  1 ;  Linear  Regression  with  Intercept 
Not  Equal  to  Zero-Linear  Momentum  Normal 
to  the  Wall 


Using  the  (/v)max  values  given  in  Table  5.3  (and  reported  in  Table  6.1  a, 
eighth  column,  in  yellow),  a  linear  regression  is  made  to  obtain  the  best-fit 
straight  line,  as  outlined  in  the  following  calculations  for  the  data  (see 
Table  D.l). 

Calculation: 


n  =  8 
s 


£  X.  =6302.63 

;=i 

8 

Y,y.  =2717.42 

i=i  ' 

8 

YJxiyi  -2236350.39 


i=l 

8 


JV  =5137616.654 


i=l 

_  6302.63 
x  - 


8 


=  787.83 


_  2717.42 
y  = - =  339.68 


Cl  / 


”X  *«•?«■  "X^X^  _  8*  {22363 50.39)-  6302.63  *  2717.42 
"X x‘  -(X^)2  8  *  5137616.654  -  (t 5302.63 f 


a0=y-ajx  =  339.68  -  0.5544  *  787.83  =  -97.09 
{Fw  )max  =  atx  +  a0  =0.5544(mv  sin  $)- 97.09 

Sr  =  X(^<  ~ao  ~alXi)2  =  41154.33 


i-l 


hi* 


41154.33 


=  82.82 


ln-2  V  8-2 
2  =  S,  -Sr  _  94098.72  -  41154.33 
S,  ~  94098.72 


0.563 


r  =  ^  r  = 0.75 


0.5544 
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The  resulting  best-fit  straight  line  using  a0  =  -97.09  and  a,  =  0.5544  is  shown 
in  Figure  D.l,  along  with  the  standard  error  of  the  estimate,  sy/x=  82.82. 


Figure  D.l .  Empirical  correlation  between  (Fh,)^,*  and  linear  momentum  normal  to 
the  wall,  straight  line  with  nonzero  intercept 


D.4  Example  2;  Linear  Regression  with  Intercept 
Equal  to  Zero-Linear  Momentum  Normal  to 
the  Wall 

Using  the  (F„,)mat  values  given  in  Table  5.3  (and  reported  in  Table  6.1a, 
eighth  column,  in  yellow),  a  linear  regression  with  the  intercept  set  equal  to  zero 
is  made  to  obtain  the  best-fit  straight  line,  as  outlined  in  the  following 
calculations  for  the  data  (see  Table  D.l). 

Calculation: 


n=  8 
8 

Y,x.  =6302.63 
1=1 
8 

=2717.42 

i=i 

8 

]Tjc  y  =2236350.39 

i=i 

8 

=5137616.654 

i=i 

_  6302.63 


*  = - : - =  787.83 

=  339.68 


8 

2717.42 


8 
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Y,x>y-i  _  2236350.39 


■-  0.435 


^xf  5137616.654 
a0  =0 

( Fw)max  =  aix  +  a0  =  0.43 5(mv  sind) 

n 

Sr  =  ^(y.  -  a,xl  Y  =  43685.76 


i=l 


f  =  °r 
'y/x 

V  n- 


43685.76 

I  8-2 


=  85.33 


2  St-Sr  94098.72  -  43685.76 

r  =— - -  = - =0.535 

S,  94098.72 


r  =  ^rz  =0.732 


The  resulting  best-fit  straight  line  using  a„  =  0  and  at  =  0.435  is  shown  in 
Figure  D.2  along  with  the  standard  error  of  the  estimate,  sy/x  =  85.33. 


Figure  D.2.  Empirical  correlation  between  (F^)^  and  linear  momentum  normal  to 
the  wall,  straight  line  with  zero  intercept 


D.5  Example  3;  Linear  Regression  with  Intercept 
Equal  to  Zero-Kinetic  Energy  Normal  to  the 
Wall 

Using  the  {Fw)max  values  given  in  Table  5.3  (and  reported  in  Table  6.1  a, 
eighth  column,  in  yellow),  a  linear  regression  with  the  intercept  set  equal  to  zero 
is  made  to  obtain  the  best-fit  straight  line,  as  outlined  in  the  following 
calculations  for  the  data  given  in  Table  D.2. 
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Table  D2 

Kinetic  Energy  Normal  to  the  Wall  and  Maximum  Force 
Normal  to  the  Wall  for  Eight  Barge  Impact  Experiments 

Experiment 

Number 

Kinetic  Energy  Normal  to  the 
Wall 

Xl 

{F v)max  Table  5.3 

Y 

29 

214.92 

286.63 

30 

233.20 

369.15 

31 

83.95 

236.20 

37 

114.27 

.327.27 

38 

134.70 

230.29 

39 

141.88 

271.07 

41 

180.59 

419.37 

42 

282.17 

577.44  ; 

Calculation: 

n-  8 
8 


=1385.68 

i=i 

8 

J^y.  =2717.42 

i=i  ‘ 

8 

=  5 1 3060.46 


i=i 

8 


JV  =271180.35 

(=i 

_  1385.68 

;t  = - =  173.21 


_  2717.42 
y=  —  =  339.68 

•'  n 
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al 


Ew  513060.46 
Ex-  ~  271180.35  ~ 


ao  =0 

i.Fw  )mav  =  aix  +  ao=  1.892(o.5mv 2 ) 


Sr=E(y>-aix,Y  =46458.89 


i=J 


Sy!x  ~ 


n  —  2 


43685.76 

8-2 


=  88.00 


. 2  S,-Sr  94098.72  -  46458.89  __ 

'  — - - - - - —  =  0.5063 

S,  94098.72 


0.7115 


The  resulting  best-fit  straight  line  using  a0  =  0  and  a,  =  1 .892  is  shown  in 
Figure  D.2  along  with  the  standard  error  of  the  estimate,  sy/x  =  88.00. 


Figure  D.3.  Empirical  correlation  between  and  kinetic  energy  normal  to  the 

wall,  straight  line  with  zero  intercept 


D.6  Summary  and  Conclusions 

D.6.1  Linear  Momentum  Normal  to  the  Wall 

After  performing  a  linear  regression  of  the  eight  data  points  obtained  from 
Table  5.3  of  the  main  text  (and  repeated  in  Table  6.1a),  a  coefficient  of 
determination  of  0.56  was  obtained  for  the  best-fit  straight  line  (with  nonzero  a0 
and  a,  values).  The  resulting  correlation  coefficient  was  0.75  (Section  D.3). 
These  results  indicate  that  56  percent  of  the  original  uncertainty  has  been 
explained  by  the  linear  model  (Fw  )max  =  0.5544  *  mvsin  9  -  97.09 . 

In  addition  to  this  model,  a  second  regression  was  determined  based  on  the 
constraint  of  having  a  value  equal  to  zero  for  the  intercept  (a0=  0).  In  this  case. 
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the  coefficient  of  determination  was  estimated  as  0.535,  and  the  resulting 
correlation  coefficient  was  0.732  (Section  D.4).  These  results  indicate  that 
53.5  percent  of  the  original  uncertainty  has  been  explained  by  the  linear  model 
(Fw)max  =  0.435*  mv  sine. 

From  a  statistical  point  of  view,  the  authors  concluded  that  both  regression 
lines  are  equally  valid.  There  is  little  difference  between  the  fit  of  both  lines 
through  the  data,  as  reflected  by  the  near  match  of  the  values  for  both  the 
coefficient  of  determination  and  the  correlation  coefficient  for  the  two  lines. 
However,  given  the  complex  nature  of  the  events/interactions  that  transpire  when 
a  barge  flotilla  impacts  a  lock  wall  and  the  gross  simplifications  introduced  by 
this  empirical  correlation,  a  straight  line  with  the  intercept  set  equal  to  zero  is 
recommended  at  this  time.  The  rationale  is  that,  from  an  engineering  point  of 
view,  if  the  velocity  and,  subsequently,  the  linear  momentum  normal  to  the  wall 
are  equal  to  zero,  there  should  be  no  barge  impact  force.  The  authors  would  also 
like  to  note  that  the  amount  of  data  (i.e.,  eight  points)  was  not  sufficient  to  allow 
for  the  fit  of  a  higher  order  polynomial  through  the  data.  Additional  barge  impact 
data,  should  these  become  available  at  a  later  date,  may  indicate  that  a  higher 
order  polynomial  may  be  warranted. 

D.6.2  Kinetic  Energy  Normal  to  the  Wall 

After  performing  a  linear  regression  of  the  eight  data  points  obtained  from 
Table  5.3  (and  repeated  in  Table  6.1a),  a  coefficient  of  determination  of  0.5063 
was  obtained  for  the  best-fit  straight  line  determined  based  on  the  constraint  of 
having  a  value  equal  to  zero  for  the  intercept  ( a0=  0).  The  resulting  correlation 
coefficient  was  0.71 15  (Section  D.5).  These  results  indicate  that  50  percent  of 
the  original  uncertainty  has  been  explained  by  the  linear  model  (Fw  ) 

=  1.892  *  (o.5  mv2). 

From  a  statistical  point  of  view,  the  authors  concluded  that  this  regression 
line  does  not  represent  the  data  points  as  well  as  the  linear  momentum  normal  to 
the  wall  regression  line  with  zero  intercept.  There  is  little  difference  between  the 
coefficients  of  determination  of  both  models.  This  difference  indicates  that  the 
linear  momentum  normal  to  the  wall  approach  (r2  =  0.535)  explains  more  of  the 
original  uncertainty  than  the  kinetic  energy  normal  to  the  wall  approach 
(/  =  0.50).  That  is,  the  closer  the  correlation  coefficient  is  to  unity,  the  better  the 
model. 

Given  the  complex  nature  of  the  events/interactions  that  transpire  when  a 
barge  flotilla  impacts  a  lock  wall  and  the  gross  simplifications  introduced  by 
these  empirical  correlations,  a  straight  line  with  the  intercept  set  equal  to  zero  is 
recommended  at  this  time.  The  rationale  is  that,  from  an  engineering  point  of 
view,  if  the  velocity  and,  subsequently,  the  kinetic  energy  normal  to  the  wall  are 
equal  to  zero,  there  should  be  no  barge  impact  force.  The  authors  would  also  like 
to  note  that  the  amount  of  data  (i.e.,  eight  points)  was  not  sufficient  to  allow  for 
the  fit  of  a  higher  order  polynomial  through  the  data.  Additional  barge  impact 
data,  should  these  become  available  at  a  later  date,  may  indicate  that  a  higher 
order  polynomial  may  be  warranted. 
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